5.2 — Frobenius Method (Regular Singular Points)

Learning Objectives
  • Classify singular points as regular or irregular
  • Apply the Frobenius method to find series solutions
  • Determine and solve the indicial equation
  • Handle all three cases of indicial root differences

Motivation

The power series method we studied in Section 5.1 works beautifully at ordinary points, but fails when we approach a singular point. What happens when we need a solution valid near a singularity?

The Frobenius method extends the power series approach by allowing an initial fractional power. Instead of:

$$y = \sum_{n=0}^{\infty} a_n x^n$$

we use:

$$y = (x - x_0)^r \sum_{n=0}^{\infty} a_n (x - x_0)^n = \sum_{n=0}^{\infty} a_n (x - x_0)^{n+r}$$

where $r$ is determined by an indicial equation. This method works at regular singular points and will give us access to solutions involving logarithms and all the classical special functions of physics and engineering.

Theory: Regular Singular Points and the Frobenius Method

Definition: Regular vs Irregular Singular Points

For the differential equation:

$$y'' + p(x)y' + q(x)y = 0$$

A point $x = x_0$ is a regular singular point if both of the following limits exist and are finite:

Regular Singular Point
$$(x - x_0)p(x) \quad \text{and} \quad (x - x_0)^2 q(x)$$

are both analytic at $x = x_0$ (or equivalently, can be expressed as convergent power series centered at $x_0$).

If at least one of these conditions fails, $x_0$ is an irregular singular point, where Frobenius typically does not apply.

The Frobenius Theorem

Theorem (Frobenius, 1873)

If $x_0$ is a regular singular point of $y'' + p(x)y' + q(x)y = 0$, then there exists at least one solution of the form:

$$y = (x - x_0)^r \sum_{n=0}^{\infty} a_n (x - x_0)^n$$

valid in some interval $0 < |x - x_0| < R$, where the exponent $r$ is determined by the indicial equation.

The Indicial Equation

Let:

$$p_0 = \lim_{x \to x_0} (x - x_0)p(x) \quad \text{and} \quad q_0 = \lim_{x \to x_0} (x - x_0)^2 q(x)$$

The indicial equation is:

Indicial Equation
$$r(r-1) + p_0 r + q_0 = 0$$

or equivalently: $r^2 + (p_0 - 1)r + q_0 = 0$

This is a quadratic equation in $r$. Its roots $r_1$ and $r_2$ (with $r_1 \geq r_2$) determine which solutions exist and whether logarithmic terms appear.

Three Cases for the Indicial Roots

Case 1: Distinct Roots Not Differing by an Integer

If $r_1 - r_2 \notin \mathbb{Z}$ (the difference is not an integer), then two independent Frobenius series solutions are:

$$y_1 = (x - x_0)^{r_1} \sum_{n=0}^{\infty} a_n^{(1)} (x - x_0)^n$$ $$y_2 = (x - x_0)^{r_2} \sum_{n=0}^{\infty} a_n^{(2)} (x - x_0)^n$$

General solution: $y = c_1 y_1 + c_2 y_2$

Case 2: Repeated Root

If $r_1 = r_2 = r$ (a double root), the second solution contains a logarithmic term:

$$y_1 = (x - x_0)^r \sum_{n=0}^{\infty} a_n (x - x_0)^n$$ $$y_2 = y_1 \ln(x - x_0) + (x - x_0)^r \sum_{n=1}^{\infty} b_n (x - x_0)^n$$

General solution: $y = c_1 y_1 + c_2 y_2$

Case 3: Roots Differing by a Positive Integer

If $r_1 - r_2 = m \in \mathbb{Z}^+$ (a positive integer), a logarithmic term may appear in the second solution:

$$y_1 = (x - x_0)^{r_1} \sum_{n=0}^{\infty} a_n (x - x_0)^n$$ $$y_2 = A y_1 \ln(x - x_0) + (x - x_0)^{r_2} \sum_{n=0}^{\infty} b_n (x - x_0)^n$$

where $A$ may be zero (no logarithm) or nonzero. The form must be determined by substitution.

