function xdot= tennis0p(t,x)
global g alpha
v= sqrt(x(3)^2+x(4)^2);
xdot = [ x(3)
         x(4)
         -alpha*0.508*x(3)*v
         -g-alpha*0.508*x(4)*v ];

