function dy=lorenzeq(t,y)
% Classic values of the parameters, Lorenz Equations

dy(1)=10*(y(2)-y(1));
dy(2)=28*y(1)-y(2)-y(1)*y(3);
dy(3)=y(1)*y(2)-(8/3)*y(3);
dy=dy(:);