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.