function ydot=lorenz(t,y,a)
ydot(1)=10*(y(2)-y(1));
ydot(2)=28*y(1)-y(2)-y(1)*y(3);
ydot(3)=y(1)*y(2)-a/3*y(3)');

