This function computes the incomplete elliptic integral D(\phi,k) which is defined through the Carlson form RD(x,y,z) by the following relation,
D(\phi,k) = (1/3)(\sin(\phi))^3 x RD (1-\sin^2(\phi), 1-k^2 \sin^2(\phi), 1).
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