MATLAB: Zeros of Bessel functions

bessel functons

Has anyone posted an m-file which computes the zeros of bessel functions, J_n and Y_n?

Why not just use FZERO?

 zroot = fzero(@(z) besselj(1, z), 3)
 besselj(1, zroot) % Should be small

I think you need to use element-wise division with the z:

a=besselj(0, z) - besselj(1, z)./z;

Download this package from FEX: https://www.mathworks.com/matlabcentral/fileexchange/48403-bessel-zero-solver and then run the following code. You will need to specify the number of roots you want

n = 2;
p = 10;
R3 = 90*10e-3;
num_roots = 10;
bessel_z = besselzero(n*p, num_roots);
alpha_k = bessel_z/R3;