Legendre Polynomials

• See Legendre polynomials/Catalogs for the first 12 Legendre polynomials.

In mathematics, the Legendre polynomials P_n(x) are orthogonal polynomials in the variable −1 ≤ x ≤ 1. Their orthogonality is with unit weight,

${\displaystyle {\displaystyle \underset{-1}{\overset{1}{\int }}}{P}_n(x)P_{n^{\prime }}(x)dx=0\text{for}n\ne n^{\prime }.}$

In physics they commonly appear as a function of a polar angle 0 ≤ θ ≤ π with x = cosθ

${\displaystyle {\displaystyle \underset{0}{\overset{\pi }{\int }}}{P}_n(\cos \theta )P_{n^{\prime }}(\cos \theta )\sin \theta d\theta =0\text{for}n\ne n^{\prime }.}$.

The polynomials as function of cosθ are part of the solution of the Laplace equation in spherical polar coordinates.

By the sequential Gram-Schmidt orthogonalization procedure applied to {1, x, x², x³, …} the n^th degree polynomial P_n can be constructed recursively. The Gram-Schmidt procedure applies to all members of the family of orthogonal polynomials, such as Hermite polynomials, Chebyshev polynomials, etc. Further, P_n(x) has in common with the other orthogonal polynomials that it has exactly n real distinct zeroes. These zeroes are used as grid points in Gauss quadrature (numerical integration) schemes.

Historical note[edit]

The polynomials were named "Legendre coefficients" by the British mathematician Isaac Todhunter^[1] in honor of the French mathematician Adrien-Marie Legendre (1752–1833), who was the first to introduce and study them. Todhunter called the functions "coefficients", instead of "polynomials", because they appear as coefficients of hⁿ in the expansion of the generating function (see below); Todhunter also introduced the notation P_n, which is still generally used.

Legendre's polynomials have been introduced by Legendre in a memoir Sur l'attraction des sphéroïdes homogènes published in the Mémoires de Mathématiques et de Physique, présentés à l'Académie royale des sciences par sçavants étrangers, Tome x, pp. 411–435, Paris, 1785.^[2]. The functions also occur in a memoir by Laplace received by the Academy in 1782, Théorie des attractions des sphéroïdes et de la figure des planètes, but the original introduction of the functions appears nevertheless to be due to Legendre, whose work was not published for several years after it was written; his memoir is mentioned as approved for publication in the report of the sittings of the Academy for 1783. Legendre himself declares that Laplace introduced the potential (i.e., generating) function, but that he himself developed the expansion. Later, the mathematicians Jacobi, Dirichlet, and Heine have written on the question of priority and agree that Legendre deserves the credit.

Rodrigues' formula[edit]

The French amateur mathematician Olinde Rodrigues (1795–1851) proved the following formula

${\displaystyle P_n(x)=\frac{1}{2^{n}n!}\frac{d^n(x^2-1{)}^n}{d{x}^n}.}$

Using the Newton binomial and the equation

${\displaystyle \frac{d^n{x}^m}{d{x}^n}=\frac{m!}{(m-n)!}{x}^{m-n},\text{for}n\le m,}$

we get the explicit expression

${\displaystyle P_n(x)=\frac{1}{2^{n}n!}{\displaystyle \sum _{k=\lceil n/2\rceil }^n}(-1{)}^{n-k}(\genfrac{}{}{0ex}{}{n}{k})\frac{(2k)!}{(2k-n)!}{x}^{2k-n}.}$

Substitution p=n-k gives this formula a slightly different appearance

${\displaystyle P_n(x)=\frac{1}{2^n}{\displaystyle \sum _{p=0}^{\lfloor n/2\rfloor }}(-1{)}^p\frac{(2n-2p)!}{p!(n-p)!(n-2p)!}{x}^{n-2p}.}$

Generating function[edit]

The coefficients of hⁿ in the following expansion of the generating function are Legendre polynomials

${\displaystyle \frac{1}{\sqrt{1-2xh+h^2}}=\sum _{n=0}{P}_n(x)h^n.}$

The expansion converges for |h| < 1. This expansion is useful in expanding the inverse distance between two points r and R

${\displaystyle \frac{1}{|\mathbf{𝐫}-\mathbf{𝐑}|}=\frac{1}{\sqrt{r^2+R^2-2rR\cos \gamma }}=\frac{1}{R}\frac{1}{\sqrt{h^2+1-2hx}},}$

where

${\displaystyle h\equiv \frac{r}{R}\text{and}x\equiv \cos \gamma \equiv \mathbf{𝐫}\cdot \mathbf{𝐑}/(rR).}$

Obviously the expansion makes sense only if R > r. The function appears in Newton's gravitational potential and in Coulomb's electrostatic potential.

Normalization[edit]

The polynomials are not normalized to unity, but

${\displaystyle {\displaystyle \underset{-1}{\overset{1}{\int }}}{P}_n(x)P_m(x)dx=\frac{2}{2n+1}{\delta }_{nm},}$

where δ_nm is the Kronecker delta.

Differential equation[edit]

The Legendre polynomials are solutions of the Legendre differential equation

${\displaystyle \frac{d}{dx}[(1-x^2)\frac{d}{dx}{P}_n(x)]+n(n+1)P_n(x)=(1-x^2)\frac{d^2{P}_n(x)}{d{x}^2}-2x\frac{d{P}_n(x)}{dx}+n(n+1)P_n(x)=0.}$

This differential equation has another class of solutions: Legendre functions of the second kind Q_n(x), which are infinite series in 1/x. These functions are of lesser importance.

Note that the differential equation has the form of an eigenvalue equation with eigenvalue −n(n+1) of the operator

${\displaystyle \frac{d}{d\cos \theta }{\sin }^2\theta \frac{d}{d\cos \theta }=\frac{1}{\sin \theta }\frac{d}{d\theta }\sin \theta \frac{d}{d\theta }}$

This operator is the θ-dependent part of the Laplace operator ∇² in spherical polar coordinates.

Properties of Legendre polynomials[edit]

Legendre polynomials have parity (−1)ⁿ under x → -x,

${\displaystyle P_n(-x)=(-1{)}^n{P}_n(x).}$

The following condition normalizes the polynomials

${\displaystyle P_n(1)=1.}$

Recurrence Relations[edit]

Legendre polynomials satisfy the recurrence relations

${\displaystyle (1-x^2)\frac{d}{dx}{P}_n=(n+1)x{P}_n-(n+1)P_{n+1}=n{P}_{n-1}-nx{P}_n}$

${\displaystyle (n+1)P_{n+1}-(2n+1)x{P}_n+n{P}_{n-1}=0}$

From these two relations follows easily

${\displaystyle \frac{d(P_{n+1}-P_{n-1})}{dx}=(2n+1)P_n.}$

External link[edit]

Weisstein, Eric W. "Legendre Polynomial." From MathWorld--A Wolfram Web Resource. [1]

Footnote[edit]