function f=dy3(t,x,a)
f=zeros(2,1);
f(1)=x(1)+x(2)-x(1)^3;
f(2)=-a*x(1);