[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
http://www.fsf.org/news/fall-2008-fundraiser