A program for numerical simulation of a whole ski race, from start to finish, is developed in MATLAB. The track is modelled by a set of cubical splines in two dimensions and can be used to simulate a track in a closed loop or with the start and finish at different locations. The forces considered in the simulations are gravitational force, normal force between snow and skis, drag force from the wind, frictional force between snow and ski and driving force from the skier. The differential equations of motion are solved from start to finish with the Runge-Kutta method. Different wind situations during the race can be modelled, as well as different glide conditions on different parts of the track. It is also possible to vary the available power during the race. The simulation program's output is the total time of the race, together with the forces and speed during different parts of the race and intermediate times at selected points. Some preliminary simulations are also presented.