Historical Context & Motivation
Throughout the history of mathematics, scientists and engineers have encountered differential equations that simply could not be solved using standard algebraic techniques. A differential equation is an equation that involves an unknown function and its derivatives — for example, an equation that says something about how fast a quantity changes. Methods like separation of variables or integrating factors work beautifully for certain simple types, but many important equations in physics and engineering resist these approaches entirely. The question arose: how can we find solutions when our usual toolkit breaks down?
The answer came from the world of power series — infinite sums of terms of the form a₀ + a₁x + a₂x² + a₃x³ + ⋯, similar to polynomials but with infinitely many terms. You've seen this idea in calculus when studying Taylor series. Mathematicians realized that even when a differential equation has no neat closed-form solution like y = e^x or y = sin(x), you can often express the solution as an infinite series. This breakthrough opened up entire families of equations to analysis.
The central question that this lesson addresses is deceptively simple: given a differential equation that we cannot solve with familiar techniques, how do we assume a power series form for the solution and systematically determine its coefficients? And crucially, how do we know when this approach is guaranteed to work? The concept of an ordinary point provides the answer.
Core Principles & Definitions
Before we can use series to solve differential equations, we need to understand a few key ideas. These definitions tell us when the method works and how to set it up properly. Think of them as a checklist you run through before starting the actual computation. Each term is defined below, and they build on one another, so it helps to read them in order.
Power Series
Analytic Function
Ordinary Point
Singular Point
Radius of Convergence
Recurrence Relation
Visualizing the Series Solution Process
The diagram below shows the overall workflow of finding a series solution near an ordinary point. Understanding this flowchart will guide you through every problem of this type. Each stage builds on the previous one, so the order matters.
Notice how the entire process is systematic. You don't need to guess or be creative — each step follows logically from the previous one. The hardest part is usually Step 3, where you align powers of (x − x₀) using index shifts (re-labeling the counting variable in a sum so that both sums contain the same power of x). We'll walk through that in detail in the worked example later.
Mathematical Framework
Let's lay out the key equations you'll use. We begin with a second-order linear ODE — an equation involving an unknown function y, its first derivative y′ (rate of change), and its second derivative y″ (rate of change of the rate of change) — and convert it into a form where the series approach can be applied. We introduce each definition carefully, since precision here is exactly what guarantees the method works.
a₂(x)y″ + a₁(x)y′ + a₀(x)y = 0, divide every term by a₂(x) first.Identifying Ordinary vs. Singular Points
Before you invest time in a series solution, you need to determine whether your chosen center x₀ is an ordinary point. This classification is the most important preliminary step. The key practical rule: write the equation in standard form (y″ + P(x)y′ + Q(x)y = 0), then check whether P(x) and Q(x) are defined and well-behaved (analytic) at x₀. For most equations you'll encounter, this simply means checking that no denominator equals zero at x₀. The diagram below compares several common differential equations and identifies the nature of key points.
The key practical takeaway from this classification is the minimum radius of convergence. When you find a series solution at an ordinary point x₀, that series is guaranteed to converge in an interval whose radius extends at least to the nearest singular point. For Equation 3 above, centering the series at x₀ = 0 gives a minimum radius of convergence of 1, since the nearest singular points are at x = −1 and x = 1 (both at distance 1 from 0).
| Equation (Standard Form) | Singular Points | x₀ = 0 Status | Min. Radius of Convergence at x₀ = 0 |
|---|---|---|---|
| y″ + xy′ + y = 0 | None | Ordinary | ∞ |
| y″ + (1/x)y′ + y = 0 | x = 0 | Singular | N/A (use Frobenius) |
| (1−x²)y″ − 2xy′ + 6y = 0 | x = ±1 | Ordinary | 1 |
| x²y″ + xy′ + (x²−1)y = 0 | x = 0 | Singular | N/A (use Frobenius) |
Worked Example: Solving y″ + y = 0 via Power Series
Let's apply the series method to an equation whose answer we already know: y″ + y = 0. We know the solutions are sin(x) and cos(x), but let's see how the series approach recovers them. This builds confidence in the method before we tackle unfamiliar equations. The notation Σ (n=0 to ∞) means we sum over all whole numbers n starting from 0; think of it as adding up infinitely many terms according to a pattern.
Strengths, Limitations & Comparisons
The power series method is a powerful tool, but like any method, it has situations where it shines and situations where it struggles. Understanding these helps you choose the right approach for each problem.
| Strengths | Limitations |
|---|---|
| Works when standard methods (undetermined coefficients, variation of parameters) fail | Only guaranteed at ordinary points; singular points require the Frobenius method |
| Produces a general solution with two arbitrary constants | Result is an infinite series, not always a recognizable closed-form function |
| Systematic and algorithmic — no guesswork needed | Algebra can become tedious for equations with non-constant P(x), Q(x) |
| Guaranteed convergence within a known radius | Series may converge slowly, requiring many terms for numerical accuracy |
| Can handle variable-coefficient equations that are common in physics | Does not directly apply to nonlinear differential equations |
Connection to Advanced Theory
The series solution at ordinary points is just the beginning of a larger family of techniques studied in Differential Equations. Understanding where this introductory method fits helps you see the bigger picture and prepares you for more advanced topics you will encounter later in the course and in upper-level mathematics and physics.
| Feature | Series at Ordinary Points (This Lesson) | Frobenius Method (Next Topic) |
|---|---|---|
| Point type | Ordinary point only | Regular singular point |
| Assumed form | y = Σ aₙ(x − x₀)ⁿ | y = xʳ Σ aₙxⁿ (extra xʳ factor, r found from an auxiliary equation) |
| Extra equation needed | None — just a recurrence relation | Indicial equation determines the exponent r before the recurrence relation can be found |
| Number of series solutions | Always two independent solutions | Sometimes only one series solution; a second may involve logarithmic terms |
| Famous examples | Airy equation, Hermite equation, Legendre equation | Bessel equation, hypergeometric equation |
Many of the most important functions in science — Bessel functions (used to model wave patterns in cylinders and heat flow), Legendre polynomials (used in gravitational field calculations), and Hermite polynomials (used in quantum mechanics to describe electron behavior) — were originally discovered through series solution techniques. These functions don't have simple formulas like sin or eˣ, but their series representations make them fully computable and well understood. Mastering the ordinary-point method is the essential first step toward working with these powerful tools in physics and engineering.
Practice Problems
Lesson Summary
This lesson introduces series solutions of differential equations, a technique for solving second-order linear ODEs whose variable coefficients rule out standard methods. When a second-order linear ODE y″ + P(x)y′ + Q(x)y = 0 has variable coefficients that prevent standard solution methods, a power series solution y = Σ aₙ(x − x₀)ⁿ can be assumed. The method works reliably at an ordinary point — a value x₀ where both P(x) and Q(x) are analytic (well-behaved, with no division-by-zero issues). To classify a point, first write the equation in standard form (leading coefficient of y″ equal to 1), then check whether P(x₀) and Q(x₀) are defined and analytic.
The solution procedure is systematic: assume a series, compute its derivatives, substitute into the ODE, align powers of x using index shifts (re-labeling the summation variable so both sums contain the same power of x), and set each coefficient of xⁿ to zero to obtain a recurrence relation. The free constants a₀ and a₁ generate two independent solutions. The series converges at least up to the nearest singular point, giving a guaranteed minimum radius of convergence. This technique is the gateway to more advanced methods — including the Frobenius method — and to important special functions such as Legendre polynomials, Bessel functions, and Hermite polynomials that appear throughout physics and engineering.