function v=lorenz(t,x)
v=[10*(x(2)-x(1)),28*x(1)-x(2)-x(1)*x(3),-(8/3)*x(3)+x(1)*x(2)];