ConceptComplete

Regular Singular Points and the Frobenius Method

When a singular point is "regular," the Frobenius method provides a systematic extension of power series techniques by allowing fractional and logarithmic terms.


Classification of Singular Points

Definition9.2Regular and irregular singular points

A singular point x0x_0 of y′′+P(x)y′+Q(x)y=0y'' + P(x)y' + Q(x)y = 0 is regular if:

(x−x0)P(x)and(x−x0)2Q(x)(x - x_0)P(x) \quad \text{and} \quad (x - x_0)^2 Q(x)

are both analytic at x0x_0. Equivalently, P(x)P(x) has at most a simple pole and Q(x)Q(x) has at most a double pole at x0x_0.

A singular point that is not regular is called irregular.

ExampleClassifying singular points
  1. Bessel equation x2y′′+xy′+(x2−ν2)y=0x^2 y'' + xy' + (x^2 - \nu^2)y = 0: P=1/xP = 1/x (simple pole), Q=(x2−ν2)/x2Q = (x^2 - \nu^2)/x^2 (double pole). So x=0x = 0 is a regular singular point.

  2. Euler equation x2y′′+αxy′+βy=0x^2 y'' + \alpha xy' + \beta y = 0: P=α/xP = \alpha/x, Q=β/x2Q = \beta/x^2. Regular singular point at x=0x = 0.

  3. The equation x3y′′+y=0x^3 y'' + y = 0: P=0P = 0, Q=1/x3Q = 1/x^3 (triple pole). So x=0x = 0 is an irregular singular point. The Frobenius method does not apply here.


The Indicial Equation

Definition9.3Indicial equation

Near a regular singular point x0x_0, write (x−x0)P(x)=∑pn(x−x0)n(x - x_0)P(x) = \sum p_n(x-x_0)^n and (x−x0)2Q(x)=∑qn(x−x0)n(x - x_0)^2 Q(x) = \sum q_n(x-x_0)^n. The indicial equation is:

r(r−1)+p0r+q0=0,r(r-1) + p_0 r + q_0 = 0,

whose roots r1,r2r_1, r_2 (with Re(r1)≥Re(r2)\text{Re}(r_1) \geq \text{Re}(r_2)) are called the exponents (or indices) at x0x_0.

Theorem9.2Frobenius theorem

Near a regular singular point x0x_0, there exists at least one solution of the form:

y1(x)=(x−x0)r1∑n=0∞an(x−x0)n,a0≠0,y_1(x) = (x - x_0)^{r_1} \sum_{n=0}^{\infty} a_n (x - x_0)^n, \quad a_0 \neq 0,

converging for 0<∣x−x0∣<R0 < |x - x_0| < R. The nature of the second solution depends on the difference r1−r2r_1 - r_2:

  1. r1−r2∉Zr_1 - r_2 \notin \mathbb{Z}: Second solution y2=(x−x0)r2∑bn(x−x0)ny_2 = (x-x_0)^{r_2}\sum b_n(x-x_0)^n.
  2. r1−r2∈Z>0r_1 - r_2 \in \mathbb{Z}_{>0}: Second solution may involve a logarithmic term: y2=Cy1ln⁡∣x−x0∣+(x−x0)r2∑bn(x−x0)ny_2 = Cy_1 \ln|x-x_0| + (x-x_0)^{r_2}\sum b_n(x-x_0)^n.
  3. r1=r2r_1 = r_2: Second solution necessarily has a logarithm: y2=y1ln⁡∣x−x0∣+(x−x0)r1∑bn(x−x0)ny_2 = y_1 \ln|x-x_0| + (x-x_0)^{r_1}\sum b_n(x-x_0)^n.

Applications

ExampleBessel equation via Frobenius

For x2y′′+xy′+(x2−ν2)y=0x^2y'' + xy' + (x^2 - \nu^2)y = 0 with ν≥0\nu \geq 0: the indicial equation at x=0x = 0 is r2−ν2=0r^2 - \nu^2 = 0, giving r1=νr_1 = \nu, r2=−νr_2 = -\nu.

Substituting y=xν∑anxny = x^\nu \sum a_n x^n yields the recurrence an+2=−an/((n+2)(n+2+2ν))a_{n+2} = -a_n / ((n+2)(n+2+2\nu)), giving the Bessel function:

Jν(x)=(x2)ν∑k=0∞(−1)kk! Γ(ν+k+1)(x2)2k.J_\nu(x) = \left(\frac{x}{2}\right)^\nu \sum_{k=0}^{\infty} \frac{(-1)^k}{k!\,\Gamma(\nu+k+1)} \left(\frac{x}{2}\right)^{2k}.

When 2ν∈Z2\nu \in \mathbb{Z}, the second solution involves logarithmic terms (Bessel functions of the second kind YνY_\nu).

RemarkThe point at infinity

The point at infinity can be studied via the substitution t=1/xt = 1/x. The equation transforms, and x=∞x = \infty corresponds to t=0t = 0. Classifying t=0t = 0 for the transformed equation determines whether x=∞x = \infty is ordinary, regular singular, or irregular singular. For example, the Bessel equation has an irregular singular point at ∞\infty, which accounts for the oscillatory behavior of Bessel functions for large xx.