Series Solutions

Learning Objectives

  • Understand how to represent solutions to differential equations as power series.
  • Determine the radius of convergence using the Ratio Test.
  • Classify ordinary and singular points for variable-coefficient DEs.
  • Apply the method of power series around an ordinary point.
  • Apply the Frobenius method around a regular singular point.

When a differential equation cannot be solved by elementary analytical methods (especially when it involves variable coefficients like x2y+y=0x^2 y'' + y = 0), we often assume the solution can be represented as an infinite power series.

Power Series Solution

We assume a solution of the form of an infinite series. Typically, we expand the series around x0=0x_0 = 0 (a Maclaurin series).

Power Series Solution Form

The general assumed form of a power series solution.

y(x)=n=0cn(xx0)ny(x) = \sum_{n=0}^{\infty} c_n (x-x_0)^n

Variables

SymbolDescriptionUnit
cnc_nThe constant coefficients of the series-
x0x_0The point around which the series is expanded-

Maclaurin Series Solution Form

The assumed form of a power series expanded around x0=0x_0 = 0.

y(x)=n=0cnxn=c0+c1x+c2x2+c3x3+y(x) = \sum_{n=0}^{\infty} c_n x^n = c_0 + c_1 x + c_2 x^2 + c_3 x^3 + \dots

Variables

SymbolDescriptionUnit
cnc_nThe constant coefficients of the series-

Convergence and the Ratio Test

A power series is only useful if it converges to a finite value. The interval of xx values for which the series converges is called the interval of convergence, centered at x0x_0, and half its length is the radius of convergence, RR.

Finding the Radius of Convergence

The primary tool to find RR is the Ratio Test. A power series an\sum a_n converges absolutely if the limit of the absolute ratio of successive terms is less than 1.

Applying this to the power series an=cn(xx0)na_n = c_n(x-x_0)^n:

  1. Set up the limit: limncn+1(xx0)n+1cn(xx0)n=limncn+1cnxx0<1\lim_{n \to \infty} \left| \frac{c_{n+1}(x-x_0)^{n+1}}{c_n(x-x_0)^n} \right| = \lim_{n \to \infty} \left| \frac{c_{n+1}}{c_n} \right| |x-x_0| < 1
  2. Solve for xx0|x-x_0|: xx0<limncncn+1=R|x-x_0| < \lim_{n \to \infty} \left| \frac{c_n}{c_{n+1}} \right| = R

If R=0R=0, it converges only at x=x0x=x_0. If R=R=\infty, it converges for all real xx.

Ratio Test for Power Series

Calculates the limit to determine absolute convergence of a power series.

limnan+1an=L<1\lim_{n \to \infty} \left| \frac{a_{n+1}}{a_n} \right| = L < 1

Variables

SymbolDescriptionUnit
ana_nThe n-th term of the series-
LLThe limit value-

Power Series Convergence

Before diving into solving DEs with series, it helps to visualize how a power series converges to a target function as more terms are added within its radius of convergence. This is the foundation of Taylor series approximations.

Interactive Simulation

Use the simulation below to explore how power series solutions converge as terms are added.

Power Series Solution

Solving y=yy' = y with y(0)=1y(0)=1 around the ordinary point x0=0x_0=0.

y(x)n=00xnn!y(x) \approx \sum_{n=0}^{0} \frac{x^n}{n!}

Notice how adding more terms increases the interval where the polynomial approximation closely matches the exact exponential solution.

1
LinearDegree 9
Exact Solution exe^x
Series Approximation
x ∈ [-3, 3]

Ordinary vs. Singular Points

Before solving a second-order linear DE, P(x)y+Q(x)y+R(x)y=0P(x)y'' + Q(x)y' + R(x)y = 0, we must check if the expansion point x0x_0 is ordinary or singular.

Classification of Points

