EllipticCK - Maple Help
For the best experience, we recommend viewing online help using Google Chrome or Mozilla Firefox.

Online Help

All Products    Maple    MapleSim


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

Calling Sequence

EllipticF(z, k)

EllipticK(k)

EllipticCK(k)

Parameters

z

-

algebraic expression (the sine of the amplitude)

k

-

algebraic expression (the modulus)

Description

  

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

(1)
  

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

(2)
> 

FunctionAdvisor( definition, EllipticCK);

EllipticCK⁡k=∫011−_α12+1⁢1+k2−1⁢_α12ⅆ_α1,with no restrictions on ⁡k

(3)
  

EllipticK, EllipticCK and EllipticF are related by

> 

FunctionAdvisor( relate, EllipticK,EllipticF);

EllipticK⁡k=EllipticF⁡1,k

(4)
> 

FunctionAdvisor( relate, EllipticK,EllipticCK);

EllipticK⁡k=EllipticCK⁡−k2+1

(5)
  

EllipticF is also identical to the InverseJacobiSN function

> 

FunctionAdvisor(relate, EllipticF, InverseJacobiSN);

EllipticF⁡z,k=InverseJacobiSN⁡z,k

(6)
  

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)
> 

(7) = convert((7), Int);

InverseJacobiAM⁡φ,k=∫0φ11−k2⁢sin⁡_θ12ⅆ_θ1

(8)
  

For instance, for -Pi/2 <= phi <= Pi/2 these two forms can be related with ease by changing variables:

> 

EllipticF(z,k);

EllipticF⁡z&comma;k

(9)
> 

(9) = convert((9), Int);

EllipticF⁡z&comma;k=∫0z1−_&alpha;12+1⁢−k2⁢_&alpha;12+1&DifferentialD;_&alpha;1

(10)
> 

{z=sin(phi), _alpha1=sin(_theta1)};     #  -1 <= z <= 1

_&alpha;1=sin⁡_&theta;1&comma;z=sin⁡φ

(11)
> 

PDEtools[dchange]((11), (10));

EllipticF⁡sin⁡φ&comma;k=∫0arcsin⁡sin⁡φcos⁡_&theta;1−sin⁡_&theta;12+1⁢1−k2⁢sin⁡_&theta;12&DifferentialD;_&theta;1

(12)
> 

simplify((12)) assuming phi in RealRange(-Pi/2, Pi/2);

EllipticF⁡sin⁡φ&comma;k=∫0φ11−k2⁢sin⁡_&theta;12&DifferentialD;_&theta;1

(13)
  

where the right-hand side is actually equal to the trigonometric form InverseJacobiAM⁡φ&comma;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⁡φ&comma;k=EllipticF⁡sin⁡π⁢12−ℜ⁡φπ+φ&comma;k−2⁢12−ℜ⁡φπ⁢EllipticK⁡k&comma;with no restrictions on ⁡φ&comma;k,InverseJacobiAM⁡φ&comma;k=EllipticF⁡sin⁡φ&comma;k&comma;−π2<ℜ⁡φ∧ℜ⁡φ<π2∨−π2=ℜ⁡φ∧0≤ℑ⁡φ∨π2=ℜ⁡φ∧ℑ⁡φ≤0

(14)
> 

FunctionAdvisor( specialize, EllipticF, InverseJacobiAM);

EllipticF⁡z&comma;k=InverseJacobiAM⁡arcsin⁡z&comma;k&comma;with no restrictions on ⁡z&comma;k

(15)

Examples

Reflection symmetry and special values for EllipticK and EllipticF

> 

FunctionAdvisor⁡special_values&comma;EllipticK

EllipticK⁡−k=EllipticK⁡k&comma;EllipticK⁡0=π2&comma;EllipticK⁡∞=0&comma;EllipticK⁡∞⁢I=0

(16)
> 

FunctionAdvisor⁡special_values&comma;EllipticF

EllipticF⁡0&comma;k=0&comma;EllipticF⁡1&comma;k=EllipticK⁡k&comma;EllipticF⁡z&comma;0=arcsin⁡z&comma;EllipticF⁡z&comma;1=arctanh⁡z&comma;EllipticF⁡z&comma;∞=0&comma;EllipticF⁡z&comma;−∞=0

(17)

Branch points for EllipticF

> 

FunctionAdvisor⁡branch_points&comma;EllipticF

EllipticF⁡z&comma;k&comma;z∈−1&comma;1&comma;−1k&comma;1k&comma;∞+∞⁢I

(18)

Branch points and the branch cut for EllipticK

> 

FunctionAdvisor⁡branch_points&comma;EllipticK

EllipticK⁡k&comma;k∈−∞−∞⁢I&comma;−1&comma;1&comma;∞+∞⁢I

(19)
> 

FunctionAdvisor⁡branch_cuts&comma;EllipticK

EllipticK⁡k&comma;k<−1∨1<k

(20)

For ℜ⁡k in the cut, so for 1<=Re⁡k<=infinity, EllipticK is continuous from below.

> 

EllipticK⁡2−1100000⁢I

EllipticK⁡2−I100000

(21)
> 

=evalf⁡

EllipticK⁡2−I100000=0.8428783289−1.078252932⁢I

(22)
> 

EllipticK⁡2

EllipticK⁡2

(23)
> 

=evalf⁡

EllipticK⁡2=0.8428751774−1.078257824⁢I

(24)
> 

EllipticK⁡2+1100000⁢I

EllipticK⁡2+I100000

(25)
> 

=evalf⁡

EllipticK⁡2+I100000=0.8428783289+1.078252932⁢I

(26)

See Also

EllipticCE

EllipticCPi

EllipticE

EllipticPi

FunctionAdvisor

InverseJacobiAM

JacobiAM

RealRange

WeierstrassP