At Sat, 29 Nov 2008 08:48:06 +1100, Sascha Merz wrote: > I still did not get the right result and found the bug. Line 489 is > const double tb = -log(y) * (b[0] * y*(b[1] + y*b[2])); > however, according to Abramowitz and Stegun, equation 17.3.34, it should be > const double tb = -log(y) * (b[0] + y*(b[1] + y*b[2])); > i. e. replace the second * by a +.
Thanks, this was recently fixed in the repository and will be in the next release: http://savannah.gnu.org/bugs/?24146 -- Brian Gough Support freedom by joining the FSF http://www.fsf.org/news/fall-2008-fundraiser _______________________________________________ Bug-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/bug-gsl