Method: Step-by-Step Solution Procedure

  1. Verify $x_0$ is a regular singular point
    Write the equation in the form $y'' + p(x)y' + q(x)y = 0$ (divide by the leading coefficient if necessary). Check that $(x - x_0)p(x)$ and $(x - x_0)^2 q(x)$ are both analytic at $x_0$.
  2. Compute the limit values $p_0$ and $q_0$
    $$p_0 = \lim_{x \to x_0} (x - x_0)p(x), \quad q_0 = \lim_{x \to x_0} (x - x_0)^2 q(x)$$
  3. Write and solve the indicial equation
    $$r(r-1) + p_0 r + q_0 = 0$$ Find the roots $r_1 \geq r_2$.
  4. Determine which case applies
    Check if $r_1 - r_2$ is an integer. This determines the form of the second solution.
  5. For the larger root $r_1$: assume the Frobenius series
    $$y = \sum_{n=0}^{\infty} a_n (x - x_0)^{n+r_1}$$ Substitute into the differential equation.
  6. Find the recurrence relation
    Collect coefficients of like powers of $(x - x_0)$ to obtain a recurrence formula relating $a_n$ to previous coefficients.
  7. Solve for coefficients
    Set $a_0 = 1$ (or any nonzero constant) and use the recurrence to find $a_1, a_2, a_3, \ldots$. This yields the first solution $y_1$.
  8. Find the second solution
    Depending on the case:
    • Case 1: Repeat steps 5–7 using $r_2$
    • Case 2: Use the formula with the logarithmic term
    • Case 3: Attempt the logarithmic form or use reduction of order

Worked Examples

Solution

Step 1: Verify $x = 0$ is regular singular

Divide by $2x^2$:

$$y'' + \frac{1}{2x}y' - \frac{x+1}{2x^2}y = 0$$

Here $p(x) = \frac{1}{2x}$ and $q(x) = -\frac{x+1}{2x^2}$.

Check: $(x - 0)p(x) = x \cdot \frac{1}{2x} = \frac{1}{2}$ ✓ and $(x - 0)^2 q(x) = x^2 \cdot \left(-\frac{x+1}{2x^2}\right) = -\frac{x+1}{2}$ ✓

Both analytic at $x = 0$, so $x = 0$ is regular singular.

Step 2: Compute $p_0$ and $q_0$

$$p_0 = \lim_{x \to 0} \frac{1}{2} = \frac{1}{2}$$ $$q_0 = \lim_{x \to 0} \left(-\frac{x+1}{2}\right) = -\frac{1}{2}$$

Step 3: Indicial equation

$$r(r-1) + \frac{1}{2}r - \frac{1}{2} = 0$$ $$r^2 - r + \frac{1}{2}r - \frac{1}{2} = 0$$ $$r^2 - \frac{1}{2}r - \frac{1}{2} = 0$$ $$2r^2 - r - 1 = 0$$ $$(2r + 1)(r - 1) = 0$$

Roots: $r_1 = 1$ and $r_2 = -\frac{1}{2}$, so $r_1 - r_2 = \frac{3}{2} \notin \mathbb{Z}$ (Case 1).

Step 4: First solution with $r_1 = 1$

Assume $y = x \sum_{n=0}^{\infty} a_n x^n = \sum_{n=0}^{\infty} a_n x^{n+1}$.

Then: $y' = \sum_{n=0}^{\infty} (n+1)a_n x^n$ and $y'' = \sum_{n=0}^{\infty} (n+1)n a_n x^{n-1}$.

Substituting into $2x^2 y'' + xy' - (x+1)y = 0$ and collecting coefficients of $x^n$ yields a recurrence relation. Setting $a_0 = 1$ and solving gives the first Frobenius series solution.

Note: The recurrence relation and explicit coefficients follow standard Frobenius methods. For brevity, we focus on the procedure and classification.

Bessel's Equation: $x^2 y'' + xy' + x^2 y = 0$

Step 1: Verify regularity

Divide by $x^2$:

$$y'' + \frac{1}{x}y' + y = 0$$

Here $p(x) = \frac{1}{x}$ and $q(x) = 1$. Check: $xp(x) = 1$ ✓ and $x^2 q(x) = x^2$ ✓. Both analytic at $x = 0$. Regular singular point ✓

Step 2: Limits

$$p_0 = \lim_{x \to 0} 1 = 1, \quad q_0 = \lim_{x \to 0} x^2 = 0$$

Step 3: Indicial equation

$$r(r-1) + r + 0 = 0$$ $$r^2 = 0$$

Double root: $r = 0$ (Case 2: repeated root with logarithmic second solution).

Step 4: First solution $y_1$ with $r = 0$

