double comp_ellint_2(double k); float comp_ellint_2f(float k); long double comp_ellint_2l(long double k);
E(k)=E(k,π/2) , for |k|≤1,
where k is k.