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:

(xx0)P(x)and(xx0)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 (xx0)P(x)=pn(xx0)n(x - x_0)P(x) = \sum p_n(x-x_0)^n and (xx0)2Q(x)=qn(xx0)n(x - x_0)^2 Q(x) = \sum q_n(x-x_0)^n. The indicial equation is:

r(r1)+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)=(xx0)r1n=0an(xx0)n,a00,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<xx0<R0 < |x - x_0| < R. The nature of the second solution depends on the difference r1r2r_1 - r_2:

  1. r1r2Zr_1 - r_2 \notin \mathbb{Z}: Second solution y2=(xx0)r2bn(xx0)ny_2 = (x-x_0)^{r_2}\sum b_n(x-x_0)^n.
  2. r1r2Z>0r_1 - r_2 \in \mathbb{Z}_{>0}: Second solution may involve a logarithmic term: y2=Cy1lnxx0+(xx0)r2bn(xx0)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=y1lnxx0+(xx0)r1bn(xx0)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.