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

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

figure(101)                 % assign a figure number

plot(xx,yy,'r-','LineWidth',2.0);  % plot the series
grid;                            % provide gridlines

% 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 bessel_gr.eps	   % Create EPS file (Figure)	
end

