function f=ejem1(t,y)
f=-9.8+0.1*y^2;