function f=f(time,y)
% note f and y are in general vectors
f=[y(2);y(1)-y(1)^3];
