Loadable Function: z = gsl_sf_ellint_RJ (arg1, arg2, arg3, arg4, mode)
Loadable Function: [z, err] = gsl_sf_ellint_RJ (…)

This routine computes the incomplete elliptic integral RJ(x,y,z,p) to the accuracy specified by the mode variable mode.

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