4.6.19Ordinary Differential Equations

Bessel's equation and Bessel functions (intro, physical relevance)

1,867 words8 min readdifficulty · medium1 backlinks

What is Bessel's equation?

WHAT the pieces mean:

  • x2yx^2 y'' and xyx y' together come from the Laplacian in polar/cylindrical coordinates.
  • The ν2-\nu^2 term comes from the angular separation constant (how many waves fit around the circle).
  • The +x2+x^2 term is what makes it oscillate — without it you'd get a power-law (Euler) equation.

WHERE it comes from (derive it, don't memorize it)

Take the 2-D Helmholtz / vibrating membrane equation 2u+k2u=0\nabla^2 u + k^2 u = 0 in polar coordinates: urr+1rur+1r2uθθ+k2u=0.u_{rr} + \frac{1}{r}u_r + \frac{1}{r^2}u_{\theta\theta} + k^2 u = 0.

Step 1 — Separate variables. Try u(r,θ)=R(r)Θ(θ)u(r,\theta) = R(r)\,\Theta(\theta). Why this step? Circular symmetry suggests radial and angular behaviour are independent.

Substitute and divide by RΘ/r2R\Theta/r^2: r2R+rRR+r2k2=ΘΘ.\frac{r^2 R'' + r R'}{R} + r^2 k^2 = -\frac{\Theta''}{\Theta}.

Step 2 — Set each side to a constant ν2\nu^2. Why? Left depends only on rr, right only on θ\theta; equal for all r,θr,\theta ⇒ both are constant.

Angular: Θ+ν2Θ=0Θ=cosνθ, sinνθ\Theta'' + \nu^2\Theta = 0 \Rightarrow \Theta = \cos\nu\theta,\ \sin\nu\theta. Single-valuedness (Θ(θ+2π)=Θ(θ)\Theta(\theta+2\pi)=\Theta(\theta)) forces ν=n\nu = n integer.

Radial: r2R+rR+(r2k2ν2)R=0.r^2 R'' + r R' + (r^2 k^2 - \nu^2)R = 0.

Step 3 — Rescale x=krx = kr. Then ddr=kddx\dfrac{d}{dr} = k\dfrac{d}{dx}, and the kk's cancel cleanly: x2R+xR+(x2ν2)R=0.x^2 R'' + x R' + (x^2 - \nu^2)R = 0. That is Bessel's equation. The 1rur\tfrac{1}{r}u_r term is the source of the xRxR' term — geometry, not magic.


Solving it: the Frobenius series (first principles)

Because x=0x=0 is a regular singular point, we use a Frobenius series y=m0amxm+sy = \sum_{m\ge0} a_m x^{m+s}.

Step 1 — Indicial equation. Plug the lowest power in. Substituting gives, for the lowest term, [s(s1)+sν2]a0=(s2ν2)a0=0s=±ν.[s(s-1)+s-\nu^2]\,a_0 = (s^2-\nu^2)a_0 = 0 \Rightarrow s = \pm\nu. Why this step? The smallest power must vanish on its own; this fixes the leading exponent.

Step 2 — Recurrence. Collecting xm+sx^{m+s}: [(m+s)2ν2]am+am2=0    am=am2(m+s)2ν2.[(m+s)^2-\nu^2]\,a_m + a_{m-2} = 0 \;\Rightarrow\; a_m = -\frac{a_{m-2}}{(m+s)^2-\nu^2}. Why? The +x2y+x^2 y shifts index by 2, linking ama_m to am2a_{m-2}. Odd terms vanish (a1=0a_1=0).

Step 3 — Take s=νs=\nu, normalize. Choosing a0=12νΓ(ν+1)a_0 = \dfrac{1}{2^\nu \Gamma(\nu+1)} gives the standard Bessel function of the first kind:

The second solution (from s=νs=-\nu, or a log term when ν\nu is an integer) is the Bessel function of the second kind Yν(x)Y_\nu(x), which blows up like lnx\ln x or xνx^{-\nu} at x=0x=0.


Key behaviours (what the graphs tell you)


Worked examples


Common mistakes


Recall Feynman: explain to a 12-year-old

Imagine ripples when you drop a stone in a round pond. On a guitar string the wiggles are even and tidy — that's sine. But on a round pond the ripple has to spread out over a bigger and bigger ring as it travels away from the center, so it gets weaker the further out it goes. Bessel functions are the special "ripple shapes" that describe exactly this: they wiggle up and down like a wave, but each wiggle is a bit smaller than the last because the circle is getting bigger. Drums, pipes, and lenses all use these round-shaped waves.


Active-recall flashcards

Bessel's equation of order ν\nu
x2y+xy+(x2ν2)y=0x^2y'' + xy' + (x^2-\nu^2)y=0
Where does the xyxy' term physically come from?
The 1rur\frac1r u_r in the polar Laplacian (area grows like rr).
Indicial roots of Bessel's equation
s=±νs=\pm\nu (from s2ν2=0s^2-\nu^2=0).
Recurrence relation for coefficients
am=am2/[(m+ν)2ν2]a_m = -\,a_{m-2}/[(m+\nu)^2-\nu^2]; odd terms vanish.
Series for Jν(x)J_\nu(x)
m0(1)mm!Γ(m+ν+1)(x/2)2m+ν\sum_{m\ge0}\frac{(-1)^m}{m!\,\Gamma(m+\nu+1)}(x/2)^{2m+\nu}
Why drop YνY_\nu on a solid disc?
YνY_\nu diverges at x=0x=0; physical displacement must stay finite.
Value J0(0)J_0(0) and Jn(0)J_n(0) for n1n\ge1
J0(0)=1J_0(0)=1, Jn(0)=0J_n(0)=0.
Large-xx behaviour of JνJ_\nu
2/(πx)cos(xνπ/2π/4)\sqrt{2/(\pi x)}\cos(x-\nu\pi/2-\pi/4) — decaying oscillation.
What sets a drum's allowed frequencies?
Zeros αν,k\alpha_{\nu,k} of JνJ_\nu, via Jν(ka)=0J_\nu(ka)=0.
J1/2(x)J_{1/2}(x) in elementary form
2/(πx)sinx\sqrt{2/(\pi x)}\,\sin x.
Why are drum overtones inharmonic?
Ratios of J0J_0 zeros (2.405, 5.520, …) are not integers.

Connections

  • Frobenius method — the engine that produced the JνJ_\nu series.
  • Regular singular points — why x=0x=0 needs Frobenius, not Taylor.
  • Separation of variables — how Bessel's eqn appears from PDEs.
  • Laplacian in polar and cylindrical coordinates — origin of the 1/r1/r.
  • Gamma function — generalizes m!m! for non-integer orders.
  • Sturm-Liouville theory — orthogonality of JνJ_\nu used in Fourier-Bessel series.
  • Legendre's equation — the spherical cousin of Bessel.

Concept Map

leads to

1/r term from area growth

angular part

single-valued forces

radial part rescale x = kr

sets nu^2 term

regular singular point at 0

lowest power

+x^2 shifts index

take s = nu, normalize

role

Circular cylindrical symmetry

Helmholtz equation in polar coords

Separate variables R times Theta

Theta'' + nu^2 Theta = 0

Integer order n

Bessel's equation

Frobenius series

Indicial equation s = plus minus nu

Recurrence links a_m to a_m-2

Bessel function J_nu

Sines and cosines for round things

Hinglish (regional understanding)

Intuition Hinglish mein samjho

Dekho, Bessel's equation tab aati hai jab problem mein round/cylindrical symmetry ho — jaise dhol ki membrane, round pipe mein heat, ya optical fibre mein light. Jab hum PDE ko separate karte hain polar coordinates mein, toh radial part clean sine-cosine nahi rehta, kyunki polar Laplacian mein ek 1rur\frac1r u_r term aata hai. Yahi 1/r1/r term equation ko x2y+xy+(x2ν2)y=0x^2y'' + xy' + (x^2-\nu^2)y=0 bana deta hai. Simple words mein: Bessel functions "round cheezon ke liye sine-cosine" hain.

Solve kaise karte? Kyunki x=0x=0 ek regular singular point hai, hum Frobenius series lagate hain. Indicial equation se s=±νs=\pm\nu milta hai, aur recurrence relation se coefficients. Isse banta hai Jν(x)J_\nu(x) (first kind) — jo center pe finite rehta hai — aur Yν(x)Y_\nu(x) (second kind) — jo x=0x=0 pe infinity tak ud jaata hai. Isliye yaad rakho: "J stays, Y flies" — solid disc problem mein YνY_\nu hamesha hata do, kyunki center pe displacement finite hona chahiye.

Physical relevance kya hai? JνJ_\nu ek damped cosine jaisa dikhta hai — oscillate karta hai par amplitude 1/x1/\sqrt{x} ke rate se ghatti hai (energy badi circle pe spread ho rahi hai). Drum ki allowed frequencies JνJ_\nu ke zeros se decide hoti hain. Aur mast baat: ye zeros (2.405, 5.520, ...) integer ratio mein nahi hote, isiliye dhol ka sound guitar string jaisa "harmonic" nahi lagta. Exam mein bas equation pehchaano, ν\nu nikaalo, aur boundary conditions (finiteness + edge =0=0) lagao — kaam ho gaya.

Go deeper — visual, from zero

Test yourself — Ordinary Differential Equations

Connections