Re: [Bug-gsl] bug in gsl_sf_synchrotron_1

[Brian Gough]

John Houck writes:
 > On Mon, Jan 16, 2006 at 14:31 +0000, Brian Gough wrote:
 > > John Houck writes:
 > >  > Here are some computed values of the function and its first
 > >  > differences in the neighborhood of this point:
 > >  > 
 > >  >     x                F(x)           delta-F(x)
 > >  > 4.212000000e-08  7.478926692e-03  -5.918277062e-07
 > >  > 4.213000000e-08  7.479518519e-03  -5.917340625e-07
 > >  > 4.214000000e-08  7.480110253e-03  -5.151888128e-07  <---
 > >  > 4.215000000e-08  7.480625442e-03  -5.915287484e-07
 > >  > 4.216000000e-08  7.481216971e-03  -5.914352158e-07

Hello,

The problem was that the taylor series needed an extra term to match
up to the chebyshev expansion at the end of its range there.  I've
added an extra term to the taylor series which removes the problem.
This is now in CVS and will be in the next release.

Thanks for the bug report.

-- 
Brian Gough
(GSL Maintainer)

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