
% RHS of Lorenz equations

function F=LorenzEqs(t,vec,param)

 sigma=param(1);
 rho=param(2);
 beta=param(3);

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

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

