% mfile for the Lorenz Attractor
% -------------------------------------------------------------------------
clear all, close all
% Initialisation 
% Define System's parameters
a=0.5; b=10; c=8/3; 
% Define ICs
x0=1; y0=1; z0=1;
% -------------------------------------------------------------------------
sim('Lorenz_EEE8013') % Simulate the system
% -------------------------------------------------------------------------
% Plot results
figure
plot(x(:,1),x(:,2)); grid on; xlabel('t, s'), ylabel('x') 
figure
plot(y(:,1),y(:,2)); grid on; xlabel('t, s'), ylabel('y') 
figure
plot(z(:,1),z(:,2)); grid on; xlabel('t, s'), ylabel('z') 
% -------------------------------------------------------------------------