EllipticF
Incomplete elliptic integral of the first kind
EllipticK
Complete elliptic integral of the first kind
EllipticCK
Complementary complete elliptic integral of the first kind
Calling Sequence
Parameters
Description
Examples
EllipticF(z, k)
EllipticK(k)
EllipticCK(k)
z
-
algebraic expression (the sine of the amplitude)
k
algebraic expression (the modulus)
EllipticF is the Incomplete Elliptic integral of the first kind and is defined by
FunctionAdvisor(definition, EllipticF);
EllipticF⁡z,k=∫0z1−_α12+1⁢−k2⁢_α12+1ⅆ_α1,with no restrictions on ⁡z,k
EllipticK and EllipticCK are respectively the Complete and the Complementary Elliptic integrals of the first kind and are defined by
FunctionAdvisor( definition, EllipticK);
EllipticK⁡k=∫011−_α12+1⁢−k2⁢_α12+1ⅆ_α1,with no restrictions on ⁡k
FunctionAdvisor( definition, EllipticCK);
EllipticCK⁡k=∫011−_α12+1⁢1+k2−1⁢_α12ⅆ_α1,with no restrictions on ⁡k
EllipticK, EllipticCK and EllipticF are related by
FunctionAdvisor( relate, EllipticK,EllipticF);
EllipticK⁡k=EllipticF⁡1,k
FunctionAdvisor( relate, EllipticK,EllipticCK);
EllipticK⁡k=EllipticCK⁡−k2+1
EllipticF is also identical to the InverseJacobiSN function
FunctionAdvisor(relate, EllipticF, InverseJacobiSN);
EllipticF⁡z,k=InverseJacobiSN⁡z,k
and therefore can be used to represent all the InverseJacobiPQ functions provided some restrictions on the function parameters hold.
Elliptic integrals and the related functions are well described in the Table of Integrals Series and Products, Gradshteyn and Ryzhik (G&R) and in the popular Handbook of Mathematical Functions edited by Abramowitz and Stegun (A&S). In A&S, these functions are expressed in terms of a parameter m, representing the square of the modulus k entering the definition of the Elliptic, JacobiPQ and InverseJacobiPQ functions in Maple and G&R. For example, the K⁡m function shown in A&S is numerically equal to the Maple EllipticK⁡m command.
It is worth noting the difference between the Legendre normal form of the Incomplete Elliptic integral of the first kind (see A&S 17.2.7), in Maple represented by EllipticF(z,k) but for the splitting of the square root in the denominator of the integrand (see definition lines above), and the normal trigonometric form of this elliptic integral (see A&S 17.2.6), in Maple represented by the InverseJacobiAM function
InverseJacobiAM(phi,k);
InverseJacobiAM⁡φ,k
(7) = convert((7), Int);
InverseJacobiAM⁡φ,k=∫0φ11−k2⁢sin⁡_θ12ⅆ_θ1
For instance, for -Pi/2 <= phi <= Pi/2 these two forms can be related with ease by changing variables:
EllipticF(z,k);
EllipticF⁡z,k
(9) = convert((9), Int);
EllipticF⁡z,k=∫0z1−_α12+1⁢−k2⁢_α12+1ⅆ_α1
{z=sin(phi), _alpha1=sin(_theta1)}; # -1 <= z <= 1
_α1=sin⁡_θ1,z=sin⁡φ
PDEtools[dchange]((11), (10));
EllipticF⁡sin⁡φ,k=∫0arcsin⁡sin⁡φcos⁡_θ1−sin⁡_θ12+1⁢1−k2⁢sin⁡_θ12ⅆ_θ1
simplify((12)) assuming phi in RealRange(-Pi/2, Pi/2);
EllipticF⁡sin⁡φ,k=∫0φ11−k2⁢sin⁡_θ12ⅆ_θ1
where the right-hand side is actually equal to the trigonometric form InverseJacobiAM⁡φ,k. The general relationship between these two forms and the restriction on the values of the parameters such that the relation is valid are given by
FunctionAdvisor( specialize, InverseJacobiAM, EllipticF);
InverseJacobiAM⁡φ,k=EllipticF⁡sin⁡π⁢12−ℜ⁡φπ+φ,k−2⁢12−ℜ⁡φπ⁢EllipticK⁡k,with no restrictions on ⁡φ,k,InverseJacobiAM⁡φ,k=EllipticF⁡sin⁡φ,k,−π2<ℜ⁡φ∧ℜ⁡φ<π2∨−π2=ℜ⁡φ∧0≤ℑ⁡φ∨π2=ℜ⁡φ∧ℑ⁡φ≤0
FunctionAdvisor( specialize, EllipticF, InverseJacobiAM);
EllipticF⁡z,k=InverseJacobiAM⁡arcsin⁡z,k,with no restrictions on ⁡z,k
Reflection symmetry and special values for EllipticK and EllipticF
FunctionAdvisor⁡special_values,EllipticK
EllipticK⁡−k=EllipticK⁡k,EllipticK⁡0=π2,EllipticK⁡∞=0,EllipticK⁡∞⁢I=0
FunctionAdvisor⁡special_values,EllipticF
EllipticF⁡0,k=0,EllipticF⁡1,k=EllipticK⁡k,EllipticF⁡z,0=arcsin⁡z,EllipticF⁡z,1=arctanh⁡z,EllipticF⁡z,∞=0,EllipticF⁡z,−∞=0
Branch points for EllipticF
FunctionAdvisor⁡branch_points,EllipticF
EllipticF⁡z,k,z∈−1,1,−1k,1k,∞+∞⁢I
Branch points and the branch cut for EllipticK
FunctionAdvisor⁡branch_points,EllipticK
EllipticK⁡k,k∈−∞−∞⁢I,−1,1,∞+∞⁢I
FunctionAdvisor⁡branch_cuts,EllipticK
EllipticK⁡k,k<−1∨1<k
For ℜ⁡k in the cut, so for 1<=Re⁡k<=infinity, EllipticK is continuous from below.
EllipticK⁡2−1100000⁢I
EllipticK⁡2−I100000
=evalf⁡
EllipticK⁡2−I100000=0.8428783289−1.078252932⁢I
EllipticK⁡2
EllipticK⁡2=0.8428751774−1.078257824⁢I
EllipticK⁡2+1100000⁢I
EllipticK⁡2+I100000
EllipticK⁡2+I100000=0.8428783289+1.078252932⁢I
See Also
EllipticCE
EllipticCPi
EllipticE
EllipticPi
FunctionAdvisor
InverseJacobiAM
JacobiAM
RealRange
WeierstrassP
Download Help Document