function Lorenz_63
%#######################################################################
%   LORENZ_63.M
%
%   DESCRIPTION
%     Lorenz_63 solves the Lorenz equations with y(t0)=y0 and 
%     x(t0)=0=z(t0) from t0=0 to tf.
%
%     w(t) = [x(t), y(t), z(t)], and the Lorenz equations are 
%       dw/dt = F(w).
%
%   AUTHOR
%     Dr. Carl Gardner, carl.gardner@asu.edu
%     Taylor Hines, taylor.hines@asu.edu
%
%#######################################################################
%
%   INITIAL PARAMETERS
%
y0 = 
tf = 
tspan = [0 tf];
w0 = [0; y0; 0];

%#######################################################################
%
%   SOLUTION
%
options = odeset('RelTol',1e-3,'AbsTol',1e-3);
%[t,w] = ode45( )

%#######################################################################
%
%   FIGURES
%
figure
%Projection on the X-Y Plane.
plot(w(:,1),w(:,2),0,0,'r.')
xlabel('X','FontSize',16)
ylabel('Y','FontSize',16)
title('X-Y Plane','FontSize',16)

figure
%Projection on the Y-Z Plane.

figure
%Plot of Y vs. T.

%#######################################################################
%
%   THE MODEL
%
    function wprime = F(t,w)
    wprime = [
        %Define the system here
             ];
    end

end