% Using ode15s to solve the ODE system.
clear all
close all
t = linspace(0,25,1e3); 
%t = linspace(0,20,1e3); 
X0 = [0,1,1.05];
[T,X]=ode15s(@lorenzatt,t,X0);
X0 = [1e-5,1,1.05];
%[T,Y]=ode15s(@lorenzatt,t,X0); % uncomment this if you want a second trajectory

hnd=plot3(X(:,1),X(:,2),X(:,3)); 
set(hnd,'linewidth',2)
%hold;plot3(Y(:,1),Y(:,2),Y(:,3),'r');hold; % for the second trajectory
%set(hnd,'linewidth',1);
xlabel('x');
ylabel('y');
zlabel('z');
axis([-20 20 -30 30 0 50]);
view(45,45)


