function f=ode2(t,y)
f1=(1-y(2).^2).*y(1)-y(2);
f2=y(1);
f=[f1;f2];

