Re: [Bug-gsl] Re: Auto-response for your message to address@hidden

[Brian Gough]

address@hidden writes:
 > #include <gsl/gsl_sf_bessel.h>
 > 
 > int main(int argc, char *argv[]) {
 > 
 >   int lmax = 94;
 >   double x = 27246.7832109820;
 >   double result_array[95];
 > 
 >   gsl_sf_bessel_jl_array(lmax, x, result_array);
 > }

Thanks, this has been added to bug number 3 "gsl_sf_bessel_Jn fails
for large argument" in the GSL BUGS file.  It is caused by a breakdown
in the continued fraction approximation for large arguments.

If you find a better algorithm please let us know.

best regards,

-- 
Brian Gough

Network Theory Ltd            
15 Royal Park                 
Bristol BS8 3AL               
United Kingdom                

Tel: +44 (0)117 3179309 
Fax: +44 (0)117 9048108              
Web: http://www.network-theory.co.uk/