freefall.m

function [v_analytical,v_terminal,t]=freefall(h,timespan)
    % help file for freefall.m
    % N is number of timesteps between 0 and 12 sec
    % v_an...
    t=[0:h:timespan]';
    c=0.25; m=60; g=9.81; v_terminal=sqrt(m*g/c);

    v_analytical = v_terminal*tanh(g*t/v_terminal);
    v_numerical=zeros(length(t),1);
    delta_time =diff(t);
    for i=1:length(t)-1
        v_numerical(i+1)=v_numerical(i)+(g-c/m*v_numerical(i)^2)*delta_time(i);
    end
    plot(t,v_numerical,'o-')
end