cc [ flags ] -I/usr/local/include file(s) -L/usr/local/lib -lmcw [ ... ] #include <mathcw.h> extern float eljkm1f (float u, float k); extern double eljkm1 (double u, double k); extern long double eljkm1l (long double u, long double k); extern __float80 eljkm1w (__float80 u, __float80 k); extern __float128 eljkm1q (__float128 u, __float128 k); extern long_long_double eljkm1ll (long_long_double u, long_long_double k); extern decimal_float eljkm1df (decimal_float u, decimal_float k); extern decimal_double eljkm1d (decimal_double u, decimal_double k); extern decimal_long_double eljkm1dl (decimal_long_double u, decimal_long_double k); extern decimal_long_long_double eljkm1dll (decimal_long_long_double u, decimal_long_long_double k);
NB: Functions with prototypes containing underscores in type names may be available only with certain extended compilers.
The companion function eljq(k) provides the inverse operation, such that
q = eljq(k), k = 1 + eljkm1(q), = eljk(q).
For q > 1/2, we have k > 0.999995, so this function computes an accurate difference of that value from one.
See M. Abramowitz & I. A. Stegun, Handbook of Mathematical Functions, Chapter 16, for definitions of the Jacobian elliptic, theta, and zeta functions.