At Mon, 25 Aug 2008 17:02:49 +0200,
Thies Heidecke wrote:
> line 489 reads:
> const double tb = -log(y) * (b[0] * y*(b[1] + y*b[2]));
>
> but it should be
> const double tb = -log(y) * (b[0] + y*(b[1] + y*b[2]));
>
> after recompiling the function now gives reasonable values.
>
Thanks for the excellent bug report. I have committed your correction
to the source repository for the next release.
http://git.savannah.gnu.org/gitweb/?p=gsl.git;a=commit;h=163bd7269d3c74af6b0a81ac82d5396a2afd3b3e
--
Brian Gough
GNU Scientific Library -
http://www.gnu.org/software/gsl/
