This routine computes the elliptic integral F(\phi,k) to the accuracy specified by the mode variable mode. Note that Abramowitz & Stegun define this function in terms of the parameter m = k^2.
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