Assume $y = \sum_{n=0}^{\infty} a_n x^n$ (ordinary power series). Substitution yields the recurrence:

$$a_n = \frac{-a_{n-2}}{n^2} \quad \text{for } n \geq 2$$

With $a_0 = 1$ and $a_1 = 0$, we get:

$$y_1(x) = J_0(x) = \sum_{n=0}^{\infty} \frac{(-1)^n x^{2n}}{2^{2n}(n!)^2} = 1 - \frac{x^2}{4} + \frac{x^4}{64} - \cdots$$

This is the Bessel function of the first kind, order 0.

Step 5: Second solution $y_2$ (Case 2)

Since we have a repeated root, the second solution is:

$$y_2(x) = y_1(x) \ln x + \sum_{n=1}^{\infty} b_n x^n$$

This involves the Bessel function of the second kind, $Y_0(x)$, which is far more complex to compute directly but is the standard companion solution to $J_0$.

Key Insight

Bessel's equation appears everywhere in physics: circular membranes, cylindrical heat flow, vibrations, electromagnetic waves, etc. The Frobenius method gives us direct access to these fundamental special functions.

Solution

Step 1: Verify regularity

Divide by $x^2$:

$$y'' + \frac{3}{x}y' + \frac{1}{x^2}y = 0$$

Here $p(x) = \frac{3}{x}$ and $q(x) = \frac{1}{x^2}$.

Check: $xp(x) = 3$ ✓ and $x^2 q(x) = 1$ ✓. Regular singular point ✓

Step 2: Limits

$$p_0 = 3, \quad q_0 = 1$$

Step 3: Indicial equation

$$r(r-1) + 3r + 1 = 0$$ $$r^2 + 2r + 1 = 0$$ $$(r + 1)^2 = 0$$

Double root: $r = -1$ (Case 2: repeated root with logarithmic second solution).

Step 4: Solution forms

For a Cauchy-Euler equation with repeated root $r$, the solutions are:

$$y_1 = x^{-1}, \quad y_2 = x^{-1} \ln x$$

General solution:

$$y(x) = c_1 x^{-1} + c_2 x^{-1} \ln x = \frac{c_1 + c_2 \ln x}{x}$$
Verification by Substitution

You can verify that both $y_1 = 1/x$ and $y_2 = (\ln x)/x$ satisfy the equation. Try differentiating and substituting!

Solution

Step 1: Verify regularity

Divide by $x^2$:

$$y'' - \frac{1+x}{x}y' + \frac{1}{x^2}y = 0$$

Here $p(x) = -\frac{1+x}{x}$ and $q(x) = \frac{1}{x^2}$.

Check: $xp(x) = -(1+x)$ ✓ and $x^2 q(x) = 1$ ✓. Regular singular point ✓

Step 2: Limits

$$p_0 = \lim_{x \to 0} (-(1+x)) = -1, \quad q_0 = \lim_{x \to 0} 1 = 1$$

Step 3: Indicial equation

$$r(r-1) - r + 1 = 0$$ $$r^2 - 2r + 1 = 0$$ $$(r - 1)^2 = 0$$

Repeated root $r_1 = r_2 = 1$, so this is Case 2. For a repeated root the logarithm is always present: only one exponent is available, so the second solution must be built as $y_2 = y_1\ln x + x\sum_{n\ge 1} b_n x^n$.

Note: Only in Case 3 ($r_1 - r_2$ a positive integer) do you have to test whether the logarithm actually appears — there its coefficient can turn out to be zero. When the roots are equal, as here, the logarithm can never be avoided.

Solution

Step 1: Verify regularity

Divide by $x^2$:

$$y'' + \frac{1}{x}y' + \left(1 - \frac{1}{4x^2}\right)y = 0$$

Here $p(x) = \frac{1}{x}$ and $q(x) = 1 - \frac{1}{4x^2}$.

Check: $xp(x) = 1$ ✓ and $x^2 q(x) = x^2 - \frac{1}{4}$ ✓. Regular singular point ✓

Step 2: Limits

$$p_0 = 1, \quad q_0 = -\frac{1}{4}$$

Step 3: Indicial equation

$$r(r-1) + r - \frac{1}{4} = 0$$ $$r^2 - \frac{1}{4} = 0$$ $$r = \pm \frac{1}{2}$$

