function f=lorentz(t,y)
f=zeros(3,1);
%global s r b
f(1)=10*(y(2)- y(1));
f(2)=-y(2)+(30-y(3))*y(1);
f(3)=-(8/3)*y(3)+y(1)*y(2);
