[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[Bug-gsl] Bug in special functions ellint.c

From: Sascha Merz
Subject: [Bug-gsl] Bug in special functions ellint.c
Date: Fri, 28 Nov 2008 07:50:11 +1100
User-agent: Thunderbird (X11/20060911)

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 Merz
PhD Student
School of Mechanical and Manufacturing Engineering
University of New South Wales, SYDNEY NSW 2052 AUSTRALIA
Ph: +61 (2) 9385 6864

reply via email to

[Prev in Thread] Current Thread [Next in Thread]