% Bessel function
clear,close all;
t=0:0.001:10;
z=0:0.01:20;
subplot(211)
J0=besselj(0,z);
J2=besselj(2,z);
J4=besselj(4,z);
J6=besselj(6,z);
J8=besselj(8,z);
plot(z,J0,'r');hold on;
plot(z,J2,'b');hold on;
plot(z,J4,'m');hold on;
plot(z,J6,'c');hold on;
plot(z,J8,'g');hold on;grid on
title('The Bessel function of the first kind')