Distinct roots: $r_1 = \frac{1}{2}$ and $r_2 = -\frac{1}{2}$. Difference $r_1 - r_2 = 1 \in \mathbb{Z}$ (Case 3: difference is a positive integer).

Step 4: Solutions

For Bessel's equation of order $\nu$, the recurrence relations yield:

$$y_1 = J_{1/2}(x) = \sqrt{\frac{2}{\pi x}}\,\sin x = \sqrt{\frac{2}{\pi}}\,x^{1/2}\left(1 - \frac{x^2}{6} + \frac{x^4}{120} - \cdots\right)$$ $$y_2 = J_{-1/2}(x) = \sqrt{\frac{2}{\pi x}}\,\cos x = \sqrt{\frac{2}{\pi}}\,x^{-1/2}\left(1 - \frac{x^2}{2} + \frac{x^4}{24} - \cdots\right)$$

Remarkable fact: Although roots differ by an integer, the second solution does not contain a logarithmic term! The general solution is:

$$y(x) = c_1 \sqrt{\frac{2}{\pi x}} \sin x + c_2 \sqrt{\frac{2}{\pi x}} \cos x$$
Warning

When $r_1 - r_2 \in \mathbb{Z}$, you cannot conclude whether a logarithm appears without checking the recurrence relation. The logarithmic term appears sometimes (Case 2 and most of Case 3) but not always (as here).

Example 6 (Case 1): $2xy'' + y' + y = 0$ — roots differing by a non-integer

Solution

Step 1: Check that $x = 0$ is a regular singular point

Divide by $2x$:

$$y'' + \frac{1}{2x}y' + \frac{1}{2x}y = 0 \quad\Longrightarrow\quad p(x) = \frac{1}{2x},\qquad q(x) = \frac{1}{2x}$$ $$xp(x) = \frac{1}{2}, \qquad x^2q(x) = \frac{x}{2}$$

Both are analytic at $x = 0$, so $x = 0$ is a regular singular point with $p_0 = \tfrac{1}{2}$ and $q_0 = 0$.

Step 2: Indicial equation

$$r(r-1) + \tfrac{1}{2}r + 0 = 0 \quad\Longrightarrow\quad r\left(r - \tfrac{1}{2}\right) = 0$$ $$r_1 = \frac{1}{2}, \qquad r_2 = 0, \qquad r_1 - r_2 = \frac{1}{2} \notin \mathbb{Z}$$

This is Case 1: two independent Frobenius series, no logarithm.

Step 3: Substitute $y = \sum_{n\ge 0} a_n x^{n+r}$ and build the recurrence

$$2xy'' = \sum_{n\ge 0} 2(n+r)(n+r-1)a_n x^{n+r-1},\qquad y' = \sum_{n\ge 0}(n+r)a_n x^{n+r-1},\qquad y = \sum_{n\ge 0} a_n x^{n+r}$$

Shift the index in the last sum so every term carries $x^{n+r-1}$, then equate coefficients:

$$a_n\Big[2(n+r)(n+r-1) + (n+r)\Big] + a_{n-1} = 0$$ $$\boxed{\;a_n\,(n+r)\big(2n+2r-1\big) = -a_{n-1},\qquad n \ge 1\;}$$

The $n = 0$ term reproduces the indicial equation $a_0\,r(2r-1) = 0$, confirming $r = \tfrac{1}{2},\,0$.

Step 4: First solution, $r_1 = \tfrac{1}{2}$

Here $(n+\tfrac{1}{2})(2n) = n(2n+1)$, so $a_n = -\dfrac{a_{n-1}}{n(2n+1)}$. With $a_0 = 1$:

$$a_1 = -\frac{1}{3},\qquad a_2 = \frac{1}{30},\qquad a_3 = -\frac{1}{630}$$ $$y_1(x) = x^{1/2}\left(1 - \frac{x}{3} + \frac{x^2}{30} - \frac{x^3}{630} + \cdots\right)$$

Step 5: Second solution, $r_2 = 0$

Now the bracket is $n(2n-1)$, so $a_n = -\dfrac{a_{n-1}}{n(2n-1)}$. With $a_0 = 1$:

$$a_1 = -1,\qquad a_2 = \frac{1}{6},\qquad a_3 = -\frac{1}{90}$$ $$y_2(x) = 1 - x + \frac{x^2}{6} - \frac{x^3}{90} + \cdots$$

General solution:

