scipy.special.ellipj(x1, x2[, out1, out2, out3, out4])
- (sn,cn,dn,ph)=ellipj(u,m) calculates the Jacobian elliptic functions of
parameter m between 0 and 1, and real u. The returned functions are
often written sn(u|m), cn(u|m), and dn(u|m). The value of ph is such
that if u = ellik(ph,m), then sn(u|m) = sin(ph) and cn(u|m) = cos(ph).