Loadable Function: [sn, cn, dn] = ellipj (u, m, err)

Compute the Jacobi elliptic functions sn, cn, dn of complex argument and real parameter.

If m is a scalar, the results are the same size as u. If u is a scalar, the results are the same size as m. If u is a column vector and m is a row vector, the results are matrices with length (u) rows and length (m) columns. Otherwise, u and m must conform and the results will be the same size.

The value of u may be complex. The value of m must be 0 <= m <= 1. .

If requested, err contains the following status information and is the same size as the result.

0. Normal return.
1. Error

   no computation, algorithm termination condition not met, return NaN.