This script uses Halley's method to compute the positive k zeros of the Bessel function of the first kind J(n,x) and second kind Y(n,x) where n is a positive number. The routine has been tested for up to k=100 and n=100.

Mostly seems to work. However, the zeros for the Bessel functions of the second kind are missing the first zero (0.893576966279). That is, the first Y_n zero it gives is actually the second, the second is actually the third, and so on.

I am using the function in Matlab R2008a
I entered the command besselzero(1/2,1,2)
and the result is different than the value from
Mathematica s BesselYZero[1/2,1]
I also tried some other combinations which seemed fine.
I just wonder if this is a special case or I should check
the numbers given by this Matlab function with some other online sources. Thanks

08 Oct 2008

Jesús Aurelio

i dont have comments jeje sorry but this informaion is good ! ! thank you!!!

02 Mar 2008

Yaniv Shaposhnik

Thank you, very useful

29 Jan 2008

Farshid Ghassemi

The algorithm is fascinating and the speed is satisfactory. But there is a small defect: For Y function, the initial guess is not correctly chosen and the routine misses the first root in some cases.

02 Oct 2007

miguel m

05 May 2007

ionel aurel

01 May 2007

Michele Campisi

It works just fine

23 Apr 2007

harrouni lyes

j'ai male a trouver la solution de l'equation dde la gaine fini et infini de la fibre vouler vous m'aider svp.

20 Mar 2007

Warren Parsons

Nice work. You've saved me a lot of time.

26 Nov 2006

Bryan Romero

works great

15 Nov 2006

R Parthasarathy

Good -- more commenting / discussion would be nice.

13 Oct 2006

Kathirvel Thiyagarajan

I am using it. But not sure about its accuracy.

18 Aug 2006

Jim Morehead

11 Jul 2006

T S

Very good. Works perfectly.
Thanks
T.S
BGU University Israel

20 Apr 2006

Hung Lin

05 Dec 2005

PinceDor Ema

French student thank you.
We hope you'll have a good life.

03 Dec 2005

Joey Pittman

I used this m-file to generate the zeros for a bessel function of the first kind and order zero and it worked just fine.

When you run the function, you should use besselzero(n,k,kind) I interpreted the documentation to suggest that the function needs only 2 inputs, but it really requires 3 to work.

04 Aug 2005

Stepan Matejka

It seems the besselzero(n,k,1) and besselzero(n,k,2) works fine. Checked with multiple plot(besselzero(n,100,kind)).
Good work ...