%% Main program: Save the program in a separate .m file and run it.
clear all; % clear all variables
t=linspace(0,50,3000)'; % time variables
y0=[-1;3;4]; % Initial conditions
[t,Y] = ode45(@fLorenz,t,y0); %Invoking built-in solver 'ode45'
plot3(Y(:,1),Y(:,2),Y(:,3));  % Plot results
grid on;