function f=torod(t,u)

par = 2;
f(1) = (par-3)*u(1)-0.25*u(2)+u(1)*(u(3)+0.2*(1-u(3)*u(3)));
f(2) = 0.25*u(1)+(par-3)*u(2)+u(2)*(u(3)+0.2*(1-u(3)*u(3)));
f(3) = par*u(3)-(u(1)*u(1)+u(2)*u(2)+u(3)*u(3));
f=f';
%y0=[0.1 0.1 0.1]; [t,y]=ode15s('torod',[0 200],y0); plot3(y(:,1),y(:,2),y(:,3));
%plot(y(:,1),y(:,2)); xlabel('x_1'); ylabel('x_2');
