Legendre Polynomials and Bessel Functions Easy Explanation

Legendre polynomials, written P(n)(x), are a family of polynomials that solve Legendre's differential equation, which appears whenever a physics problem has spherical symmetry. Bessel functions, written J(n)(x), solve Bessel's equation, which appears whenever a problem has cylindrical or circular symmetry.

These functions matter because nature is full of spheres and cylinders: the electric potential of a charged sphere, the gravitational field of a planet, the vibration of a drumhead, and the radiation pattern of an antenna are all described by expansions in Legendre polynomials or Bessel functions.

Key Concepts

• Legendre's equation: (1 minus x squared) y double-prime minus 2x y prime + n(n+1) y = 0, with integer n giving the polynomial P(n)(x).
• First Legendre polynomials: P0(x) = 1, P1(x) = x, P2(x) = (3x squared minus 1)/2, P3(x) = (5x cubed minus 3x)/2.
• Bessel's equation: x squared y double-prime + x y prime + (x squared minus n squared) y = 0, with solutions J(n)(x).
• Orthogonality: the integral from minus 1 to 1 of P(m)(x) P(n)(x) dx equals zero when m and n differ, which lets any function on that interval be expanded in Legendre polynomials.

Formulas

• Rodrigues formula: P(n)(x) = 1 / (2 to the power n times n factorial) times the n-th derivative of (x squared minus 1) to the power n.
• Recurrence relation: (n+1) P(n+1)(x) = (2n+1) x P(n)(x) minus n P(n minus 1)(x).
• Orthogonality: integral from minus 1 to 1 of P(m) P(n) dx = 2 / (2n+1) when m = n, and 0 otherwise.
• Expansion coefficients: a(n) = (2n+1)/2 times the integral from minus 1 to 1 of f(x) P(n)(x) dx.

Worked Example

Expand f(x) = x squared on the interval from minus 1 to 1 as a sum of Legendre polynomials.

1. Write x squared = a0 P0 + a1 P1 + a2 P2, since x squared is even and degree 2, so higher and odd terms vanish.
2. Use the coefficient formula with P0 = 1: a0 = (1/2) times the integral from minus 1 to 1 of x squared dx = (1/2) times (2/3) = 1/3.
3. For a2, use P2 = (3x squared minus 1)/2: a2 = (5/2) times the integral of x squared times (3x squared minus 1)/2 dx from minus 1 to 1.
4. The integral equals (1/2) times [ (6/5) minus (2/3) ] = (1/2) times (8/15) = 4/15, and (5/2) times (4/15) = 2/3.
5. Check: (1/3)P0 + (2/3)P2 = 1/3 + x squared minus 1/3 = x squared. Exact.

Common Misconceptions

• "Special functions are just fancy names for approximations." Legendre polynomials are exact polynomials, and the Bessel series converges for every finite x; they are exact solutions, not approximations.
• "Any normalization works." The convention P(n)(1) = 1 is universal in physics; other normalizations silently break recurrences and tables.
• "J(n) is the only Bessel solution." Bessel's equation is second order, so it also has the Neumann functions Y(n), which blow up at x = 0 and are needed whenever the origin is excluded, like waves outside a cylinder.

Real-World Applications

• Multipole expansions in electrostatics and gravity write potentials as sums of P(n)(cos theta), separating monopole, dipole, quadrupole parts.
• Drum and loudspeaker diaphragm vibrations use J0 and J1 zeros to find the allowed frequencies.

Legendre polynomials and Bessel functions are the natural language of spherical and cylindrical problems; their recurrences, zeros, and orthogonality turn hard geometries into routine algebra.