$$\boxed{y(x) = c_1\,x^{1/2}\left(1 - \frac{x}{3} + \frac{x^2}{30} - \cdots\right) + c_2\left(1 - x + \frac{x^2}{6} - \cdots\right)}$$
Why is there no logarithm? The factor $(n+r)(2n+2r-1)$ never vanishes for $n \ge 1$ at either root, so the recurrence never breaks down. That is precisely what $r_1 - r_2 \notin \mathbb{Z}$ guarantees.
Example 7 (Case 2): $xy'' + y' - y = 0$ — repeated root, logarithm required

Solution

Step 1: Indicial equation

Dividing by $x$ gives $p(x) = 1/x$ and $q(x) = -1/x$, so $xp = 1$ and $x^2q = -x$ are analytic: $x = 0$ is regular singular with $p_0 = 1$, $q_0 = 0$.

$$r(r-1) + r = r^2 = 0 \quad\Longrightarrow\quad r_1 = r_2 = 0$$

A repeated root — Case 2 — so the second solution must contain $\ln x$.

Step 2: Recurrence and first solution

With $y = \sum a_n x^{n+r}$ the coefficient of $x^{n+r-1}$ gives

$$a_n(n+r)^2 = a_{n-1} \quad\xrightarrow{\;r\,=\,0\;}\quad a_n = \frac{a_{n-1}}{n^2}$$

Taking $a_0 = 1$ gives $a_n = \dfrac{1}{(n!)^2}$, so

$$y_1(x) = \sum_{n=0}^{\infty}\frac{x^n}{(n!)^2} = 1 + x + \frac{x^2}{4} + \frac{x^3}{36} + \frac{x^4}{576} + \cdots$$

Step 3: Second solution with the logarithm

$$y_2(x) = y_1(x)\ln x + \sum_{n=1}^{\infty} b_n x^n$$

Substituting, the $\ln x$ terms cancel automatically (because $y_1$ solves the equation) and what remains is

$$2y_1'(x) + \sum_{n\ge 1}\Big[n^2 b_n - b_{n-1}\Big]x^{n-1} = 0 \quad\Longrightarrow\quad \boxed{\;n^2 b_n = b_{n-1} - 2n\,a_n\;}$$

with $b_0 = 0$. Running it:

$$b_1 = -2a_1 = -2,\qquad 4b_2 = b_1 - 4a_2 = -3 \Rightarrow b_2 = -\frac{3}{4},$$ $$9b_3 = b_2 - 6a_3 = -\frac{11}{12} \Rightarrow b_3 = -\frac{11}{108},\qquad b_4 = -\frac{25}{3456}$$

In closed form $b_n = -\dfrac{2H_n}{(n!)^2}$, where $H_n = 1 + \tfrac{1}{2} + \cdots + \tfrac{1}{n}$ is the $n$-th harmonic number:

$$\boxed{y_2(x) = y_1(x)\ln x - 2\sum_{n=1}^{\infty}\frac{H_n}{(n!)^2}x^n = y_1(x)\ln x - \left(2x + \frac{3x^2}{4} + \frac{11x^3}{108} + \cdots\right)}$$
Sanity check: substituting the truncated $y_2$ back into $xy'' + y' - y$ leaves only the terms beyond the truncation order — every coefficient up to that order cancels exactly.
Example 8 (Case 3): $xy'' - y = 0$ — roots differing by an integer, logarithm unavoidable

Solution

Step 1: Indicial equation

Here $p(x) = 0$ and $q(x) = -1/x$, so $xp = 0$ and $x^2q = -x$: regular singular with $p_0 = q_0 = 0$.

$$r(r-1) = 0 \quad\Longrightarrow\quad r_1 = 1,\quad r_2 = 0,\qquad r_1 - r_2 = 1 \in \mathbb{Z}^+$$

This is Case 3. Always start from the larger root.

Step 2: Larger root $r_1 = 1$

The coefficient of $x^{n+r-1}$ gives $a_n(n+r)(n+r-1) = a_{n-1}$, so with $r = 1$:

$$a_n = \frac{a_{n-1}}{n(n+1)} \quad\Longrightarrow\quad a_n = \frac{1}{n!\,(n+1)!}$$ $$y_1(x) = x + \frac{x^2}{2} + \frac{x^3}{12} + \frac{x^4}{144} + \cdots = \sum_{n=0}^{\infty}\frac{x^{n+1}}{n!\,(n+1)!}$$

Step 3: Try the smaller root $r_2 = 0$ — and watch it fail

