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 EPSMO1RU2.

  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