  function LorenzCheck()
% Checks out LorenzZones
sigma = 10; rho = 28; beta = 8/3; eta = sqrt(beta*(rho-1));
u = [rho-1;eta;eta];           v = [rho-1;-eta;-eta];
% Starting point and total duration...
y0 = v + [0;0;3];
tFinal = 50;
% The simulation...
[Enter1,Exit1,Enter2,Exit2] = LorenzZones(y0,tFinal);
clc
fprintf('  Enter1   Exit1  Enter2  Exit2\n--------------------------------------\n')
L1 = length(Enter1);
L2 = length(Exit1);
L3 = length(Enter2);
L4 = length(Exit2);
for k=1:max([L1 L2 L3 L4])
    if k<=L1
        fprintf('  %6.3f',Enter1(k))
    else
        fprintf('        ')
    end
        if k<=L2
        fprintf('  %6.3f',Exit1(k))
    else
        fprintf('        ')
    end
    if k<=L3
        fprintf('  %6.3f',Enter2(k))
    else
        fprintf('        ')
    end
    if k<=L4
        fprintf('  %6.3f',Exit2(k))
    else
        fprintf('        ')
    end
    fprintf('\n')
end
% comet3(yOut(:,1),yOut(:,2),yOut(:,3))
%shg

% P2: 2 points
% LorenzZones: 3 points
