— Function: double
gsl_sf_ellint_Kcomp (
double k, gsl_mode_t mode)
— Function: int
gsl_sf_ellint_Kcomp_e (
double k, gsl_mode_t mode, gsl_sf_result * result)
These routines compute the complete elliptic integral K(k) to the
accuracy specified by the mode variable mode.