Bessel Functions — Mathematics Knowledge System

Overview

Bessel functions of the first kind Jₙ(x) and second kind Yₙ(x) are the two linearly independent solutions of Bessel's differential equation x²y″ + xy′ + (x²−ν²)y = 0. They arise naturally when solving the Laplace or Helmholtz equation in cylindrical coordinates by separation of variables, and appear throughout physics and engineering: heat flow in a cylinder, vibrations of a circular drum, electromagnetic wave propagation in waveguides, and signal processing. Jₙ is bounded at x=0 for integer n; Yₙ is unbounded at x=0.

Intuition

Bessel functions are to circular/cylindrical geometry what sines and cosines are to rectangular geometry: they are the natural oscillatory basis for radial problems. Just as sin and cos solve y″+y=0 (the harmonic oscillator with constant coefficients), Jₙ and Yₙ solve the radial part of the Helmholtz equation. At large x, Jₙ(x) ≈ √(2/πx) cos(x − nπ/2 − π/4), showing the characteristic decaying oscillation. The zeros of Jₙ serve as the eigenvalues for boundary value problems on disks.

Formal Definition

Definition

Bessel's equation of order ν and its two canonical solutions:

x^2 y'' + x y' + (x^2 - \nu^2)y = 0

Bessel's differential equation

J_\nu(x) = \sum_{m=0}^{\infty} \frac{(-1)^m}{m!\,\Gamma(m+\nu+1)}\left(\frac{x}{2}\right)^{2m+\nu}

Bessel function of the first kind

Y_\nu(x) = \frac{J_\nu(x)\cos(\nu\pi) - J_{-\nu}(x)}{\sin(\nu\pi)}

Bessel function of the second kind (Neumann function)

J_n(x) = \frac{1}{\pi}\int_0^{\pi} \cos(n\theta - x\sin\theta)\,d\theta

Integral representation (integer n)

Worked Examples

For ν=0: J₀(x) = ∑_{m=0}^∞ (−1)^m/(m!)² (x/2)^{2m}.
$J_0(x) = \sum_{m=0}^{\infty} \frac{(-1)^m}{(m!)^2}\left(\frac{x}{2}\right)^{2m}$
First terms: m=0 gives 1; m=1 gives −(x/2)²/1 = −x²/4; m=2 gives (x/2)⁴/4 = x⁴/64.
$J_0(x) = 1 - \frac{x^2}{4} + \frac{x^4}{64} - \frac{x^6}{2304} + \cdots$
Compare with cos(x) = 1 − x²/2 + x⁴/24 − … The differences grow with order; J₀ decays slower for large x.
$J_0(x) \approx \sqrt{\frac{2}{\pi x}}\cos\!\left(x - \frac{\pi}{4}\right) \text{ for large } x$

Answer: J₀(x) = 1 − x²/4 + x⁴/64 − x⁶/2304 + …, a decaying oscillatory function.

Practice Problems

Difficulty 7/10

Verify that J₀ and J₁ satisfy J₀′(x) = −J₁(x).

Difficulty 8/10

What boundary condition do the zeros of J₀ satisfy, and how are they used in solving the heat equation on a disk?

Difficulty 8/10

Determine the general solution to Bessel's equation of order 1 near x=0.

Common Mistakes

Common Mistake

Bessel functions are polynomials.

Bessel functions are infinite series and are transcendental functions, not polynomials (unlike Legendre polynomials which are polynomial).

Common Mistake

J₋ₙ = Jₙ for all ν.

For integer n, J₋ₙ(x) = (−1)ⁿ Jₙ(x); for non-integer ν, J_{−ν} and Jᵥ are linearly independent.

Quiz

Summary

Bessel's equation x²y″+xy′+(x²−ν²)y=0 has solutions Jᵥ (first kind) and Yᵥ (second kind).
Jᵥ(x) = ∑ (−1)^m/[m!Γ(m+ν+1)] (x/2)^{2m+ν} converges for all x.
For large x, Jᵥ(x) ≈ √(2/πx) cos(x − νπ/2 − π/4).
Recurrence: J_{n−1}(x) + J_{n+1}(x) = (2n/x) Jₙ(x).
Zeros of Jₙ are eigenvalues in cylindrical boundary value problems.

References

WebsiteWikipedia — Bessel function