%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%    Draw a picture of the Slow Manifold.
%
%                SLOWMAN
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Y = linspace(-0.5,1,151);
Z = linspace(0,1,101);

YY = 0*Y;

plot(Y,YY,'LineWidth',1);  % Plot Y-axis
hold on
plot([0,0],[0,1],'LineWidth',1);  % Plot Z-axis

SS = 0.5*Y.^2;
plot(Y,SS,'LineWidth',2)   %  Plot Slow Manifold

plot([0.65,0.65],[0,0.55]);  % Plot vertical line


textY=text(1,0,'Y');
set(textY,'FontSize',22,'FontWeight','bold','FontAngle','italic');
textZ=text(0,1,'\bf Z');
set(textZ,'FontSize',22,'FontWeight','bold','FontAngle','italic');
textS=text(1,0.5,'\bf S');
set(textS,'FontSize',22,'FontWeight','bold','FontAngle','italic');

plot(0.65,0,'O'); plot(0.65,0,'+'); plot(0.65,0,'x');
textL=text(0.65+0.02,0-0.03,'L');
set(textL,'FontSize',16,'FontWeight','bold');
plot(0.65,0.5*(0.65)^2,'O');
plot(0.65,0.5*(0.65)^2,'+');
plot(0.65,0.5*(0.65)^2,'x');
textN=text(0.65+0.02,0.5*(0.65)^2-0.02,'N');
set(textN,'FontSize',16,'FontWeight','bold');
plot(0.65,0.55,'O'); plot(0.65,0.55,'+'); plot(0.65,0.55,'x');
textA=text(0.65+0.02,0.55-0.02,'A');
set(textA,'FontSize',16,'FontWeight','bold');

axis('off')

print -deps slowman.eps
hold off

return