function long_jump=f(t,y,p)

%This function will attempt to find the trajectory of a long jumper using
%the differential equations provided in the online textbook problem #7.16

g=9.8;      %acceleration due to gravity
c=0.72;     %drag coefficient
m=80;       %mass of long jumper in kg
s=0.5;      %effective surface area of jumper

long_jump=[-g*cos(y(1))/y(2);-(c*p*s/(2*m))*(y(2))^2*((cos(y(1)))^2+sin(y(1))^2)-g*sin(y(1));y(2)*cos(y(1));y(2)*sin(y(1))];