global g c;
g=9.814;
c=0.3;
r0=[0;10.;10.;0.];
tt=0:0.01:3 ;
options = odeset('RelTol',1e-8,'AbsTol',1e-8,'MaxStep',1.);
[t,r]=ode45(@motion,tt,r0,options);

plot(r(:,1),r(:,2))