With $r = 0$ the recurrence reads $a_n\,n(n-1) = a_{n-1}$. At the critical index $n = r_1 - r_2 = 1$:

$$a_1\cdot 1\cdot 0 = a_0 \quad\Longrightarrow\quad 0 = a_0,$$

contradicting $a_0 \ne 0$. No pure Frobenius series exists at $r_2$, so the logarithm is forced.

Step 4: Logarithmic second solution

$$y_2(x) = y_1(x)\ln x + \sum_{n=0}^{\infty} b_n x^{n}$$

Substituting and matching the coefficient of $x^k$ gives

$$\boxed{\;k(k+1)\,b_{k+1} = b_k - (2k+1)a_k\;}$$

At $k = 0$ this forces $b_0 = a_0 = 1$; $b_1$ is free (it only adds a multiple of $y_1$), so take $b_1 = 0$. Then

$$2b_2 = b_1 - 3a_1 = -\frac{3}{2} \Rightarrow b_2 = -\frac{3}{4},\qquad 6b_3 = b_2 - 5a_2 = -\frac{7}{6} \Rightarrow b_3 = -\frac{7}{36},\qquad b_4 = -\frac{35}{1728}$$ $$\boxed{y_2(x) = y_1(x)\ln x + 1 - \frac{3}{4}x^2 - \frac{7}{36}x^3 - \frac{35}{1728}x^4 - \cdots}$$
Contrast with Example 5. Bessel's equation of order $\tfrac{1}{2}$ also has $r_1 - r_2 = 1$, yet there the critical equation was $0\cdot a_1 = 0$ — satisfied automatically, leaving $a_1$ free and producing a second Frobenius series with no logarithm. In Case 3 you must always run the recurrence down to $n = r_1 - r_2$ and read off what the critical equation actually says.

📝 Exam-Style Practice Problems

These problems are similar in style and difficulty to past exam questions. Click each problem to reveal the step-by-step solution.

Practice 1: Classify $x = 2$ and $x = -2$ for $(x^2-4)^2y'' + (x-2)y' + y = 0$

Solution

Step 1: Put the equation in standard form

$$y'' + \underbrace{\frac{x-2}{(x-2)^2(x+2)^2}}_{p(x)}y' + \underbrace{\frac{1}{(x-2)^2(x+2)^2}}_{q(x)}y = 0, \qquad p(x) = \frac{1}{(x-2)(x+2)^2}$$

Both $x = 2$ and $x = -2$ are singular points.

Step 2: Test $x = 2$

$$(x-2)p(x) = \frac{1}{(x+2)^2}, \qquad (x-2)^2q(x) = \frac{1}{(x+2)^2}$$

Both are analytic at $x = 2$ (value $\tfrac{1}{16}$), so $x = 2$ is a regular singular point, with $p_0 = q_0 = \tfrac{1}{16}$.

Step 3: Test $x = -2$

$$(x+2)p(x) = \frac{1}{(x-2)(x+2)}$$

This blows up as $x \to -2$, so it is not analytic there.

$$\boxed{x = 2 \text{ is regular singular};\qquad x = -2 \text{ is irregular singular}}$$
Exam tip: a single failing test is enough. You do not need to check $q$ once $(x-x_0)p(x)$ has already failed to be analytic.
Practice 2: Find the indicial roots and recurrence for $3xy'' + y' - y = 0$

Solution

Step 1: Verify the singular point and read off $p_0, q_0$

Dividing by $3x$: $p(x) = \dfrac{1}{3x}$, $q(x) = -\dfrac{1}{3x}$, so $xp = \tfrac{1}{3}$ and $x^2q = -\tfrac{x}{3}$ are analytic: regular singular, $p_0 = \tfrac{1}{3}$, $q_0 = 0$.

Step 2: Indicial equation

$$r(r-1) + \tfrac{1}{3}r = 0 \quad\Longrightarrow\quad r\left(r - \tfrac{2}{3}\right) = 0 \quad\Longrightarrow\quad r_1 = \frac{2}{3},\; r_2 = 0$$

Since $r_1 - r_2 = \tfrac{2}{3} \notin \mathbb{Z}$, this is Case 1: two Frobenius series, no logarithm.

Step 3: Recurrence

Substituting $y = \sum a_n x^{n+r}$ and collecting the coefficient of $x^{n+r-1}$:

