clear				% Clear previous definitions
figure(1)           % Assign figure number
clf                 % Clear previous figures
hold off            % Start with fresh graph

mytitle = 'Boundary Layer Problem'; 	% Title
xlab = '$x$'; 		                % X-label
ylab = '$y$'; 	                    % Y-label
xx = linspace(0,1,500);
ep = 0.05;
y1 = (exp(-xx) - exp(-xx/ep))./(exp(-1) - exp(-1/ep));
y2 = exp(1-xx);
y3 = exp(1) - exp(1-xx/ep);
plot(xx,y1,'r-','LineWidth',1.5);  
hold on                           
plot(xx,y2,'r--','LineWidth',1.5);
plot(xx,y3,'m--','LineWidth',1.5);
grid                                % Adds Gridlines
h =legend('Exact', 'Outer Approx',...
    'Inner Approx','Location','southeast');
set(h,'Interpreter','latex')
%text(6,-1.25,'$\varepsilon = 0.2$','FontSize',14,...
%    'interpreter','latex');
xlim([0 1]);
ylim([0 3]);
fontlabs = 'Times New Roman';   % Font in labels
xlabel(xlab,'FontSize',14,'FontName',fontlabs,...
    'interpreter','latex');  
ylabel(ylab,'FontSize',14,'FontName',fontlabs,...
    'interpreter','latex');  
title(mytitle,'FontSize',16,'FontName',...
    'Times New Roman','interpreter','latex'); 
set(gca,'FontSize',12);         % Axis tick font size
print -depsc bvplayio_plot.eps		% Create figure as EPS file

