** 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 **

** **