Power Series Solutions
Power Series Solutions
1
Chapter 5
Note: These lecture notes aim to present a clear and crisp presentation of some topics in Ordinary
Differential Equations. Comments/suggestions are welcome on the e-mail: [email protected] to
Dr. Suresh Kumar.
Similarly, by Cauchy’s root test the power series (5.1) converges if lim |an |1/n |x − x0 | < 1,
n→∞
that is, if |x − x0 | < R where R = lim |an |−1/n .
n→∞
∞
X
Ex. xn (R = 1. So the power series converges for −1 < x < 1.
n=0
∞
X xn
Ex. (R = ∞. So the power series converges for all x.)
n=0
n!
2
∞
X
Ex. n!xn (R = 0. So the power series converges only for x = 0.)
n=0
Now suppose that the power series (5.1) converges to f (x) for |x − x0 | < R, that is,
∞
X
f (x) = an (x − x0 )n = a0 + a1 (x − x0 ) + a2 (x − x0 )2 + a3 (x − x0 )3 + ........ (5.2)
n=0
Then it can be proved that f (x) possesses derivatives of all orders in |x − x0 | < R. Also, the series
can be differentiated termwise in the sense that
∞
X
f ′ (x) = nan (x − x0 )n−1 = a1 + 2a2 (x − x0 ) + 3a3 (x − x0 )2 + ........,
n=1
∞
X
′′
f (x) = n(n − 1)an (x − x0 )n−2 = 2a2 + 3.2a3 (x − x0 ) + ........,
n=2
and so on, and each of the resulting series converges for |x − x0 | < R. The successive differentiated
series suggest that an = f (n) (0)/n!. Also, the power series (5.2) can be integrated termwise provided
the limits of integration lie inside the interval of convergence.
X∞
If we have another power series bn (x − x0 )n converging to g(x) for |x − x0 | < R, that is,
n=0
∞
X
g(x) = bn (x − x0 )n = b0 + b1 (x − x0 ) + b2 (x − x0 )2 + b3 (x − x0 )3 + ........, (5.3)
n=0
then (5.2) and (5.3) can be added or subtracted termwise, that is,
∞
X
f (x) ± g(x) = (an ± bn )(x − x0 )n = (a0 ± b0 ) + (a1 ± b1 )(x − x0 ) + (a2 ± b2 )(x − x0 )2 + ........
n=0
3
that converges to f (x). On the other hand, if a convergent power series is given, then it is not
always possible to find/recognize its sum function. In fact, very few power series have sums that
are elementary functions.
∞
X f (n) (x0 )
If the power series (x − x0 )n converges to f (x) for all values of x in some neigh-
n=0
n!
bourhood of x0 (open interval containing x0 ), then f (x) is said to be analytic at x0 and the power
series is called Taylor series of f (x) at x0 . Notice that f (x) is analytic at each point in the interval
∞
X f (n) (x0 )
of convergence (x0 − R, x0 + R) of the power series (x − x0 )n .
n=0
n!
4
5.2 Power series solution
The exact methods that we have learned in Chapter 2 and Chapter 3 are applicable to only se-
lected class of DE. There are DE such as the Bessel DE, which can not be solved by exact methods.
Solutions of such DE can be found in the form of power series. We start with a simple example.
∞
X ∞
X
n−1
nan x − an xn = 0, (5.6)
n=1 n=0
which must be an identity in x since (5.4) is, by assumption, a solution of the given DE. So
coefficients of all powers of x must be zero. In particular, equating to 0 the coefficient of xn−1 , the
lowest degree term in x, we obtain
1
nan − an−1 = 0 or an = an−1 .
n
Substituting n = 1, 2, 3, ...., we get
a1 = a0 ,
1 1
a2 = a1 = a0 ,
2 2!
1 1
a3 = a2 = a0 ,
3 3!
and so on. Plugging the values of a1 , a2 , ..... into (5.4), we get
1 1
y = a0 + a0 x + a0 x2 + a0 x3 + ........,
2! 3!
x2 x3
= a0 1 + x + + + ..............
2! 3!
2 3
Let us examine the validity of this solution. We know that the power series 1 + x + x2! + x3! +
.............. converges for all x. It implies that the term by term differentiation carried out in (5.5)
is valid for all x. Similarly, the difference
of the two series (5.4)and (5.5) considered in (5.6) is
2 3
valid for all x. It follows that y = a0 1 + x + x2! + x3! + .............. is a valid solution of the given
x2 x3
DE for all x. Also, we know that ex = 1 + x + 2!
+ 3!
+ ............... So y = a0 ex is general solution
of the DE y ′ − y = 0, as expected.
5
5.2.1 Ordinary and regular singular points
Consider a second order homogeneous LDE
If the functions p(x) and q(x) are analytic at x = x0 , then x0 is called an ordinary point of the DE
(5.7).
If p(x) and/or q(x) fail to be analytic at x0 but (x − x0 )p(x) and (x − x0 )2 q(x) are analytic at
x0 , then x0 is called a regular singular point.
If x0 is neither an ordinary point nor a regular singular point, then it is called an irregular
singular point.
Caution: In the corresponding video lecture, limit approach is used to explain the nature of
points, which is not given in the text book. So DO NOT follow the limit approach in solving the
problems. To show analyticity of a function f (x) at a point x = x0 , express it in the form of a
power series in x − x0 :
f (x) = a0 + a1 (x − x0 ) + a2 (x − x0 )2 + ..........
Here p(x) = 1/x and q(x) = 2/x2 . These are not defined at x = 0, and hence not analytic at
x = 0. So x = 0 is not an ordinary point of the given DE.
Next,
So xp(x) and x2 q(x) both are analytic at x = 0. Hence x = 0 is a regular singular point of the
given DE.
Here p(x) = sin x/x2 , q(x) = 1/x2 . These are clearly not analytic at x = 0. But
x2 q(x) = 1
This shows that xp(x) and x2 q(x) are analytic at x = 0. Hence x = 0 is a regular singular point.
6
Ex. Test the nature of x = 0 for the DE x3 y ′′ + x4 y ′ + y = 0.
Here q(x) = 1/x3 . It is not defined at x = 0, and hence not analytic at x = 0. So x = 0 is not
an ordinary point of the given DE.
Next, x2 q(x) = 1/x is also not defined at x = 0, and hence not analytic at x = 0. So x = 0 is
not a regular singular point of the given DE. It is an irregular singular point.
The above theorem asserts that there exists a unique power series solution of the form
∞
X
y(x) = an (x − x0 )n = a0 + a1 (x − x0 ) + a2 (x − x0 )2 + a3 (x − x0 )3 + ........,
n=0
about the ordinary point x0 satisfying the initial conditions y(x0 ) = a0 and y ′ (x0 ) = a1 . The con-
stants a2 , a3 and so on are determined in terms of a0 or a1 as illustrated in the following examples.
7
1 1
a3 = a1 = a1 ,
3.2 3!
1 1
a4 = a2 = a0 ,
4.3 4!
1 1
a5 = a4 = a1 ,
5.4 5!
and so on. Plugging the values of a3 , a4 , a5 ..... into (5.8), we get
1 1 1 1
y = a0 + a1 x + a0 x2 + a1 x3 + a0 x4 + a1 x5 + ........,
2! 3! 4! 5!
1 2 1 4 1 3 1 5
= a0 1 + x + x + .............. + a1 x + x + x + .............. ,
2! 4! 3! 5!
the required power series solution of the given DE. We know that (ex + e−x )/2 = 1 + 2!1 x2 + 4!1 x4 +
.............. and (ex − e−x )/2 = x + 3!1 x3 + 5!1 x5 + ................ So the power series solution becomes
y = c1 ex + c2 e−x , where c1 = (a0 + a1 )/2 and c1 = (a0 − a1 )/2, which is the same solution of
y ′′ − y = 0 as we obtain by exact method.
Substituting the power series solution (5.9) into the given DE, we get
∞
X ∞
X ∞
X
2 n−2 n−1
(1 + x ) an n(n − 1)x +x an nx − an xn = 0.
n=2 n=1 n=0
∞
X ∞
X
=⇒ an n(n − 1)xn−2 + an [n(n − 1) + n − 1]xn = 0.
n=2 n=0
∞
X ∞
X
n−2
=⇒ an n(n − 1)x + an (n − 1)(n + 1)xn = 0.
n=2 n=0
8
a5 = 0,
3 3
a6 = − a4 = a0 ,
6 6.4.2
and so on.
Plugging the values of a2 , a3 , a4 , a5 , a6 and so on into (5.9), we get
1 1 3
y = a0 + a1 x + a0 x2 + 0.x3 − a0 x4 + 0.x5 + a0 x6 + ........,
2 4.2 6.4.2
1 2 1 4 3 6
= a0 1 + x − x + x .............. + a1 x,
2 4.2 6.4.2
the required power series solution of the given differential equation.
Theorem: If x0 is a regular singular point of a DE y ′′ +p(x)y ′ +q(x)y = 0, then there exists at least
X∞
one power series solution of the form y = an (x − x0 )n+r (a0 ̸= 0), where r is some root of the
n=0
quadratic equation r2 + (p0 − 1)r + q0 = 0, known as indicial equation, with p0 = lim (x − x0 )p(x)
x→x0
and q0 = lim (x − x0 )2 q(x).
x→x0
Remark:P∞ The above n+r theorem by Frobenius guarantees at least one power series solution of the
form n=0 an (x − x0 ) (a0 ̸= 0) of the DE y ′′ + p(x)y ′ + q(x)y = 0, which we call Frobenius
solution. If the roots of the indicial equation do not differ by an integer, we get two LI Frobenious
solutions. In case, there exists only one Frobenious solution, it corresponds to larger root of the
indicial equation. The other LI solution depends on the nature of roots of the indicial equation as
illustrated in the following examples.
a0 (r − 1)(2r + 1) = 0 or (r − 1)(2r + 1) = 0.
9
Therefore, roots of the indicial equation are r = 1, −1/2, which do not differ by an integer. So we
shall get two LI Frobenius solutions.
Next equating to 0 the coefficient of xr+1 , we find
where n = 2, 3, 4....
For r = 1, we have
1 1 1 1 1
an = an−2 , a2 = a0 , a3 = a1 = 0, a4 = a2 = a0 , .......
n(2n + 3) 2.7 3.9 4.11 2.7.4.11
x2 x4
y 1 = a0 x 1 + + + ....... ,
2.7 2.7.4.11
x2 x4
−1/2
y 2 = a0 x 1+ + + ....... .
2.1 2.1.4.5
Ex. Find power series solutions of xy ′′ + y ′ − xy = 0 about x = 0.
Sol. Here x = 0 is a regular singular point of the given DE. So there exists at least one Frobenius
solution of the form
∞
X
y= an xn+r = xr (a0 + a1 x + a2 x2 + a3 x3 + ........). (5.12)
n=0
a0 r2 = 0 or r2 = 0.
Therefore, roots of the indicial equation are r = 0, 0, which are equal. So we shall get only one
Frobenious series solution.
Next equating to 0 the coefficient of xr , we find
a1 (r + 1)2 = 0 or a1 = 0 for r = 0.
10
Now equating to 0 the coefficient of xn+r−1 , we have the recurrence relation
1
an (n + r)2 − an−2 = 0 or an = an−2 .
(n + r)2
where n = 2, 3, 4....
Therefore, we have
1 1 1 1
a2 = a0 , a3 = a1 = 0, a4 = a2 = a0 , .......
(r + 2)2 (r + 3)2 (r + 4) 2 (r + 2) (r + 4)2
2
x2 x4
r
y = a0 x 1 + + + ....... (5.14)
(r + 2)2 (r + 2)2 (r + 4)2
Taking r = 0, we get the following Frobenius solution
x2 x4
y1 = a0 1 + 2 + 2 2 + .......
2 2 .4
Ex. Find power series solutions of x(1 + x)y ′′ + 3xy ′ + y = 0 about x = 0.
Sol. Here x = 0 is a regular singular point of the given DE. So there exists at least one Frobenius
solution of the form
∞
X
y= an xn+r = xr (a0 + a1 x + a2 x2 + a3 x3 + ........). (5.15)
n=0
a0 r(r − 1) = 0 or r(r − 1) = 0.
Therefore, roots of the indicial equation are r = 0, 1, which differ by an integer. So we shall get
only one Frobenious solution and that corresponds to the larger root r = 1.
Now equating to 0 the coefficient of xn+r−1 , we have the recurrence relation
n+r
an (n + r − 1) + an−1 (n + r) = 0 or an = − an−1 .
n+r−1
where n = 1, 2, 3, 4....
Therefore, we have
r+1 r+2 r+3
a1 = − a0 , a2 = − a0 , a3 = a0 , .......
r r r
For r = 1, we get a1 = −2a0 , a2 = 3a0 , a3 = −4a0 , ... So the Frobenious series solution is
11
5.3 Gauss’s Hypergeometric Equation
A DE of the form
a0 r(c + r − 1) = 0 or r(c + r − 1) = 0.
Therefore, roots of the indicial equation are r = 0, 1 − c. Now comparing the coefficient of xn+r−1 ,
we have the recurrence relation
(a + n − 1 + r)(b + n − 1 + r)
an (n+r)(c+n+r−1)−an−1 (n−1+r+a)(n−1+r+b) = 0 or an = an−1 .
(n + r)(c + n − 1 + r)
where n = 1, 2, 3, 4....
For r = 0, we have
(a + n − 1)(b + n − 1) a.b (a + 1)(b + 1) a(a + 1)b(b + 1)
an = an−1 , a1 = a0 a2 = a1 = a0 , .......
n(c + n − 1) 1.c 2(c + 1) 1.2c(c + 1)
So the Frobenius solution corresponding to r = 0 reads as
a.b a(a + 1)b(b + 1) 2
y = a0 1 + x+ x + .......... ,
1.c 1.2c(c + 1)
This series with a0 = 1 is called hypergeometric series and is denoted by F (a, b, c, x). Thus,
∞
X a(a + 1)...(a + n − 1)b(b + 1)...(b + n − 1)
F (a, b, c, x) = 1 + xn .
n=1
n!c(c + 1)...(c + n − 1)
F (1, b, b, x) = 1 + x + x2 + .........
the familiar geometric series. Thus, F (a, b, c, x) generalizes the geometric series. That is why it is
named as hypergeometric series. Further, we find
an+1 (a + n)(b + n)
lim |x| = lim |x| = |x| ,
n→∞ an n→∞ (n + 1)(c + n)
12
provided c is not zero or negative integer. Therefore, F (a, b, c, x) is analytic function-called the
hypergeometric function-on the interval |x| < 1. It is the simplest particular solution of the
hypergeometric equation.
Next we find the series solution corresponding to the indicial root r = 1 − c. The series solution
in this case is given by
where the constants a1 , a2 and so on can be determined using the recurrence relation. Alternatively,
we substitute y = x1−c z into the given DE (5.18) and obtain
z = F (a − c + 1, b − c + 1, 2 − c, x).
y = x1−c F (a − c + 1, b − c + 1, 2 − c, x).
where the prime denotes the derivative with respect to t. It is a hypergeometric equation with c
replaced by a + b − c + 1. So its solution with t replaced by 1 − x in view of (5.22) reads as
where
x−A
t=
B−A
13
and F , G, H are certain combinations of the constants in (5.23). The primes in (5.24) denote
the derivatives with respect to t. This is a hypergeometric equation with a, b and c defined by
F = −c, G = (a + b + 1) and H = ab. Therefore, (5.24) can be solved in terms of hypergeoetric
function near t = 0 and t = 1. It follows that (5.23) can be solved in terms of the same function
near x = A and x = B.
Remark: Most of the familiar functions in elementary analysis can be expressed in terms of
hypergeometric function.
(i) (1 + x)n = F (−n, b, b, −x)
(ii) log(1 + x) = xF (1, 1, 2, −x)
(iii) sin−1 x = xF (1/2, 1/2, 3/2, x2 )
(iv) ex = lim F (a, b, a, x/b)
b→∞
14
5.4 Misc. Practice Problems
1. Find power series solution of of the Legendre’s DE (1 − x2 )y ′′ − 2xy ′ + n(n + 1)y = 0 about
x = 0, where n is a constant.
Sol. We have
We observe that x = 0 is an ordinary point of (5.26). So there exists a series solution of the
form
∞
X
y= ak xk = a0 + a1 x + a2 x2 + a3 x3 + ......... (5.27)
k=0
(n − k + 2)(n + k − 1)
ak k(k − 1) + ak−2 (n − k + 2)(n + k − 1) = 0 or ak = − ak−2 .
k(k − 1)
∴ a2 = − n(n+1)
2!
a0 , a3 = (n−1)(n+2)
3!
a1 a4 = (n−2)n(n+1)(n+3)
4!
a0 , a5 = (n−3)(n−1)(n+2)(n+4)
5!
a1 , .......
Substituting these values into (5.27), we obtain the general solution of (5.26) as
y = c1 y 1 + c2 y 2
where
n(n + 1) 2 (n − 2)n(n + 1)(n + 3) 4
y1 = a0 1− x + x − ......... ,
2! 4!
(n − 1)(n + 2) 3 (n − 3)(n − 1)(n + 2)(n + 4) 5
y 2 = a1 x− x + x − ......... .
3! 5!
Note: Here y1 reduces to a polynomial of degree n when n is an even integer. Likewise,
y2 reduces to a polynomial of degree n when n is an odd integer. If we fix the values of a0
and a1 demanding the value of the polynomials to be 1 at x = 1, then such polynomials,
denoted by Pn (x), are called Legendre polynomials. We will study more details about
these polynomials in Chapter 8.
15
Corresponding to r = 2, we get
2 3 2 3 4
y = a0 x 1 − x + x − ............ .
10 56
It is the general solution. Notice that the solution corresponding to r = 2 is a particular case
of the general solution.
We see that xp(x) and x2 q(x) are analytic at x = 0. Therefore, x = 0 is a regular singular
point of the given differential equation, and hence there exists at least one Frobenius solution
of the form
∞
X
y= an xn+r = xr (a0 + a1 x + a2 x2 + a3 x3 + ........), (a0 ̸= 0). (1)
n=0
16
∞
X ∞
X
n+r
=⇒ an [2(n + r)(n + r − 1) − (n + r) + 1]x + an .2(n + r)xn+r+1 = 0.
n=0 n=0
∞
X ∞
X
n+r
=⇒ an (n + r − 1)(2n + 2r − 1)x + an .2(n + r)xn+r+1 = 0.
n=0 n=0
a0 (r − 1)(2r − 1) = 0 or (r − 1)(2r + 1) = 0.
Therefore, roots of the indicial equation are r = 1, 1/2 where the smaller root is r = 1/2.
Now equating to 0 the coefficient of xn+r , we have the recurrence relation
−2
an (n + r − 1)(2n + 2r − 1) + an−1 .2(n + r − 1) = 0 or an = an−1 ,
2n + 2r − 1
provided n + r − 1 ̸= 0.
For r = 1/2, we have
1 1 1 1 1
an = − an−1 , a1 = −a0 , a2 = − a1 = a0 , a3 = − a2 = − a0 , .......
n 2 2 3 3.2
Plugging these values into (1), we get
x2 x3
1/2 x
y 2 = a0 x 1− + − + ....... ,
1! 2! 3!
Python programs link: Some python programs for finding power series solutions are available
at the link: https://github.com/surgithub/M3/blob/master/chap5.ipynb
17