This routine computes the incomplete elliptic integral RD(x,y,z) 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
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