Jcd Jcd_sub Jcn Jcn_sub Jcs Jcs_sub Jdc Jdc_sub Jdn Jdn_sub Jds Jds_sub Jelli Jnc Jnc_sub Jnd Jnd_sub Jns Jns_sub Jsc Jsc_sub Jsd Jsd_sub Jsn Jsn_sub Jtheta1 Jtheta1_sub Jtheta1t Jtheta1t_sub Jtheta2 Jtheta2_sub Jtheta2t Jtheta2t_sub Jtheta3 Jtheta3_sub Jtheta3m1 Jtheta3m1_sub Jtheta3m1t Jtheta3m1t_sub Jtheta3t Jtheta3t_sub Jtheta4 Jtheta4_sub Jtheta4m1 Jtheta4m1_sub Jtheta4m1t Jtheta4m1t_sub Jtheta4t Jtheta4t_sub

Jelli() Sub Jelli ( U As  Double, K As  Double, Sn As  Double, Cn As  Double, Dn As  Double, Optional Info As  Long  ) Jacobi elliptic functions sn(u, k), cn(u,k), dn(u, k) Purpose This routine computes the Jacobi elliptic functions sn(u, k), cn(u,k) and dn(u, k). The functions are defined as follows, given: u = ∫ 1/√(1 - k^2(sin(t))^2) dt [0, φ] The angle φ is called the amplitude and: sn(u, k) = sinφ cn(u, k) = cosφ dn(u, k) = √(1 - k^2(sin(φ))^2) Parameters [in] U Argument u. [in] K Argument k. [out] Sn sn(u, k) value. [out] Cn cn(u, k) value. [out] Dn dn(u, k) value. [out] Info (Optional) = 0: Successful exit. = 1: Floating point range error. Reference boost/math/special_functions