clc, clear all;
H0=2000;
t0=0;
v0=0;
v=v0;
h=H0;
t=t0;
T=t;
H=h;
while h>0
   U=[h;v];
   U=U+rksteg(t,U,0.1, @dhojdDt);
   h=U(1);
   v=U(2);
   t=t+0.1;
   H=[H h];
   T=[T t];
   
end
plot(T,H)