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
|
|
EllipticF(z, k)
EllipticK(k)
EllipticCK(k)
|
|
Parameters
|
|
z
|
-
|
algebraic expression (the sine of the amplitude)
|
k
|
-
|
algebraic expression (the parameter)
|
|
|
|
|
Description
|
|
|
EllipticF is the Incomplete Elliptic integral of the first kind and is defined by
|
>
|
FunctionAdvisor(definition, EllipticF);
|
| (1) |
|
EllipticK and EllipticCK are respectively the Complete and the Complementary Elliptic integrals of the first kind and are defined by
|
>
|
FunctionAdvisor( definition, EllipticK);
|
| (2) |
>
|
FunctionAdvisor( definition, EllipticCK);
|
| (3) |
|
EllipticK, EllipticCK and EllipticF are related by
|
>
|
FunctionAdvisor( relate, EllipticK,EllipticF);
|
| (4) |
>
|
FunctionAdvisor( relate, EllipticK,EllipticCK);
|
| (5) |
>
|
FunctionAdvisor(relate, EllipticF, InverseJacobiSN);
|
| (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 function shown in A&S is numerically equal to the Maple 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);
|
| (7) |
>
|
(7) = convert((7), Int);
|
| (8) |
|
For instance, for -Pi/2 <= phi <= Pi/2 these two forms can be related with ease by changing variables:
|
| (9) |
>
|
(9) = convert((9), Int);
|
| (10) |
>
|
{z=sin(phi), _alpha1=sin(_theta1)}; # -1 <= z <= 1
|
| (11) |
>
|
PDEtools[dchange]((11), (10));
|
| (12) |
>
|
simplify((12)) assuming phi in RealRange(-Pi/2, Pi/2);
|
| (13) |
|
where the right-hand side is actually equal to the trigonometric form . 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);
|
| (14) |
>
|
FunctionAdvisor( specialize, EllipticF, InverseJacobiAM);
|
| (15) |
|
|
Examples
|
|
Reflection symmetry and special values for EllipticK and EllipticF
>
|
|
| (16) |
>
|
|
| (17) |
Branch points for EllipticF
>
|
|
| (18) |
Branch points and the branch cut for EllipticK
>
|
|
| (19) |
>
|
|
| (20) |
For in the cut, so for , EllipticK is continuous from below.
>
|
|
| (21) |
>
|
|
| (22) |
>
|
|
| (23) |
>
|
|
| (24) |
>
|
|
| (25) |
>
|
|
| (26) |
>
|
|
|
|