double ellint_1(double k, double phi); float ellint_1f(float k, float phi); long double ellint_1l(long double k, long double phi);
F(k,ϕ)=∫ϕ0dθ√1−k2sin2θ , for |k|≤1,
where k is k and φ is phi.