function y=lorenz(t,x)
y(1)=10*(x(2)-x(1));
y(2)=28*x(1)-x(2)-x(1)*x(3);
y(3)=x(1)*x(2)-8*x(3)/3;
end
