Landen's transformation

Last updated October 06, 2024

Landen's transformation is a mapping of the parameters of an elliptic integral, useful for the efficient numerical evaluation of elliptic functions. It was originally due to John Landen and independently rediscovered by Carl Friedrich Gauss.^[1]

Statement

The incomplete elliptic integral of the first kind $F$ is

F(\varphi \setminus \alpha )=F(\varphi ,\sin \alpha )=\int _{0}^{\varphi }{\frac {d\theta }{\sqrt {1-(\sin \theta \sin \alpha )^{2}}}},

where $\alpha$ is the modular angle. Landen's transformation states that if $\alpha _{0}$ , $\alpha _{1}$ , $\varphi _{0}$ , $\varphi _{1}$ are such that $(1+\sin \alpha _{1})(1+\cos \alpha _{0})=2$ and $\tan(\varphi _{1}-\varphi _{0})=\cos \alpha _{0}\tan \varphi _{0}$ , then^[2]

{\begin{aligned}F(\varphi _{0}\setminus \alpha _{0})&=(1+\cos \alpha _{0})^{-1}F(\varphi _{1}\setminus \alpha _{1})\\&={\tfrac {1}{2}}(1+\sin \alpha _{1})F(\varphi _{1}\setminus \alpha _{1}).\end{aligned}}

Landen's transformation can similarly be expressed in terms of the elliptic modulus $k=\sin \alpha$ and its complement $k'=\cos \alpha$ .

Complete elliptic integral

In Gauss's formulation, the value of the integral

I=\int _{0}^{\frac {\pi }{2}}{\frac {1}{\sqrt {a^{2}\cos ^{2}(\theta )+b^{2}\sin ^{2}(\theta )}}}\,d\theta

is unchanged if $a$ and $b$ are replaced by their arithmetic and geometric means respectively, that is

a_{1}={\frac {a+b}{2}},\qquad b_{1}={\sqrt {ab}},

I_{1}=\int _{0}^{\frac {\pi }{2}}{\frac {1}{\sqrt {a_{1}^{2}\cos ^{2}(\theta )+b_{1}^{2}\sin ^{2}(\theta )}}}\,d\theta .

Therefore,

I={\frac {1}{a}}K\left({\frac {\sqrt {a^{2}-b^{2}}}{a}}\right),

I_{1}={\frac {2}{a+b}}K\left({\frac {a-b}{a+b}}\right).

From Landen's transformation we conclude

K\left({\frac {\sqrt {a^{2}-b^{2}}}{a}}\right)={\frac {2a}{a+b}}K\left({\frac {a-b}{a+b}}\right)

and $I_{1}=I$ .

Proof

The transformation may be effected by integration by substitution. It is convenient to first cast the integral in an algebraic form by a substitution of $\theta =\arctan(x/b)$ , $d\theta =(\cos ^{2}(\theta )/b)dx$ giving

I=\int _{0}^{\frac {\pi }{2}}{\frac {1}{\sqrt {a^{2}\cos ^{2}(\theta )+b^{2}\sin ^{2}(\theta )}}}\,d\theta =\int _{0}^{\infty }{\frac {1}{\sqrt {(x^{2}+a^{2})(x^{2}+b^{2})}}}\,dx

A further substitution of $x=t+{\sqrt {t^{2}+ab}}$ gives the desired result

{\begin{aligned}I&=\int _{0}^{\infty }{\frac {1}{\sqrt {(x^{2}+a^{2})(x^{2}+b^{2})}}}\,dx\\&=\int _{-\infty }^{\infty }{\frac {1}{2{\sqrt {\left(t^{2}+\left({\frac {a+b}{2}}\right)^{2}\right)(t^{2}+ab)}}}}\,dt\\&=\int _{0}^{\infty }{\frac {1}{\sqrt {\left(t^{2}+\left({\frac {a+b}{2}}\right)^{2}\right)\left(t^{2}+\left({\sqrt {ab}}\right)^{2}\right)}}}\,dt\end{aligned}}

This latter step is facilitated by writing the radical as

