double ellint_2(double k, double phi); float ellint_2f(float k, float phi); long double ellint_2l(long double k, long double phi);
E(k,ϕ)=∫ϕ0√1−k2sin2θdθ , for |k|≤1,
where k is k and φ is phi.