function f=fmrusa(t,z)
global c q m;

h=forma(z(1));

aux=1+h(2)^2;

f(1)=z(2);
f(2)=-9.8*h(2)/aux-h(2)*h(3)*z(2)^2/aux;
end