Loadable Function: z = gsl_sf_ellint_D (arg1, arg2, mode)
Loadable Function: [z, err] = gsl_sf_ellint_D (…)

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

0 = GSL_PREC_DOUBLE

Double-precision, a relative accuracy of approximately 2e-16.

1 = GSL_PREC_SINGLE

Single-precision, a relative accuracy of approximately 1e-7.

2 = GSL_PREC_APPROX

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