Jacobi's elliptic functions sn(u|m), cn(u|m), dn(i|m) and φ(u|m) are defined by using the incomplete elliptic integral of the first kind:
JacobianEllipticFunctions subroutine calculates these functions using arithmetic-geometric mean algorithm.
This article is licensed for personal use only.
ALGLIB Project offers you two editions of ALGLIB:
ALGLIB Free Edition:
delivered for free
offers full set of numerical functionality
extensive algorithmic optimizations
no low level optimizations
non-commercial license
ALGLIB Commercial Edition:
flexible pricing
offers full set of numerical functionality
extensive algorithmic optimizations
high performance (SMP, SIMD)
commercial license with support plan
Links to download sections for Free and Commercial editions can be found below: