This routine computes the incomplete elliptic integral \Pi(\phi,k,n) to the accuracy specified by the mode variable mode. Note that Abramowitz & Stegun define this function in terms of the parameters m = k^2 and \sin^2(\alpha) = k^2, with the change of sign n \to -n.
err contains an estimate of the absolute error in the value z.
The argument mode must be an integer corresponding to
Double-precision, a relative accuracy of approximately 2e-16
.
Single-precision, a relative accuracy of approximately 1e-7
.
Approximate values, a relative accuracy of approximately 5e-4
.
This function is from the GNU Scientific Library, see http://www.gnu.org/software/gsl/ for documentation.
Package: gsl