h0=100;
tmax=5;
tmin=0;

[t,y]=ode45('odef',[tmin tmax],[h0;0]);
plot(t,y(:,1),'k')
hold on
[t,y]=ode45('gravf',[tmin tmax],[h0;0]);
plot(t,y(:,1),'r.')
hold off