$$a_n\Big[3(n+r)(n+r-1) + (n+r)\Big] = a_{n-1} \quad\Longrightarrow\quad a_n(n+r)(3n+3r-2) = a_{n-1}$$

Step 4: Both solutions

For $r_1 = \tfrac{2}{3}$: $a_n = \dfrac{a_{n-1}}{n(3n+2)}$, giving $a_1 = \tfrac{1}{5}$, $a_2 = \tfrac{1}{80}$, $a_3 = \tfrac{1}{2640}$.

For $r_2 = 0$: $a_n = \dfrac{a_{n-1}}{n(3n-2)}$, giving $a_1 = 1$, $a_2 = \tfrac{1}{8}$, $a_3 = \tfrac{1}{168}$.

$$\boxed{y = c_1x^{2/3}\left(1 + \frac{x}{5} + \frac{x^2}{80} + \frac{x^3}{2640} + \cdots\right) + c_2\left(1 + x + \frac{x^2}{8} + \frac{x^3}{168} + \cdots\right)}$$
Practice 3 (Case 2): Solve $xy'' + (1-x)y' - y = 0$ about $x = 0$

Solution

Step 1: Indicial equation

$p(x) = \dfrac{1-x}{x}$ and $q(x) = -\dfrac{1}{x}$, so $xp = 1-x$ and $x^2q = -x$ are analytic: $p_0 = 1$, $q_0 = 0$.

$$r(r-1) + r = r^2 = 0 \quad\Longrightarrow\quad r_1 = r_2 = 0 \quad(\text{Case 2})$$

Step 2: Recurrence for the first solution

Substituting $y = \sum a_n x^{n+r}$, the coefficient of $x^{n+r-1}$ gives

$$(n+r)^2a_n = (n+r-1)a_{n-1} + a_{n-1} = (n+r)a_{n-1} \quad\Longrightarrow\quad a_n = \frac{a_{n-1}}{n+r}$$

At $r = 0$: $a_n = \dfrac{a_{n-1}}{n} = \dfrac{1}{n!}$, so

$$y_1(x) = \sum_{n=0}^{\infty}\frac{x^n}{n!} = e^{x}$$

Check directly: $xe^x + (1-x)e^x - e^x = 0$. ✓

Step 3: Second solution by reduction of order

With $p(x) = \dfrac{1-x}{x}$ we get $\displaystyle\int p\,dx = \ln x - x$, so $e^{-\int p\,dx} = \dfrac{e^{x}}{x}$ and

$$y_2 = y_1\int \frac{e^{-\int p\,dx}}{y_1^2}dx = e^{x}\int \frac{e^{x}}{x\,e^{2x}}dx = e^{x}\int\frac{e^{-x}}{x}\,dx$$

Expanding $\dfrac{e^{-x}}{x} = \dfrac{1}{x} - 1 + \dfrac{x}{2} - \dfrac{x^2}{6} + \cdots$ and integrating term by term:

$$\int\frac{e^{-x}}{x}dx = \ln x + \sum_{n=1}^{\infty}\frac{(-1)^n x^n}{n\,n!} = \ln x - x + \frac{x^2}{4} - \frac{x^3}{18} + \cdots$$ $$\boxed{y_2(x) = e^{x}\ln x - x - \frac{3x^2}{4} - \frac{11x^3}{36} - \frac{25x^4}{288} - \cdots}$$

exactly the Case 2 form $y_2 = y_1\ln x + \sum_{n\ge1} b_n x^n$ with $b_1 = -1$, $b_2 = -\tfrac{3}{4}$, $b_3 = -\tfrac{11}{36}$, $b_4 = -\tfrac{25}{288}$.

Practice 4 (Case 3, no logarithm): Solve $xy'' + 2y' + xy = 0$

Solution

Step 1: Indicial equation

$xp = 2$ and $x^2q = x^2$ are analytic, so $p_0 = 2$, $q_0 = 0$:

$$r(r-1) + 2r = r(r+1) = 0 \quad\Longrightarrow\quad r_1 = 0,\; r_2 = -1,\qquad r_1 - r_2 = 1 \in \mathbb{Z}^+$$

Case 3 — the logarithm may be needed. We must check the critical index.

Step 2: Recurrence

Collecting the coefficient of $x^{n+r-1}$ in $xy'' + 2y' + xy$:

