function f=fEXS711(t,U,k)
global g m;
f=[U(2);
    (k*U(2)^2-m*g)/m];