
function F=Rossler(t,vec)
 a=0.2;
 b=0.2;
 c=5.0;
 x=vec(1);
 y=vec(2);
 z=vec(3);
 F(1)= -y-z;
 F(2)= x+a*y;
 F(3)= b+z*(x-c);
 F=transpose(F);

