global AA BB RR
AA = 10; BB = 8/3; RR = 28;
[t, u] = ode45('lor', [0 7], [1 2 3]);

plot(t, u, 'LineWidth', 2); grid
title('A solution to Lorenz system')
xlabel('t')
%gtext('x')
%gtext('y')
%gtext('z')
fprintf('Press <Enter> to continue\n')
pause

clf;
[t, u] = ode45('lor', [0 20], [1 2 3], odeset('OutputFcn','odeplot'));
pause

plot3(u(:,1), u(:,2), u(:,3))
grid
