[Top][All Lists]

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

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

From: Brian Gough
Subject: Re: [Bug-gsl] Bug in special functions ellint.c
Date: Fri, 28 Nov 2008 14:29:44 +0000
User-agent: Wanderlust/2.14.0 (Africa) Emacs/22.2 Mule/5.0 (SAKAKI)

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

reply via email to

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