1 Definition and basic properties
Legendre polynomials are a sequence of polynomials usually denoted by \(P_0(x), P_1(x), P_2(x), \dots\), defined on the interval \([-1,1]\). They form a classical orthogonal family and appear naturally in problems where rotational symmetry is present. Because they are orthogonal with respect to the constant weight function, they provide an especially convenient basis for representing functions on a finite interval.
A standard characterization is that each \(P_n\) is a polynomial of degree \(n\) satisfying Legendre’s differential equation. This differential equation, together with a normalization such as \(P_n(1)=1\), determines the family uniquely. The sequence begins with simple low-degree polynomials and grows in a highly regular pattern that is captured by recurrence relations and generating formulas.
1.1 Legendre’s differential equation
Legendre polynomials satisfy \[ (1-x^2)y''-2xy'+n(n+1)y=0, \] where \(n\) is a nonnegative integer. For each such \(n\), there is a polynomial solution of degree \(n\). The equation has regular singular points at \(x=\pm 1\), and the polynomial solutions are the ones that remain finite on the interval \([-1,1]\).
This differential equation is central to the theory because it identifies Legendre polynomials as eigenfunctions of a Sturm–Liouville problem. That viewpoint explains their orthogonality and many of their analytic properties.
1.2 Rodrigues’ formula
A compact formula for Legendre polynomials is Rodrigues’ formula: \[ P_n(x)=\frac{1}{2^n n!}\frac{d^n}{dx^n}(x^2-1)^n. \] It gives an explicit way to construct each polynomial by differentiating a simple algebraic expression. The formula is often used to prove orthogonality, derive recurrence relations, and establish other identities.
Because \((x^2-1)^n\) vanishes to high order at \(x=\pm 1\), repeated differentiation produces a polynomial with well-controlled endpoint behavior. This makes Rodrigues’ formula one of the most useful starting points for computations and proofs.
1.3 Recurrence relations
Legendre polynomials satisfy several recurrence relations that allow higher-degree polynomials to be built from lower-degree ones. These relations are efficient computationally and reveal the internal structure of the family.
1.3.1 Three-term recurrence
The most important recurrence is \[ (n+1)P_{n+1}(x)=(2n+1)xP_n(x)-nP_{n-1}(x), \] valid for \(n\ge 1\). This three-term relation is typical of orthogonal polynomial systems. It permits stable iterative generation of the sequence and is widely used in numerical work.
1.3.2 Derivative relations
Differentiation produces further identities connecting neighboring orders. One useful form is \[ (1-x^2)P_n'(x)=n\bigl(P_{n-1}(x)-xP_n(x)\bigr). \] Such formulas link the derivative of a Legendre polynomial to lower-degree terms and help in studying endpoint behavior, differential equations, and related functions.
1.4 Initial polynomials
The first few Legendre polynomials are \[ P_0(x)=1,\quad P_1(x)=x,\quad P_2(x)=\tfrac12(3x^2-1), \] \[ P_3(x)=\tfrac12(5x^3-3x),\quad P_4(x)=\tfrac18(35x^4-30x^2+3). \] These examples illustrate the alternating structure and the normalization \(P_n(1)=1\). They also show that \(P_n\) contains only powers of the same parity as \(n\).
2 Orthogonality and normalization
One of the defining features of Legendre polynomials is their orthogonality on \([-1,1]\) with respect to the constant weight. This property makes them especially effective for expansions and approximation.
2.1 Orthogonality on [-1,1]
For \(m\neq n\), \[ \int_{-1}^1 P_m(x)P_n(x)\,dx=0. \] When \(m=n\), the integral is positive and gives the square norm of the polynomial. Orthogonality follows from the differential equation and the self-adjoint form of the associated operator.
This relation means that different Legendre polynomials behave like perpendicular vectors in a function space. As a result, coefficients in Legendre expansions can be extracted cleanly by projection.
2.2 Norm formulas
The \(L^2\)-norm of \(P_n\) on \([-1,1]\) is \[ \int_{-1}^1 P_n(x)^2\,dx=\frac{2}{2n+1}. \] This simple closed form is one reason Legendre polynomials are so convenient in analysis. It also leads directly to normalized variants and Fourier-like coefficient formulas.
2.3 Normalized Legendre polynomials
A normalized version is often defined by \[ \widehat{P}_n(x)=\sqrt{\frac{2n+1}{2}}\,P_n(x), \] so that \[ \int_{-1}^1 \widehat{P}_m(x)\widehat{P}_n(x)\,dx=\delta_{mn}. \] Normalized polynomials are useful in Hilbert-space settings, where orthonormal bases simplify expansions and error estimates.
3 Generating functions and explicit formulas
Legendre polynomials admit compact generating and explicit representations. These formulas connect the family to combinatorics, special functions, and coefficient extraction methods.
3.1 Generating function
The generating function is \[ \frac{1}{\sqrt{1-2xt+t^2}}=\sum_{n=0}^\infty P_n(x)t^n, \]
| valid for \( | t | <1\). Expanding the left-hand side in powers of \(t\) produces the entire sequence of polynomials. |
|---|
This representation is extremely useful because it encodes all Legendre polynomials in a single expression. It also leads to addition formulas, explicit coefficients, and many identities by formal manipulation.
3.2 Explicit polynomial expressions
Each \(P_n(x)\) can be written as a finite sum: \[ P_n(x)=\sum_{k=0}^{\lfloor n/2\rfloor} (-1)^k \frac{(2n-2k)!}{2^n k!(n-k)!(n-2k)!}x^{n-2k}. \] This formula makes the polynomial structure transparent. The powers of \(x\) decrease by steps of two, reflecting the parity of the index.
3.3 Hypergeometric representation
Legendre polynomials can also be expressed using the Gauss hypergeometric function: \[ P_n(x)={}_2F_1\!\left(-n,n+1;1;\frac{1-x}{2}\right). \] This form places them within the broader theory of special functions. It is especially valuable for analytic continuation, transformation formulas, and comparisons with related orthogonal polynomial families.
4 Structural properties
Legendre polynomials possess a rich internal structure that is visible in symmetry, zeros, and coefficient patterns. These properties are important both theoretically and numerically.
4.1 Parity
Legendre polynomials satisfy \[ P_n(-x)=(-1)^n P_n(x). \] Thus even-indexed polynomials are even functions and odd-indexed polynomials are odd functions. This parity property explains why only matching powers appear in each explicit polynomial.
4.2 Zeros and interlacing
For \(n\ge 1\), \(P_n\) has exactly \(n\) simple real zeros, all lying in the open interval \((-1,1)\). Moreover, the zeros of \(P_n\) and \(P_{n+1}\) interlace. This means that between any two consecutive zeros of one polynomial lies a zero of the next.
Such behavior is characteristic of orthogonal polynomials and is closely tied to approximation theory and quadrature. The simplicity and location of the zeros are especially important for numerical methods.
4.3 Leading coefficients
The leading coefficient of \(P_n\) is \[ \frac{(2n)!}{2^n(n!)^2}. \] This coefficient grows rapidly with \(n\) and determines the polynomial’s highest-order behavior. It is often used in normalization and asymptotic estimates.
4.4 Addition and symmetry formulas
Legendre polynomials satisfy addition-type identities that relate values at composite arguments and angles. In angular form, such formulas often appear through the cosine of the angle between vectors. They reflect the rotational symmetry underlying many applications.
Symmetry relations also link \(P_n\) to expressions involving rotations and coordinate changes. These identities are especially useful in spherical problems and in the theory of harmonic expansions.
5 Legendre series
Functions on \([-1,1]\) can be expanded in a series of Legendre polynomials, analogous to Fourier series. Such expansions are fundamental in approximation and in solving differential equations.
5.1 Expansion in Legendre polynomials
A suitable function \(f\) may be written as \[ f(x)\sim \sum_{n=0}^\infty a_n P_n(x). \] The sequence \(\{P_n\}\) forms a basis for many function spaces, so the series can represent broad classes of functions. The approximation improves as more terms are included.
5.2 Coefficient formulas
When the expansion is valid in an \(L^2\) sense, the coefficients are \[ a_n=\frac{2n+1}{2}\int_{-1}^1 f(x)P_n(x)\,dx. \] This formula follows directly from orthogonality. It allows each coefficient to be computed independently, which is a major advantage in analysis and numerical practice.
5.3 Convergence properties
The convergence of a Legendre series depends on the smoothness of the underlying function. Smooth functions often admit rapid convergence, while functions with limited regularity converge more slowly. Under suitable assumptions, the series converges in mean square and may converge pointwise away from singularities or discontinuities.
5.4 Parseval-type identities
For square-integrable functions, Legendre coefficients satisfy an energy identity analogous to Parseval’s theorem: \[
| \int_{-1}^1 | f(x) | ^2\,dx |
|---|---|---|
| =\sum_{n=0}^\infty \frac{2}{2n+1} | a_n | ^2, |
\] when the expansion is taken in the standard Legendre basis. This relation quantifies how the function’s total size is distributed among its polynomial modes.
6 Connections with other special functions
Legendre polynomials sit at the center of a larger network of special functions. Many related families can be obtained by generalization or specialization.
6.1 Associated Legendre functions
Associated Legendre functions arise by differentiating Legendre polynomials and extending the resulting expressions to higher orders. They are commonly denoted \(P_n^m(x)\) and appear when separating variables in spherical coordinates.
These functions preserve much of the structure of the original family while introducing an additional index. They are essential in harmonic analysis on the sphere and in many physical applications.
6.2 Spherical harmonics
Spherical harmonics are built from associated Legendre functions and trigonometric factors. In problems with spherical symmetry, Legendre polynomials appear in the axisymmetric case, where the dependence reduces to a single angular variable.
The appearance of Legendre polynomials in spherical harmonics explains their central role in angular decompositions, multipole expansions, and solutions of Laplace-type equations in spherical coordinates.
6.3 Gegenbauer and Jacobi polynomials
Legendre polynomials are a special case of Jacobi polynomials: \[ P_n(x)=P_n^{(0,0)}(x), \] and also of Gegenbauer polynomials after a change of parameters. This places them within the broad Askey scheme of orthogonal polynomials.
These relationships allow many identities for Legendre polynomials to be derived from more general formulas. They also clarify how the Legendre family fits into a hierarchy of classical orthogonal systems.
6.4 Bessel function relations
In certain asymptotic and integral contexts, Legendre polynomials are related to Bessel functions. Such connections often arise when analyzing large-degree behavior or when approximating angular integrals.
Although Legendre and Bessel functions belong to different classes of special functions, they frequently meet in harmonic analysis and in the study of wave phenomena.
7 Applications
Legendre polynomials are widely used in mathematics, physics, and computational methods. Their orthogonality and recurrence properties make them practical tools for both symbolic and numerical work.
7.1 Solving boundary value problems
When a differential equation is posed on an interval or in spherical geometry, Legendre polynomials often arise as separated solutions. They are especially useful in boundary value problems with symmetry around an axis.
Because the polynomials satisfy a Sturm–Liouville equation, they fit naturally into eigenfunction expansions. This makes them effective for representing solutions and enforcing boundary conditions.
7.2 Potential theory
In potential theory, Legendre polynomials appear in expansions of harmonic functions. They are commonly used to express potentials generated by axisymmetric sources.
Their role here is closely tied to the geometry of the sphere. By decomposing a potential into Legendre modes, one can isolate contributions of different angular complexity.
7.3 Numerical quadrature
Legendre polynomials are a foundation of high-accuracy quadrature rules on \([-1,1]\). Their zeros and orthogonality are especially valuable for constructing efficient integration formulas.
7.3.1 Gauss-Legendre quadrature
Gauss-Legendre quadrature uses the zeros of \(P_n\) as nodes. With suitable weights, it integrates polynomials of degree up to \(2n-1\) exactly. This is one of the most effective classical methods for numerical integration on a finite interval.
7.3.2 Spectral methods
In spectral methods, Legendre polynomials serve as basis functions for approximating solutions of differential equations. Their global support and orthogonality make them well suited for smooth problems requiring high precision.
These methods often combine Legendre expansions with matrix formulations of differential operators. The result is a powerful framework for solving ordinary and partial differential equations.
8 Advanced topics
More advanced studies of Legendre polynomials focus on asymptotic formulas, integral representations, operator-theoretic interpretations, and identities that connect neighboring degrees.
8.1 Asymptotic behavior
For large \(n\), Legendre polynomials exhibit oscillatory behavior in the interior of \((-1,1)\) and more delicate scaling near the endpoints. Asymptotic formulas describe their size and phase in these regimes.
Such results are important in approximation theory and in numerical analysis, where understanding the behavior of high-degree polynomials is essential for stability and error estimates.
8.2 Integral representations
Legendre polynomials admit several integral representations, some of which are useful for deriving estimates or proving identities. These formulas often involve trigonometric or contour-integral expressions.
Integral forms can highlight symmetry and can be adapted to analytic continuation. They provide another perspective on the same family of functions.
8.3 Christoffel-Darboux formula
The Christoffel-Darboux formula gives a compact expression for sums of products of orthogonal polynomials. For Legendre polynomials, it provides an efficient way to evaluate partial sums and kernel functions.
This identity is important in approximation theory, quadrature, and the analysis of projection operators. It also clarifies the reproducing behavior of Legendre expansions.
8.4 Differential operator viewpoint
Legendre polynomials may be viewed as eigenfunctions of a second-order differential operator of Sturm–Liouville type. In that setting, the operator is self-adjoint on \([-1,1]\) with an appropriate domain and boundary behavior.
This viewpoint unifies many properties of the polynomials, including orthogonality, completeness, and the structure of the spectrum. It also connects Legendre polynomials to broader themes in spectral theory and functional analysis.