We are creating the world's most trusted encyclopedia and knowledge base.
Once you join us and log in, you'll be able to edit this page instantly!

Legendre polynomials

From Citizendium, the Citizens' Compendium

Jump to: navigation, search

Image:Statusbar1.png
Main Article
Talk
Definition [?]
Related Articles  [?]
Bibliography  [?]
External Links  [?]
 
This is a draft article, under development. These unapproved articles are subject to a disclaimer.

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


\int_{-1}^{1} P_{n}(x) P_{n'}(x) dx = 0\quad \hbox{for}\quad n\ne n'.

The polynomials are named after the French mathematician Legendre (1752–1833).

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

 
\int^{\pi}_{0} P_{n}(\cos\theta) P_{n'}(\cos\theta) \sin\theta \;d\theta  = 0\quad \hbox{for}\quad n\ne n'.
.

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 nth degree polynomial Pn 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, Pn(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.

Contents

Rodrigues' formula

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


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

Using the Newton binomial and the equation


\frac{d^n x^m}{dx^n}  = \frac{m!}{(m-n)!} x^{m-n}, \quad\hbox{for}\quad n\le m,

we get the explicit expression


P_n(x) = \frac{1}{2^n \, n!}\sum_{k=\lceil n/2 \rceil}^n (-1)^{n-k} {n \choose k}\frac{(2k)!}{(2k-n)!} x ^{2k-n} .

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


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

Generating function

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


\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


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

where


h\equiv \frac{r}{R}\quad \hbox{and}\quad x \equiv \cos\gamma \equiv \mathbf{r}\cdot\mathbf{R}/(rR).

Obviously the expansion makes sense only if R > r.

Normalization

The polynomials are not normalized to unity, but


\int_{-1}^{1} P_{n}(x) P_{m}(x) dx = \frac{2}{2n+1} \delta_{n m},

where δnm is the Kronecker delta.

Differential equation

The Legendre polynomials are solutions of the Legendre differential equation


\frac{d}{dx} \left[ (1-x^2) \frac{d}{dx} P_n(x) \right] + n(n+1)P_n(x)  
= (1-x^2) \frac{d^2 P_n(x)}{dx^2} - 2 x \frac{dP_n(x)}{dx} + n(n+1) P_n(x) = 0.

This differential equation has another class of solutions: Legendre functions of the second kind Qn(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


\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

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

P_n(-x) = (-1)^n P_n(x). \,

The following condition normalizes the polynomials

P_n(1) = 1. \,

Recurrence Relations

Legendre polynomials satisfy the recurrence relations


(1-x^2) \frac{d}{dx} P_n = (n+1) xP_n - (n+1)  P_{n+1} = n P_{n-1} - nx P_n
 
(n+1) P_{n+1} - (2n+1) x P_n + n P_{n-1} = 0 \,

From these two relations follows easily


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

External link

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

Views
Personal tools