double ellint_3(double k, double nu, double phi); float ellint_3f(float k, float nu, float phi); long double ellint_3l(long double k, long double nu, long double phi);
Π(ν,k,ϕ)=∫ϕ0dθ(1−νsin2θ)√1−k2sin2θ , for |k|≤1,
where ν is nu, k is k, and φ is phi.