function yprime = Jump(t,y,g,m,c,s,rho)

yprime = [y(4)*cos(y(3));
    y(4)*sin(y(3));
    -g/y(4)*cos(y(3));
    -c*rho*s/2*(y(4))^2/m-g*sin(y(3))];


