function f=f711(t,u)
global k
m=75;
g=9.81;
f=[u(2) (k/m)*u(2)^2-g]';