clear; clc;
v0 = 0;
t0 = 0;
tF = 20;                % end time
timespan = [ t0 tF ];   % only two elements, first and last
[t_soln, v_soln] = ode45('parachuteODE', timespan, v0);
plot( t_soln, v_soln,'.' ); 