First, rewrite the DE in standard form: y+p(x)y+q(x)y=0y'' + p(x)y' + q(x)y = 0, where p(x)=Q(x)/P(x)p(x) = Q(x)/P(x) and q(x)=R(x)/P(x)q(x) = R(x)/P(x).

  • Ordinary Point: If both p(x)p(x) and q(x)q(x) are analytic (have convergent Taylor series) at x0x_0. Practically, this means P(x0)0P(x_0) \neq 0. Use the standard power series method.
  • Singular Point: If P(x0)=0P(x_0) = 0 (meaning p(x)p(x) or q(x)q(x) is undefined at x0x_0).

    • Regular Singular Point: If the singularities are "mild." Specifically, if (xx0)p(x)(x-x_0)p(x) and (xx0)2q(x)(x-x_0)^2q(x) are both analytic at x0x_0. Use the Frobenius Method.
    • Irregular Singular Point: If the conditions for a regular singular point fail. Solutions behave wildly, and these methods generally do not apply.

Method of Power Series (Around an Ordinary Point)

  1. Assume Solution: Write y=n=0cnxny = \sum_{n=0}^\infty c_n x^n.
  2. Differentiate: Find y=n=1ncnxn1y' = \sum_{n=1}^\infty n c_n x^{n-1} and y=n=2n(n1)cnxn2y'' = \sum_{n=2}^\infty n(n-1) c_n x^{n-2}.
  3. Substitute: Plug y,y,yy, y', y'' into the DE.
  4. Shift Indices: Adjust the summation indices (nkn \to k) so all terms involve the exact same power of xx (e.g., xkx^k) and start at the same lower limit.
  5. Find Recurrence Relation: Combine the sums. By the identity property of series, the coefficient of every power xkx^k must equal zero. This gives an algebraic equation for ck+2c_{k+2} in terms of ckc_k and ck1c_{k-1}, etc.
  6. Solve: Use the recurrence relation to find c2,c3,c4c_2, c_3, c_4\dots in terms of arbitrary constants c0c_0 and c1c_1 (which correspond to initial conditions).

Method of Frobenius (Around a Regular Singular Point)

If x0=0x_0 = 0 is a regular singular point, standard power series fail. Instead, the solution takes the form of a Frobenius series, which adds an arbitrary power rr, where c00c_0 \neq 0 and rr is a constant to be determined (the indicial root).

Frobenius Series Solution

The assumed solution form around a regular singular point.

y(x)=n=0cnxn+r=xrn=0cnxny(x) = \sum_{n=0}^{\infty} c_n x^{n+r} = x^r \sum_{n=0}^{\infty} c_n x^n

Variables

SymbolDescriptionUnit
cnc_nThe constant coefficients of the series-
rrThe indicial root-

Applying the Method of Frobenius

  1. Assume Solution: Write y=n=0cnxn+ry = \sum_{n=0}^{\infty} c_n x^{n+r}.
  2. Find Derivatives: Differentiate twice: y=n=0(n+r)cnxn+r1y' = \sum_{n=0}^{\infty} (n+r) c_n x^{n+r-1} and y=n=0(n+r)(n+r1)cnxn+r2y'' = \sum_{n=0}^{\infty} (n+r)(n+r-1) c_n x^{n+r-2}.
  3. Substitute into DE: Substitute into the differential equation, distribute any xx coefficients, and shift indices to align the powers of xx.
  4. Indicial Equation: Look at the term with the lowest power of xx (usually xr2x^{r-2} or xrx^r) and set its coefficient to zero. Since we assumed c00c_0 \neq 0, this yields a quadratic equation in rr, known as the indicial equation.
  5. Solve for r: Find the roots r1,r2r_1, r_2 (indicial roots). The form of the second linearly independent solution depends heavily on whether r1r2r_1 - r_2 is an integer or zero.
Key Takeaways
  • Power Series Method: Assume y=cnxny = \sum c_n x^n to solve variable-coefficient DEs around an Ordinary Point (P(x0)0P(x_0) \neq 0).
  • Radius of Convergence (RR): Use the Ratio test to determine the interval xx0<R|x - x_0| < R where the infinite series approximation is mathematically valid.
  • Recurrence Relation: The core result of the method. It defines how to calculate higher-order coefficients cnc_n based on previous ones, usually leaving c0c_0 and c1c_1 as arbitrary constants.
  • Frobenius Method: Required when expanding around a Regular Singular Point (P(x0)=0P(x_0) = 0, but singularities are mild). Assumes y=xrcnxny = x^r \sum c_n x^n and requires solving an indicial equation for rr.