% Bessel functions J_0(z) and Y_0(z)

z = linspace(0,20,500);

j0 = besselj(0,z);
y0 = bessely(0,z);

plot(z,j0,'b-','LineWidth',1.5);
hold on
plot(z,y0,'r-','LineWidth',1.5);
grid;

text(2,0.75,'$J_0(z)$','FontSize',18,'Color','b','Interpreter','latex'); 
text(1,-0.75,'$Y_0(z)$','FontSize',18,'Color','r','Interpreter','latex');

xlim([0,20]);
ylim([-1,1]);
xlabel('$z$','FontSize',16,'FontName',fontlabs, ...
    'interpreter','latex');  
ylabel('','FontSize',16,'FontName',fontlabs, ...
    'interpreter','latex'); 
axis on; 
set(gca,'FontSize',16);

print -depsc bess_jy0.eps
