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 multithreading
-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:
