<aside> 📌 Below is the report I wrote for the final project of the numerical methods class I took fall of my Sophomore year. The project was to use numerical methods in MATLAB to model and solve a real-world physics problem


ezgif.com-gif-maker (1).gif

Name: Jesse Gilbert

Due: December 21, 2021

Class: Numerical Methods (ES-55)

Final Project


This program simulates a skier skiing down a slope with varied terrain and stops when the skier reaches the end. The program uses ODEs to solve for the skier’s position and velocity over time, taking into account factors such as air resistance, friction with the snow, and the varied downward force the skier experiences with different hill curvatures.

An animation is generated to show the skier as they go down the hill. The ski slope is made up of a straight line that has a slope of -1 and a user-defined height, and then a parabola with a user-defined shape.


  1. Factor of area of skier (2 doubles the area, .5 halves the area) (default=1)
  2. Factor of coefficient of kinetic friction of snow (2 doubles it simulating slushing conditions for example, .5 halves the it simulating icy conditions for example) (default=1)
  3. Two attributes of the ski hill
    1. the starting height of the hill (default=5)
    2. a shape factor of the parabolic bottom part of the hill (default=1)


The program does the following:

Numerical Methods

  1. ODEs (using ode45, with events)
  2. Twice differentiation of ski slope (using diff)
  3. Interpolation (using spline)
    1. Of derivatives of ski slope within ode45 evaluation
    2. Of position of skier for animation
Powered by Fruition