function delta =dhojdDt(t,U)
m=75;
g=9.81;
if t<30
    k=m*g/2500;
else
    k=m*g/25;
end
delta=[U(2); k*(U(2))^2/m-g];