**First‐order equations**. The validity of term‐by‐term differentiation of a power series within its interval of convergence implies that first‐order differential equations may be solved by assuming a solution of the form

substituting this into the equation, and then determining the coefficients *c _{n} *.

**Example 1**: Find a power series solution of the form

for the differential equation

Substituting

into the differential equation yields

Now, write out the first few terms of each series,

and combine like terms:

Since the pattern is clear, this last equation may be written as

*In order for this equation to hold true for all x, every coefficient on the left‐hand side must be zero*. This means *c* _{1} = 0, and for all *n* ≥ 2,

This last equation defines the **recurrence relation** that holds for the coefficients of the power series solution:

Since there is no constraint on *c* _{0}, *c* _{0} is an arbitrary constant, and it is already known that *c* _{1} = 0. The recurrence relation above says *c* _{2} = ½ *c* _{0} and *c* _{3} = ⅓ *c* _{1}, which equals 0 (because *c* _{1} does). In fact, it is easy to see that every coefficient *c _{n} *with

*n*odd will be zero. As for

*c*

_{4}, the recurrence relation says

and so on. Since all *c _{n} *with

*n*odd equal 0, the desire power series solution is therefore

Note that the general solution contains one parameter ( *c* _{0}), as expected for a first‐order differential equation. This power series is unusual in that it is possible to express it in terms of an elementary function. Observe:

It is easy to check that *y* = *c* _{0} *e* ^{x2 }/ ^{2} is indeed the solution of the given differential equation, *y*′ = *xy*. Remember: Most power series cannot be expressed in terms of familiar, elementary functions, so the final answer would be left in the form of a power series.

**Example 2**: Find a power series expansion for the solution of the IVP

Substituting

into the differential equation yields

or, collecting all the terms on one side,

Writing out the first few terms of the series yields

or, upon combining like terms,

Now that the pattern is clear, this last equation can be written

*In order for this equation to hold true for all x, every coefficient on the left‐hand side must be zero*. This means

The last equation defines the recurrence relation that determines the coefficients of the power series solution:

The first equation in (*) says *c* _{1} = *c* _{0}, and the second equation says *c* _{2} = ½(1 + *c* _{1}) = ½(1 + *c* _{0}). Next, the recurrence relation says

and so on. Collecting all these results, the desired power series solution is therefore

Now, the initial condition is applied to evaluate the parameter *c* _{0}:

Therefore, the power series expansion for the solution of the given IVP is

If desired, it is possible to express this in terms of elementary functions. Since

equation (**) may be written

which does indeed satisfy the given IVP, as you can readily verify.

**Second‐order equations**. The process of finding power series solutions of homogeneous second‐order linear differential equations is more subtle than for first‐order equations. Any homogeneous second‐order linear differential equation may be written in the form

If both coefficient functions *p* and *q* are analytic at *x* _{0}, then *x* _{0} is called an **ordinary point** of the differential equation. On the other hand, if even one of these functions fails to be analytic at *x* _{0}, then *x* _{0} is called a **singular point**. Since the method for finding a solution that is a power series in *x* _{0} is considerably more complicated if *x* _{0} is a singular point, attention here will be restricted to power series solutions at ordinary points.

**Example 3**: Find a power series solution in *x* for the IVP

Substituting

into the differential equation yields

The solution may now proceed as in the examples above, writing out the first few terms of the series, collecting like terms, and then determining the constraints on the coefficients from the emerging pattern. Here's another method.

The first step is to re‐index the series so that each one involves *x ^{n} *. In the present case, only the first series must be subjected to this procedure. Replacing

*n*by

*n*+ 2 in this series yields

Therefore, equation (*) becomes

The next step is to rewrite the left‐hand side in terms of a *single* summation. The index *n* ranges from 0 to ∞ in the first and third series, but only from 1 to ∞ in the second. Since the common range of all the series is therefore 1 to ∞, the single summation which will help replace the left‐hand side will range from 1 to ∞. Consequently, it is necessary to first write (**) as

and then combine the series into a single summation:

*In order for this equation to hold true for all x, every coefficient on the left‐hand side must be zero*. This means 2 *c* _{2} + *c* _{0} = 0, and for *n* ≥ 1, the following recurrence relation holds:

Since there is no restriction on *c* _{0} or *c* _{1}, these will be arbitrary, and the equation 2 *c* _{2} + *c* _{0} = 0 implies *c* _{2} = −½ *c* _{0}. For the coefficients from *c* _{3} on, the recurrence relation is needed:

The pattern here isn't too difficult to discern: *c _{n} *= 0 for all odd

*n*≥ 3, and for all even

*n*≥ 4,

This recurrence relation can be restated as follows: for all *n* ≥ 2,

The desired power series solution is therefore

As expected for a second‐order differential equation, the general solution contains two parameters ( *c* _{0} and *c* _{1}), which will be determined by the initial conditions. Since *y*(0) = 2, it is clear that *c* _{0} = 2, and then, since *y*′(0) = 3, the value of *c* _{1} must be 3. The solution of the given IVP is therefore

**Example 4**: Find a power series solution in *x* for the differential equation

Substituting

into the given equation yields

or

Now, all series but the first must be re‐indexed so that each involves *x ^{n} *:

Therefore, equation (*) becomes

The next step is to rewrite the left‐hand side in terms of a *single* summation. The index *n* ranges from 0 to ∞ in the second and third series, but only from 2 to ∞ in the first and fourth. Since the common range of all the series is therefore 2 to ∞, the single summation which will help replace the left‐hand side will range from 2 to ∞. It is therefore necessary to first write (**) as

and then combine the series into a single summation:

Again, in order for this equation to hold true for all *x*, every coefficient on the left‐hand side must be zero. This means *c* _{1} + 2 *c* _{2} = 0, 2 *c* _{2} + 6 *c* _{3} = 0, and for *n* ≥ 2, the following recurrence relation holds:

Since there is no restriction on *c* _{0} or *c* _{1}, these will be arbitrary; the equation *c* _{1} + 2 *c* _{2} = 0 implies *c* _{2} = −½ *c* _{1}, and the equation 2 *c* _{2} + 6 *c* _{3} = 0 implies *c* _{3} = −⅓ *c* _{2} = −⅓(‐½ *c* _{1}) = ⅙ *c* _{1}. For the coefficients from *c* _{4} on, the recurrence relation is needed:

The desired power series solution is therefore

Determining a specific pattern to these coefficients would be a tedious exercise (note how complicated the recurrence relation is), so the final answer is simply left in this form.