function f=f(time,y)
f=[y(2);y(1)-y(1)^3];
