function dw = Lorentz(t,w,sigma,rho,beta)

dw(1) = sigma*(w(2) - w(1));
dw(2) = w(1)*(rho - w(3)) - w(2);
dw(3) = w(1)*w(2) - beta*w(3);

end