Dear Sir/Madam, 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.
Best regards, Sascha -- Sascha Merz PhD Student School of Mechanical and Manufacturing Engineering University of New South Wales, SYDNEY NSW 2052 AUSTRALIA Ph: +61 (2) 9385 6864 _______________________________________________ Bug-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/bug-gsl
