function bessel2_plot(Lx,res,tol)
% Create a plot using the bessel function

xx = linspace(0,Lx,res); % x points of evaluation
yy = bessel(xx,0,tol);   % run bessel series program  

figure(102)                 % assign a figure number

plot(xx,yy,'k-','LineWidth',2.0);  % plot the series
hold on
grid;      % provide gridlines
for m = 1:3
    if (m == 1)
        lsty = 'b-';
    elseif (m == 2)
        lsty = 'r-';
    else
        lsty = 'm-';
    end
    yy = bessel(xx,m,tol);
    plot(xx,yy,lsty,'LineWidth',2.0);
end
h = legend('$J_0(x)$', '$J_1(x)$', '$J_2(x)$', '$J_3(x)$',...
    'Location','northeast');
 set(h,'Interpreter','latex')
    
% Set up fonts and labels for the Graph
fontlabs = 'Times New Roman';  
xlabel('$x$','FontSize',16,'FontName',fontlabs, ...
    'interpreter','latex');  
ylabel('$J_m(x)$','FontSize',16,'FontName',fontlabs, ...
    'interpreter','latex'); 
mytitle = 'Bessel Plot';
title(mytitle,'FontSize',16,'FontName', ...
    'Times New Roman','interpreter','latex'); 
set(gca,'FontSize',16);        

print -depsc bessel2_gr.eps	   % Create EPS file (Figure)	
end

