function dxdt=osc(t,x)
r2=1-x(1)^2-x(2)^2;
dxdt=[x(2)+x(1)*r2;-x(1)+x(2)*r2];
