At Fri, 28 Nov 2008 07:50:11 +1100, Sascha Merz wrote: > there is a bug in the file ellint.c, lines 486, 487, function > gsl_sf_ellint_Kcomp_e. The elliptic integral is not computed correctly > near 1.0, since the coefficients in the arrays a and b are wrong. They > are not as given in Abramowitz and Stegun, equation 17.3.33. Also > computation of ta and tb is not as in Abramowitz and Stegun.
Thanks for the email, the values are from 17.3.34 which is a more accurate version of 17.3.33 -- I will update the comment. -- Brian Gough (GSL Maintainer) 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
