

% To invoke do: [T,Y]=threeplot('lorenzmio',0.1,0.1,0.1,0,30);

function F=lorenzmio(t,vec)

 sigma=10.0;
 rho=28.0;
 beta=2.66666666;

 x=vec(1);
 y=vec(2);
 z=vec(3);

 F(1)= sigma* (y -x);
 F(2)= rho* x - y - x*z;
 F(3)= -beta* z + x*y;

 F=transpose(F);

