Legendre form
In mathematics, the Legendre forms of elliptic integrals are a canonical set of three elliptic integrals to which all others may be reduced. Legendre chose the name elliptic integrals because[1] the second kind gives the arc length of an ellipse of unit semi-major axis and eccentricity [math]\displaystyle{ \scriptstyle{k} }[/math] (the ellipse being defined parametrically by [math]\displaystyle{ \scriptstyle{x = \sqrt{1 - k^{2}} \cos(t)} }[/math], [math]\displaystyle{ \scriptstyle{y = \sin(t)} }[/math]).
In modern times the Legendre forms have largely been supplanted by an alternative canonical set, the Carlson symmetric forms. A more detailed treatment of the Legendre forms is given in the main article on elliptic integrals.
Definition
The incomplete elliptic integral of the first kind is defined as,
- [math]\displaystyle{ F(\phi,k) = \int_0^\phi \frac{1}{\sqrt{1 - k^2 \sin^2(t)}} dt, }[/math]
the second kind as
- [math]\displaystyle{ E(\phi,k) = \int_0^\phi \sqrt{1 - k^2 \sin^2(t)}\,dt, }[/math]
and the third kind as
- [math]\displaystyle{ \Pi(\phi,n,k) = \int_0^\phi \frac{1}{(1 - n \sin^2(t))\sqrt{1 - k^2 \sin^2(t)}}\,dt. }[/math]
The argument n of the third kind of integral is known as the characteristic, which in different notational conventions can appear as either the first, second or third argument of Π and furthermore is sometimes defined with the opposite sign. The argument order shown above is that of Gradshteyn and Ryzhik[2] as well as Numerical Recipes.[3] The choice of sign is that of Abramowitz and Stegun[4] as well as Gradshteyn and Ryzhik,[2] but corresponds to the [math]\displaystyle{ \scriptstyle{\Pi(\phi,-n,k)} }[/math] of Numerical Recipes.[3]
The respective complete elliptic integrals are obtained by setting the amplitude, [math]\displaystyle{ \scriptstyle{\phi} }[/math], the upper limit of the integrals, to [math]\displaystyle{ \scriptstyle{\pi / 2} }[/math].
The Legendre form of an elliptic curve is given by
- [math]\displaystyle{ y^2 = x(x - 1)(x - \lambda) }[/math]
Numerical evaluation
The classic method of evaluation is by means of Landen's transformations. Descending Landen transformation decreases the modulus [math]\displaystyle{ \scriptstyle{k} }[/math] towards zero, while increasing the amplitude [math]\displaystyle{ \scriptstyle{\phi} }[/math]. Conversely, ascending transformation increases the modulus towards unity, while decreasing the amplitude. In either limit of [math]\displaystyle{ \scriptstyle{k} }[/math] approaching zero or one, the integral is readily evaluated.
Most modern authors recommend evaluation in terms of the Carlson symmetric forms, for which there exist efficient, robust and relatively simple algorithms. This approach has been adopted by Boost C++ Libraries, GNU Scientific Library and Numerical Recipes.[3]
References
- ↑ Gratton-Guinness, Ivor (1997). The Fontana History of the Mathematical Sciences. Fontana Press. pp. 308. ISBN 0-00-686179-2.
- ↑ 2.0 2.1 "8.1: Special Functions: Elliptic Integrals and Functions" (in Russian). Tablitsy integralov, summ, rjadov i proizvedenii (5 ed.). Moscow: Nauka. 1971.
- ↑ 3.0 3.1 3.2 William H. Press; Saul A. Teukolsky; William T. Vetterling; Brian P. Flannery (1992). "Chap. 6.11 Special Functions: Elliptic Integrals and Jacobian Functions". Numerical Recipes in C (2 ed.). Cambridge University Press. pp. 261–271. ISBN 0-521-43108-5. https://archive.org/details/numericalrecipes0865unse/page/261.
- ↑ Abramowitz, Milton; Stegun, Irene Ann, eds (1983). "Chapter 17: Elliptic Integrals". Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables. Applied Mathematics Series. 55 (Ninth reprint with additional corrections of tenth original printing with corrections (December 1972); first ed.). Washington D.C.; New York: United States Department of Commerce, National Bureau of Standards; Dover Publications. pp. 589, 589–628. LCCN 65-12253. ISBN 978-0-486-61272-0. http://www.math.sfu.ca/~cbm/aands/page_589.htm.
See also
 |  |
