% 4-5-2008
% Lorentz Attractor example

clear all;
close all;

% Parameter values for our ODE
sigma = 10;
beta = 8/3;
rho = 28;

% ODE solver parameters
y0 = [0 1 0];     % Initial conditions
tspan = [-10 10];  % time interval we would like to solve for

% Solve our ode using 4th Order Runge-Kutta
[t,y] = ode45('lorentz_rhs',tspan,y0,[],sigma,rho,beta);

plot3(y(:,1),y(:,2),y(:,3));
box on;
xlabel('x'); ylabel('y'); zlabel('z');