Jacobi elliptic functions SN, CN and DN.
[Sn,Cn,Dn] = ELLIPJ(U,M) returns the values of the Jacobi
elliptic functions SN, CN and DN, evaluated at argument U
and parameter M. As currently implemented, M is limited
to 0 < M < 1.
ELLIPJ(U,M) are accurate to EPS .
U and M must either be matrices of the same size or either
one may be a scalar.
Be sure you don't confuse the modulus K with the parameter M -
they are related in the following way: M = K^2