Re: [Bug-gsl] bug in gsl_sf_bessel_jl

[Brian Gough]

Mario Santos writes:
 > There is a possible big problem with the spherical bessel routine:
 > gsl_sf_bessel_jl(L, x) When running with L<1000 I can see a second
 > peak at x=1001.  For instance, try L=100 and plot the output from
 > x=90 to x=1100...

Thanks for the bug report, I've found the problem at L=100,x=1000. It
is now fixed in CVS and will be in the next release.

-- 
Brian Gough
(GSL Maintainer)

Network Theory Ltd,
Commercial support for GSL --- http://www.network-theory.co.uk/gsl/