function dx = my_model(t,x)

dx = zeros(2,1);  %column vector
dx(1) = x(2);
dx(2) = cos(x(1))-0.5*x(2)-x(1)-0.02*x(1)^3;



