4.6.20Ordinary Differential Equations

Legendre's equation and Legendre polynomials (intro)

1,709 words8 min readdifficulty · medium2 backlinks

WHAT is Legendre's equation?

WHY does x=cosθx=\cos\theta appear? In spherical coordinates the angular equation has θ\theta as the variable; substituting x=cosθx=\cos\theta turns it into the form above, with x[1,1]x\in[-1,1]. The endpoints x=±1x=\pm1 are the poles of the sphere — and they are singular points of the ODE.


HOW to derive it (the self-adjoint form first)


HOW to solve: the power-series (Frobenius) method

WHY a series? x=0x=0 is an ordinary point (the coefficient of yy'', which is 1x21-x^2, is nonzero there), so a plain Taylor series y=amxmy=\sum a_m x^m converges and works.

Assume y=m=0amxm.y = \sum_{m=0}^{\infty} a_m x^m.

Then y=mamxm1y' = \sum m a_m x^{m-1} and y=m(m1)amxm2y'' = \sum m(m-1)a_m x^{m-2}. Plug in:

(1x2)m(m1)amxm22xmamxm1+n(n+1)amxm=0.(1-x^2)\sum m(m-1)a_m x^{m-2} - 2x\sum m a_m x^{m-1} + n(n+1)\sum a_m x^m = 0.

Why this step? We re-index each sum to the same power xmx^m:

  • m(m1)amxm2(m+2)(m+1)am+2xm\sum m(m-1)a_m x^{m-2} \to \sum (m+2)(m+1)a_{m+2}x^{m} (shift mm+2m\to m+2).
  • The rest already carry xmx^m.

Collecting the coefficient of xmx^m and setting it to 00: (m+2)(m+1)am+2m(m1)am2mam+n(n+1)am=0.(m+2)(m+1)a_{m+2} - m(m-1)a_m - 2m\,a_m + n(n+1)a_m = 0.

Group the ama_m terms: m(m1)2m+n(n+1)=n(n+1)m(m+1)-m(m-1)-2m+n(n+1) = n(n+1)-m(m+1). So


Building the first few Legendre polynomials

We normalize with the convention Pn(1)=1P_n(1)=1.

Figure — Legendre's equation and Legendre polynomials (intro)

Two short-cut formulas (derived, not dumped)


Quick parity & value facts

  • Pn(x)=(1)nPn(x)P_n(-x)=(-1)^n P_n(x) — even nn even function, odd nn odd.
  • Pn(1)=1P_n(1)=1, Pn(1)=(1)nP_n(-1)=(-1)^n.
  • Degree of PnP_n is exactly nn; leading coefficient (2n)!2n(n!)2\dfrac{(2n)!}{2^n (n!)^2}.



Recall Feynman: explain to a 12-year-old

Imagine you're describing the temperature all over a beach ball using only smooth "shape patterns." There's a flat pattern (same everywhere), a top-vs-bottom pattern, a "ring" pattern, and so on. Legendre polynomials are exactly those building-block patterns. The rule for finding them is one tidy equation, and the cool trick is: only the patterns that don't go crazy at the two poles of the ball are allowed. Add up the right amounts of these patterns and you can describe any smooth temperature on the ball.


Flashcards