{\sqrt {(x^{2}+a^{2})(x^{2}+b^{2})}}=2x{\sqrt {t^{2}+\left({\frac {a+b}{2}}\right)^{2}}}

and the infinitesimal as

dx={\frac {x}{\sqrt {t^{2}+ab}}}\,dt

so that the factor of $x$ is recognized and cancelled between the two factors.

Arithmetic-geometric mean and Legendre's first integral

If the transformation is iterated a number of times, then the parameters $a$ and $b$ converge very rapidly to a common value, even if they are initially of different orders of magnitude. The limiting value is called the arithmetic-geometric mean of $a$ and $b$ , $\operatorname {AGM} (a,b)$ . In the limit, the integrand becomes a constant, so that integration is trivial

I=\int _{0}^{\frac {\pi }{2}}{\frac {1}{\sqrt {a^{2}\cos ^{2}(\theta )+b^{2}\sin ^{2}(\theta )}}}\,d\theta =\int _{0}^{\frac {\pi }{2}}{\frac {1}{\operatorname {AGM} (a,b)}}\,d\theta ={\frac {\pi }{2\operatorname {AGM} (a,b)}}

The integral may also be recognized as a multiple of Legendre's complete elliptic integral of the first kind. Putting $b^{2}=a^{2}(1-k^{2})$

I={\frac {1}{a}}\int _{0}^{\frac {\pi }{2}}{\frac {1}{\sqrt {1-k^{2}\sin ^{2}(\theta )}}}\,d\theta ={\frac {1}{a}}F\left({\frac {\pi }{2}},k\right)={\frac {1}{a}}K(k)

Hence, for any $a$ , the arithmetic-geometric mean and the complete elliptic integral of the first kind are related by

K(k)={\frac {\pi }{2\operatorname {AGM} (1,{\sqrt {1-k^{2}}})}}

By performing an inverse transformation (reverse arithmetic-geometric mean iteration), that is

a_{-1}=a+{\sqrt {a^{2}-b^{2}}}\,

b_{-1}=a-{\sqrt {a^{2}-b^{2}}}\,

\operatorname {AGM} (a,b)=\operatorname {AGM} \left(a+{\sqrt {a^{2}-b^{2}}},a-{\sqrt {a^{2}-b^{2}}}\right)\,

the relationship may be written as

K(k)={\frac {\pi }{2\operatorname {AGM} (1+k,1-k)}}\,

which may be solved for the AGM of a pair of arbitrary arguments;

\operatorname {AGM} (u,v)={\frac {\pi (u+v)}{4K\left({\frac {u-v}{v+u}}\right)}}.

Related Research Articles

In mathematics, the arithmetic–geometric mean of two positive real numbers $x$ and $y$ is the mutual limit of a sequence of arithmetic means and a sequence of geometric means. The arithmetic–geometric mean is used in fast algorithms for exponential, trigonometric functions, and other special functions, as well as some mathematical constants, in particular, computing $π$ .

In integral calculus, an elliptic integral is one of a number of related functions defined as the value of certain integrals, which were first studied by Giulio Fagnano and Leonhard Euler. Their name originates from their originally arising in connection with the problem of finding the arc length of an ellipse.

In mechanics and geometry, the 3D rotation group, often denoted SO(3), is the group of all rotations about the origin of three-dimensional Euclidean space $under the operation of composition.$

In mathematics and physical science, spherical harmonics are special functions defined on the surface of a sphere. They are often employed in solving partial differential equations in many scientific fields. The table of spherical harmonics contains a list of common spherical harmonics.

In mathematics, the inverse trigonometric functions are the inverse functions of the trigonometric functions, under suitably restricted domains. Specifically, they are the inverses of the sine, cosine, tangent, cotangent, secant, and cosecant functions, and are used to obtain an angle from any of the angle's trigonometric ratios. Inverse trigonometric functions are widely used in engineering, navigation, physics, and geometry.

In mathematics, the Clausen function, introduced by Thomas Clausen, is a transcendental, special function of a single variable. It can variously be expressed in the form of a definite integral, a trigonometric series, and various other forms. It is intimately connected with the polylogarithm, inverse tangent integral, polygamma function, Riemann zeta function, Dirichlet eta function, and Dirichlet beta function.

