double comp_ellint_1(double k); float comp_ellint_1f(float k); long double comp_ellint_1l(long double k);
K(k)=F(k,π/2) , for |k|≤1,
where k is k.