s=10;
r=28;
b=8/3;
[T,Y]=ode45(@lorentz,[0 50],[3; 2; 7]);
%plot(T, Y(:,1));
plot3(Y(:,1),Y(:,2),Y(:,3))