function dydt = lorenzD(t,y,s,b,r)
dydt = [s*(y(2)-y(1));r*y(1)-y(2)-y(1)*y(3);-b*y(3)+y(1)*y(2)];