Π(n,φ,k)=∫0φdθ(1−nsin2θ)1−k2sin2θ