% Lorenz equations
function dxdt = lorenz(x,t)
   % Parameters for the Lorenz system
   global global_beta global_sigma global_rho
      
   % ODEs
   dxdt(1) = global_sigma*( x(2) - x(1) ); 
   dxdt(2) = global_rho*x(1) - x(1)*x(3) - x(2);
   dxdt(3) = x(1)*x(2) - global_beta*x(3);



