function dxdt = lorenz(t,x)
sig=10.0;
rho=28;
bet=8/3;
dxdt = [sig*(x(2)-x(1));x(1)*(rho-x(3))-x(2);x(1)*x(2)-bet*x(3)];


