F(φ,k)=sinφRF(cos2φ,1−k2sin2φ,1)E(φ,k)=sinφRF(cos2φ,1−k2sin2φ,1)−13k2sin3φRD(cos2φ,1−k2sin2φ,1)Π(n,φ,k)=sinφRF(cos2φ,1−k2sin2φ,1)+13nsin3φRJ(cos2φ,1−k2sin2φ,1,1−nsin2φ)