In probability theory, the Borel–Kolmogorov paradox is a paradox relating to conditional probability with respect to an event of probability zero. It is named after Émile Borel and Andrey Kolmogorov.

In mathematics, theta functions are special functions of several complex variables. They show up in many topics, including Abelian varieties, moduli spaces, quadratic forms, and solitons. As Grassmann algebras, they appear in quantum field theory.

In mathematics, the Jacobi elliptic functions are a set of basic elliptic functions. They are found in the description of the motion of a pendulum, as well as in the design of electronic elliptic filters. While trigonometric functions are defined with reference to a circle, the Jacobi elliptic functions are a generalization which refer to other conic sections, the ellipse in particular. The relation to trigonometric functions is contained in the notation, for example, by the matching notation $for . The Jacobi elliptic functions are used more often in practical problems than the Weierstrass elliptic functions as they do not require notions of complex analysis to be defined and/or understood. They were introduced by Carl Gustav Jakob Jacobi. Carl Friedrich Gauss had already studied special Jacobi elliptic functions in 1797, the lemniscate elliptic functions in particular, but his work was published much later.$

In trigonometry, tangent half-angle formulas relate the tangent of half of an angle to trigonometric functions of the entire angle.

In mathematics, the Legendre chi function is a special function whose Taylor series is also a Dirichlet series, given by

In Bayesian statistics, the Jeffreys prior is a non-informative prior distribution for a parameter space. Named after Sir Harold Jeffreys, its density function is proportional to the square root of the determinant of the Fisher information matrix:

In calculus, the Leibniz integral rule for differentiation under the integral sign, named after Gottfried Wilhelm Leibniz, states that for an integral of the form $where and the integrands are functions dependent on the derivative of this integral is expressible as where the partial derivative indicates that inside the integral, only the variation of with is considered in taking the derivative.$

<span class="mw-page-title-main">Multiple integral</span> Generalization of definite integrals to functions of multiple variables

In mathematics (specifically multivariable calculus), a multiple integral is a definite integral of a function of several real variables, for instance, $f (x, y)$ or $f (x, y, z)$ .

There are several equivalent ways for defining trigonometric functions, and the proofs of the trigonometric identities between them depend on the chosen definition. The oldest and most elementary definitions are based on the geometry of right triangles and the ratio between their sides. The proofs given in this article use these definitions, and thus apply to non-negative angles not greater than a right angle. For greater and negative angles, see Trigonometric functions.

A ratio distribution is a probability distribution constructed as the distribution of the ratio of random variables having two other known distributions. Given two random variables X and Y, the distribution of the random variable Z that is formed as the ratio Z = X/Y is a ratio distribution.

In mathematics, vector spherical harmonics (VSH) are an extension of the scalar spherical harmonics for use with vector fields. The components of the VSH are complex-valued functions expressed in the spherical coordinate basis vectors.

In calculus, the integral of the secant function can be evaluated using a variety of methods and there are multiple ways of expressing the antiderivative, all of which can be shown to be equivalent via trigonometric identities,

References

↑ Gauss, C. F.; Nachlass (1876). "Arithmetisch geometrisches Mittel, Werke, Bd. 3". Königlichen Gesell. Wiss., Göttingen: 361–403.
↑ Abramowitz, Milton; Stegun, Irene Ann, eds. (1983) [June 1964]. Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables. Applied Mathematics Series. Vol. 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. ISBN 978-0-486-61272-0. LCCN 64-60036. MR 0167642. LCCN 65-12253.

Louis V. King On The Direct Numerical Calculation Of Elliptic Functions And Integrals (Cambridge University Press, 1924)

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[1] Gauss, C. F.; Nachlass (1876). "Arithmetisch geometrisches Mittel, Werke, Bd. 3". Königlichen Gesell. Wiss., Göttingen: 361–403.

[2] Abramowitz, Milton; Stegun, Irene Ann, eds. (1983) [June 1964]. Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables. Applied Mathematics Series. Vol. 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. ISBN 978-0-486-61272-0. LCCN 64-60036. MR 0167642. LCCN 65-12253.

[1]

[2]