Legendre and Bessel Equations: Special Functions Intro

Legendre's equation, (1−x2)y″ − 2xy′ + n(n+1)y = 0, and Bessel's equation, x2y″ + xy′ + (x2 − ν2)y = 0, are the two most famous differential equations whose solutions are special functions: the Legendre polynomials Pn(x) and the Bessel functions Jν(x). These functions are not built from sines or exponentials; they are defined by the equations themselves, exactly as sin x could be defined by y″ + y = 0.

They matter because separation of variables in spherical coordinates produces Legendre's equation (angular part) and in cylindrical coordinates produces Bessel's equation (radial part).

Key Concepts

• Legendre polynomials: when n is a non-negative integer, one solution is a polynomial of degree n: P0 = 1, P1 = x, P2 = (3x2−1)/2, P3 = (5x3−3x)/2.
• Bessel functions: Jν(x) oscillates like a decaying cosine, with amplitude falling off as 1/√x; ν is the order.
• Orthogonality: ∫−11 PnPm dx = 0 for n ≠ m, and Bessel functions are orthogonal with weight x — so they serve as bases for expansions, like Fourier series.
• Rodrigues' formula: Pn(x) = (1/(2nn!))·dn/dxn[(x2−1)n] generates every Legendre polynomial by differentiation.

Formulas

• Legendre: (1−x2)y″ − 2xy′ + n(n+1)y = 0. Bessel: x2y″ + xy′ + (x2−ν2)y = 0.
• Rodrigues: Pn(x) = (1/(2nn!)) dn/dxn (x2−1)n.
• Bessel series: Jν(x) = Σm=0∞ (−1)m/(m!·Γ(m+ν+1))·(x/2)2m+ν.
• Orthogonality: ∫−11 Pn(x)Pm(x)dx = 0 (n ≠ m); = 2/(2n+1) (n = m).

Worked Example

Build P2(x) from Rodrigues' formula and verify it satisfies Legendre's equation with n = 2. P2 = (1/(22·2!))·d2/dx2[(x2−1)2] = (1/8)·d2/dx2(x4 − 2x2 + 1). First derivative: 4x3 − 4x. Second: 12x2 − 4. Dividing by 8: P2(x) = (3x2−1)/2. Check: y = (3x2−1)/2, y′ = 3x, y″ = 3. Then (1−x2)(3) − 2x(3x) + 2·3·(3x2−1)/2 = 3 − 3x2 − 6x2 + 9x2 − 3 = 0. It satisfies the equation exactly, confirming the polynomial solution for n = 2.

Common Misconceptions

• Bessel functions are just fancy sines. They oscillate, but with decaying amplitude (~1/√x) and unevenly spaced zeros — genuinely different functions with their own tables.
• Legendre gives polynomials for every n. Only for non-negative integers n; other values give infinite series that blow up at x = ±1 and are usually rejected in physics.
• Special functions are artificial inventions. They are as concrete as sin x — sin x is simply the special function of y″ + y = 0 that students meet first.

Real-World Applications

• Multipole expansions of gravitational and electrostatic potentials use Legendre polynomials — Earth's slightly non-spherical gravity field is a Pn series.
• Quantum angular momentum: spherical harmonics, built on associated Legendre functions, give atomic orbital shapes.
• Bessel functions describe vibrating circular membranes (drum modes), microwave waveguides, and diffraction through circular apertures.

Legendre's and Bessel's equations are the canonical singular ODEs of spherical and cylindrical symmetry, solved by Frobenius series into the Legendre polynomials and Bessel functions. orthogonality lets them expand arbitrary functions like a Fourier series. Together they are the working vocabulary of potentials, vibrations, and quantum angular momentum.