Legendre's equation (standard form)
(1x2)y2xy+n(n+1)y=0(1-x^2)y'' - 2xy' + n(n+1)y = 0
Legendre's equation in self-adjoint form
ddx[(1x2)y]+n(n+1)y=0\frac{d}{dx}[(1-x^2)y'] + n(n+1)y = 0
Recurrence relation for series coefficients
am+2=m(m+1)n(n+1)(m+2)(m+1)ama_{m+2} = \dfrac{m(m+1)-n(n+1)}{(m+2)(m+1)}a_m
Why does the series terminate for integer nn?
At m=nm=n the numerator m(m+1)n(n+1)=0m(m+1)-n(n+1)=0, so an+2=0a_{n+2}=0 and all higher same-parity coefficients vanish, leaving a degree-nn polynomial.
P0,P1,P2,P3P_0,P_1,P_2,P_3
1,  x,  12(3x21),  12(5x33x)1,\;x,\;\tfrac12(3x^2-1),\;\tfrac12(5x^3-3x)
Rodrigues' formula
Pn(x)=12nn!dndxn(x21)nP_n(x)=\dfrac{1}{2^n n!}\dfrac{d^n}{dx^n}(x^2-1)^n
Orthogonality relation
11PmPndx=22n+1δmn\int_{-1}^1 P_m P_n\,dx = \dfrac{2}{2n+1}\delta_{mn}
Normalization convention
Pn(1)=1P_n(1)=1 (and Pn(1)=(1)nP_n(-1)=(-1)^n)
Parity of PnP_n
Pn(x)=(1)nPn(x)P_n(-x)=(-1)^n P_n(x)
What is the second solution and why discard it?
Qn(x)Q_n(x), the Legendre function of the second kind; it is unbounded (log singularity) at x=±1x=\pm1, so it's physically rejected.
Why is x=0x=0 an ordinary point but x=±1x=\pm1 singular?
Coefficient (1x2)(1-x^2) of yy'' is nonzero at 00 but vanishes at ±1\pm1 (regular singular points = the poles of the sphere).

Connections

  • Power Series / Frobenius Method — how we built y=amxmy=\sum a_m x^m.
  • Sturm-Liouville Theory — self-adjoint form ⇒ orthogonality and real eigenvalues.
  • Laplace's Equation in Spherical Coordinates — where x=cosθx=\cos\theta and Legendre's equation are born.
  • Fourier Series — same idea, orthogonal basis expansion, polynomials instead of sines.
  • Associated Legendre Functions & Spherical Harmonics — the next step (m0m\ne0 angular modes).
  • Generating Functions112xt+t2=Pn(x)tn\frac{1}{\sqrt{1-2xt+t^2}}=\sum P_n(x)t^n links to the multipole expansion.

Concept Map

leads to

becomes

transforms

rewritten as

Sturm-Liouville gives

justifies

solved by

yields

when n integer

forces

gives

property of

Spherical symmetry physics

Angular equation in theta

Substitute x = cos theta

Legendre equation

Self-adjoint form

Orthogonal solutions

x=0 ordinary point

Power series y = sum a_m x^m

Recurrence relation

Integer n

Series terminates

Legendre polynomial P_n

Hinglish (regional understanding)

Intuition Hinglish mein samjho

Dekho, Legendre's equation tab aati hai jab aap koi physics problem solve karte ho jiska shape sphere (gol ball) jaisa hota hai — jaise planet ki gravity, ya kisi charge ka electric potential. Spherical coordinates use karte ho to angular part hamesha is ek hi equation mein convert ho jaata hai: (1x2)y2xy+n(n+1)y=0(1-x^2)y'' - 2xy' + n(n+1)y = 0, jahan x=cosθx=\cos\theta hota hai. Yaani xx humesha 1-1 se 11 ke beech, aur x=±1x=\pm1 matlab ball ke do poles.

Solve kaise karte hain? Hum y=amxmy=\sum a_m x^m power series maan lete hain aur substitute karke ek recurrence relation nikaalte hain: am+2=m(m+1)n(n+1)(m+2)(m+1)ama_{m+2}=\frac{m(m+1)-n(n+1)}{(m+2)(m+1)}a_m. Sabse beautiful baat — jab nn ek poora integer hota hai, to m=nm=n par numerator zero ho jaata hai, series ruk jaati hai, aur humein ek saaf-suthra polynomial mil jaata hai. Isiliye P0=1P_0=1, P1=xP_1=x, P2=12(3x21)P_2=\frac12(3x^2-1), etc. Jo doosra solution QnQ_n hota hai woh poles par phat jaata hai (unbounded), isliye usko physics mein chhod dete hain.

Kyun important hai? Kyunki yeh polynomials orthogonal hote hain: 11PmPndx=22n+1δmn\int_{-1}^1 P_m P_n\,dx = \frac{2}{2n+1}\delta_{mn}. Iska matlab inko Fourier series ki tarah use karke kisi bhi function ko inke combination mein likh sakte ho. Yeh property self-adjoint (Sturm-Liouville) form se aati hai. So short mein: ek equation yaad karo, terminate hone wali trick samjho, aur poora 3D physics ka angular part aapke haath mein.

Go deeper — visual, from zero

Test yourself — Ordinary Differential Equations

Connections