$$a_n(n+r)(n+r+1) + a_{n-2} = 0 \quad\Longrightarrow\quad a_n = \frac{-a_{n-2}}{(n+r)(n+r+1)}$$

Step 3: The smaller root $r_2 = -1$ — the critical test

At $r = -1$ the factor is $(n-1)n$, which vanishes at $n = 1$. The equation at $n = 1$ reads

$$a_1\cdot 0\cdot 1 + a_{-1} = 0 \quad\Longrightarrow\quad 0 = 0,$$

which is automatically satisfied. So $a_1$ is free, no logarithm is needed, and the single Frobenius series at $r_2 = -1$ already carries both solutions.

Step 4: Sum the series

With $a_0 = 1$, $a_1 = 0$: $a_2 = -\tfrac{1}{2}$, $a_4 = \tfrac{1}{24}$, giving

$$y = x^{-1}\left(1 - \frac{x^2}{2} + \frac{x^4}{24} - \cdots\right) = \frac{\cos x}{x}$$

With $a_0 = 0$, $a_1 = 1$ the odd branch gives $x^{-1}\left(x - \tfrac{x^3}{6} + \tfrac{x^5}{120} - \cdots\right) = \dfrac{\sin x}{x}$.

$$\boxed{y(x) = c_1\frac{\cos x}{x} + c_2\frac{\sin x}{x}}$$
Slick check: setting $u = xy$ turns the equation into $u'' + u = 0$, whose solutions are $u = \cos x$ and $u = \sin x$ — exactly the answer above.

Self-Assessment Quiz (scored)

Test your understanding of the Frobenius method with these multiple-choice questions. After answering all questions, your score will be displayed.

0 / 6 correct
Question 1: Which of the following is a regular singular point?
Question 2: For the indicial equation $r^2 - 3r + 2 = 0$, what are the roots?
Question 3: If the indicial roots are $r_1 = 3$ and $r_2 = 1/2$, which case applies?
Question 4: For Bessel's equation of order 0, $x^2 y'' + xy' + x^2 y = 0$, the indicial equation is:
Question 5: For the Cauchy-Euler equation $x^2 y'' + 3xy' + y = 0$ with repeated root $r = -1$, the general solution is:
Question 6: If $r_1 - r_2 = 0$ (repeated root), which solution form is guaranteed?

Final Score: 0 / 6

Review any incorrect answers above and refer to the Theory and Examples sections for clarification.

Practice Problems

Test your understanding with these multiple-choice questions. Click a choice to see if you're correct.

Q1: Classify the point $x = 0$ for $x^3y'' + x^2y' + y = 0$.

Q2: What is the indicial equation of $2x^2y'' + 3xy' - (1+x)y = 0$ at $x = 0$?

Q3: The indicial roots of an equation are $r_1 = 4$ and $r_2 = 2$. What can you conclude?

Q4: For $2xy'' + y' + y = 0$ the recurrence is $a_n(n+r)(2n+2r-1) = -a_{n-1}$. At the root $r = \tfrac{1}{2}$, what is $a_1/a_0$?

Q5: Which situation guarantees that $\ln(x - x_0)$ appears in the second solution?

Q6: Bessel's equation of order $0$ gives the recurrence $a_n = -a_{n-2}/n^2$ with $a_0 = 1$, $a_1 = 0$. What is $a_4$?

Quick Reference Card

Regular Singular Point Test

For $y'' + p(x)y' + q(x)y = 0$ at $x = x_0$:

(x-x₀)p(x) and (x-x₀)²q(x) must be analytic at $x_0$.

Indicial Equation Form

$$r(r-1) + p_0 r + q_0 = 0$$

where $p_0 = \lim_{x \to x_0}(x-x_0)p(x)$ and $q_0 = \lim_{x \to x_0}(x-x_0)^2 q(x)$

Case 1: Distinct Roots

If $r_1 - r_2 \notin \mathbb{Z}$:

Two Frobenius series, both analytic, no logarithms.

Case 2: Repeated Root

If $r_1 = r_2 = r$:

$y_2 = y_1\ln(x-x_0) + \ldots$ (logarithm always present)

Case 3: Integer Difference

If $r_1 - r_2 = m \in \mathbb{Z}^+$:

Logarithm may or may not appear; use recurrence to check.

Frobenius Series Form

$$y = (x-x_0)^r \sum_{n=0}^{\infty} a_n(x-x_0)^n$$

Exponent $r$ determined by indicial equation roots.