Advanced Engineering Mathematics: Power Series and Power Series Solutions of Differential Equations
Advanced Engineering Mathematics: Power Series and Power Series Solutions of Differential Equations
a , a , a , ⋯ , a
in a given order. Each of a , a , a and so on represents a number. These are the terms of the sequence.
An infinite sequence of numbers is a function whose domain is the set of positive integers. For example,
the function associated to the sequence
sends 1 to a = 2, 2 to a = 4 and so on. The general behavior of this sequence is described by the
formula
a = 2n
Sequences can be described by writing rules that specify their terms or by listing terms.
Convergence and Divergence of Sequences. Sometimes, the numbers in a sequence approach a single
value as the index n increases. This happens in the sequence
1 1 1 1
1, , , , ⋯ , , ⋯
2 3 4 n
1 2 3 4 1
0, , , , , ⋯ ,1 − , ⋯
2 3 4 5 n
whose terms approach one. Such sequences are said to converge to that value. On the other hand,
sequences like
have terms that get larger than any number as n increases; such sequences are said to diverge. Still there
are sequences like
that terms oscillate between two values and does not converge to a single value. The sequence a
converges to the number L if to every positive number ϵ there corresponds an integer N such that for all n,
n > N and |a − L| < ϵ , where L is the limit of the sequence a . If no such L exists, then the
sequence is said to diverge.
A series (or for this matter, an infinite series) is the sum of an infinite sequence of numbers
a + a + a + ⋯ + a + ⋯
s = a + a + a + ⋯ + a
is an ordinary finite sum and is called the nth partial sum. As n gets larger, we expect the partial sums to
get closer and closer to a limiting value in the same sense that the terms of a sequence approach a limit.
For example, the series
1 1 1 1
1 + + + + +⋯
2 4 8 16
s = 1
1 3
s = 1 + =
2 2
1 1 7
s = 1 + + =
2 4 4
1 1 1 15
s' = 1 + + + =
2 4 8 8
⋮
1 1 1 1 1
s = 1 + + + + ⋯ + ) = 2 − )
2 4 8 2 2
The partial sums form a sequence whose nth term is
1
s = 2 −
2)
This sequence of partial sums converges to 2, because
1
lim =0
→. 2)
Thus the sum of the sequence /1 + + ' + 0 + 1 + ⋯ 2 is 2.In general, the series
a + a + a + ⋯ + a + ⋯
has partial sums that form the terms of the sequences . If the latter sequence converges to a limit L,
then L is the sum of the series. If the sequence of the partial sums of the series does not converge, then
the series diverges.
Example 4.1
1. Find the sum of the geometric series a + ar + ar + ar + ⋯ + ar ) + ⋯ = ∑. 5 ar
)
.
Determine the necessary conditions for this series to converge and the value to which this series
converges.
2. Using the results above, determine the sum of a geometric series whose first term is 1/9 and a common
ratio 1/3.
)8 9
3. Using again the results of the first problem, determine the sum ∑.
5: '8
.
Answers:
;)<8
1. The nth partial sum is s = )< , r ≠ 1. The series converges when |r| < 1 and this converges to
;
s = )<, otherwise, the series diverges.
2. s = 1
3. s = 4
4. ∑. = 1, ∑.
1
5 5 ) = 3
? n = 1 + 4 + 9 + ⋯ + n + ⋯
5
diverges because the partial sums grow beyond every number L. Thus we can observe that if the partial
sum s = ∑. 5 a converges, then a → 0. Thus, we can test if the series diverges based on the above
observation: that the partial sum s = ∑.
5 a diverges if lim→. a fails to exist or is different from
zero.
Example 4.2
Determine the following series if it diverges using the nth term test for divergence.
1. ∑.
5 n
2. ∑.
5
3. ∑.
5−1
Answers:
1. The series diverges because n → ∞.
→1
2. The series diverges because
3. The series diverges because lim→. −1 does not exist.
?a + b = ? a + ? b = A + B (4.1a)
?a − b = ? a − ? b = A − B (4.1b)
? ka = k ? a = kA (4.1c)
is divergent, even though → 0. Thus the nth term test for divergence fails in this case. The reason why it
diverges is because there is no upper bound for its partial sums.
If we let the sequence as a . Suppose that a = fn, where f is a continuous, positive, decreasing
.
function of x for all x ≥ N (N is a positive integer. Then the series ∑.
5H a and the integral IH fx dx
both converge or both diverge.
Example 4.3
Test the following series for convergence using the integral test.
1. ∑.
5
2. ∑.
5 K
Answers:
1. The series diverges.
2. The series converges.
The ratio test. Let ∑ a be a series with positive terms and suppose that
a
lim =p
→. a
(4.2)
then,
Example 4.4
Investigate the convergence of the following series, and if the series is convergent, find its sum.
8 9
1. ∑.
5: 8
!
2. ∑.
5 !!
'8 !!
3. ∑.
5 !
Answers:
1. The series converges. s = .
2. The series diverges.
3. The series diverges.
The root test. Let ∑ a be a series with a ≥ 0 for n ≥ N and suppose that
lim 8Na = p
→.
(4.3)
then,
Example 4.5
Investigate the convergence of the following series using root test.
K
1. ∑.
5 8
8
2. ∑.
5 K
3. ∑.
5 OP
Answers:
1. The series converges.
2. The series diverges.
3. The series converges.
Alternating series, absolute and conditional convergence. A series in which the terms are alternately
positive and negative is an alternating series. For example, the alternating harmonic series
1 1 1 1 −1
1 − + − + − ⋯+ +⋯
2 3 4 5 n
is convergent, as well as the alternating geometric series
1 1 1 −1 4
−2 + 1 − + − + ⋯ + +⋯
2 4 8 2
1 − 2 + 3 − 4 + 5 − 6 + ⋯ + −1 n + ⋯
however, is divergent because the nth term does not approach zero. In general, the alternating series
3. u → 0.
A series ∑ a converges absolutely (is absolutely convergent) if the corresponding series of absolute
values, ∑|a |, converges. In effect, a series that converge but does not converge absolutely converges
conditionally.
The absolute convergence test. If ∑|a | converges, then ∑ a converges, although the converse of this
statement is generally not true (take a look at the harmonic series).
Example 4.6
Investigate the convergence of the following series.
1. ∑.
5−1
K
) 8RS
2. ∑.
5 T
for p > 0, p > 1, and 0 < p ≤ 1.
Answers:
1. The series converges absolutely.
2. For p > 0, the series converges. For p > 1, the series converges absolutely. For 0 < p ≤ 1, the
series converges conditionally.
W
3. a = O1 + P
5. ∑.
'
5 ')'
7. ∑.
X
5
8. ∑.
5 K
9. ∑.
5 O1 −
P
10. ∑.
5 O − K P
11. ∑.
! X
5 !
12. ∑.
!
5
13. ∑.
5−1
O:P
14. ∑.
5−1
ln O1 + P
15. ∑.
5 n e
)
? c x = c: + c x + c x + ⋯ + c x + ⋯ (4.4)
5:
in which the center a and the coefficients c: , c , c , ⋯ , c , ⋯ are constants. For 4.4, taking all the
coefficients to be 1 gives the geometric power series
? x = 1 + x + x + ⋯ + x + ⋯
5:
This series converges to )[ for |x| < 1. We express this fact by writing
1
= 1 + x + x + ⋯ + x + ⋯
1−x (4.6)
for −1 < x < 1. Up to now, we have used the equation 4.6 as a formula for the sum of the series on the
right. At this point, we wish to emphasize the fact that the polynomial on the right P x approximates the
function on the left.
1 1 1
1 − x − 2 + x − 2 + ⋯ + ]− ^ x − 2 + ⋯
2 4 2
2 1 1 1
= 1 − x − 2 + x − 2 + ⋯ + ]− ^ x − 2 + ⋯
x 2 4 2 (4.7)
? a x = a: + a x + a x + ⋯ + a x + ⋯ (4.8)
5:
converges for x = c ≠ 0, then it converges absolutely for all x with |x| < c. If the series diverges for
x = d, then it diverges for all x with |x| > d.
Radius of convergence of power series. The convergence of the series ∑ c x − a is described by
one of the following three possibilities:
1. There is a positive number R such that the series diverges for x with |x − a| > R but converges
absolutely for x with |x − a| < R. The series may or may not converge at either of the endpoints x = a −
R and x = a + R.
The number R is called the radius of convergence of the power series and the interval of radius R centered
at x = a is called the interval of convergence. The interval of convergence may be open, closed, or half-
open, depending on the particular series.
1. Use the ratio test (or the nth root test) to find the interval where the series converges absolutely.
Ordinarily, this is an open interval |x − a| < R or a − R < x < a + R.
2. If the interval of absolute convergence is finite, test for convergence at each end point. Use comparison
test, integral test, or the alternating series test.
3. If the interval of absolute convergence is a − R < x < a + R, the series diverges for |x − a| > R (it
does not even converge conditionally), because the nth term does not approach zero for those values of x.
Example 4.7
For what values of x do the following power series converge?
8
1. ∑.
5−1
) [
[K8RS
2. ∑.
5−1
)
)
[8
3. ∑.
5: !
4. ∑.
5: n! x
Answers:
1. −1 < x ≤ 1
2. −1 ≤ x ≤ 1
3. All x
4. All x except x = 0
Term-by-term differentiation. If ∑ c x − a converges for a − R < x < a + R for some R > 0, it
defines a function f
.
fx = ? c x − a
(4.9a)
5:
Such a function f has derivatives of all orders inside the interval of convergence. We can obtain the
derivatives by differentiating the original series term-by-term.
.
and so on. Each of these derived series converges at every interior point of the interval of convergence of
the original series.
Example 4.8
Obtain the first and second derivative of the function
.
1
fx = = 1 + x + x + x + ⋯ + x + ⋯ = ? x
1−x
5:
Answers:
.
1
f x =
a
= 1 + 2x + 3x + 4x + ⋯ + nx ) + ⋯ = ? nx )
1 − x
5
.
aa
2
f x = = 2 + 6x + 12x + ⋯ + nn − 1x ) + ⋯ = ? nn − 1x )
1 − x
5
fx = ? c x − a
(4.9a)
5:
with
.
x − a
? c
n+1
5:
Example 4.9
Investigate the following series
1.
x x9
fx = x − + − ⋯
5 5
for −1 ≤ x ≤ 1.
2.
ft = 1 − t + t − t + ⋯
Answer:
1. Using term-by-term differentiation and integration of the resulting series, we find that
x x9
fx = x − + − ⋯ = tan) x
5 5
x x x'
ln1 + x = x − + − +⋯
2 3 4
. . .
e? a x f e?
b x f
= ? c x (4.11)
5: 5: 5:
1. ∑.
5:x + 5
2. ∑.
5:−1 4x + 1
[)8
3. ∑.
5
[8
4. ∑.
5:
) 8 [8
5. ∑.
5
B. For the following series, find the series’ interval of convergence and within this interval, the sum of the
series as a function of x.
[)K8
6. ∑.
5 '
[K8
7. ∑.
5: g8
8. ∑.
√[
5: O − 1P
9. ∑.
5:ln x
[K
10. ∑.
5: O
P
This section shows how functions that are infinitely differentiable generate power series called Taylor
series. This is when we think that a function fx has derivatives of all orders on an interval I can be
expressed as a power series on that interval.
If we have
.
f a a = c
f aa a = 1 ∙ 2c
f a = n! c
With this, a general formula for the nth term of the coefficients of the power series can be deduced, that is
f a
c =
n!
(4.12)
which is the Taylor series expansion of fx about x = a. When a = 0, the series becomes
.
f d 0 d f aa 0 f 0
fx = ? x = f0 + f a 0x + x + ⋯+ x +⋯
k! 2! n!
(4.14)
d5:
The right hand side of equations 4.13 and 4.14 are called the Taylor (Maclaurin) polynomials of order n of
the function fx.
These polynomials provided a linear approximation of the function fx in the neighborhood of a (or 0)
Example 4.10
1. Find the Taylor series generated by fx = [ at a = 2. Where, if anywhere, does the series converge to
[
.
2. Find the Taylor series expansion, and the Taylor polynomials generated by fx = e[ at x = 0.
3. Find the Taylor series and Taylor polynomials generated by fx = cos x at x = 0.
Answer:
[) [)K [) 8
1. fx = − K + m − ⋯ + −1 8nS + ⋯. The series converges for 0 < x < 4.
[o
2. The Taylor series expansion of fx is fx = ∑.
d5: d! . The Taylor polynomials generated by the
function is
x x
P x = 1 + x + + ⋯
2 n!
) o [Ko
3. The Taylor series generated by function f at x = 0 is ∑.
d5: d!
. The 2nth order (as well as the
2n + 1th order) polynomial of cos x is then
x x' x
P x = P x = 1 − + − ⋯ + −1
2! 4! 2n!
2. fx = e[/
3. fx = [
4. fx = )[
5. sin 3x
6. cosh x
7. sinh x
8. x ' − 2x − 5x + 4
9. x + 1
[
10. sin
12. fx = 2x + 2x + 3x − 8, a = 1
15. fx =
[K
,a=1
17. fx = e[ , a = 2
18. fx = 2[ , a = 1
C. Find the first and second order Taylor polynomial approximation of the following functions at x = 0.
19. fx = lncos x
4.4 Power Series Solution of Differential Equations: Solutions Near an Ordinary Point
the point x = x: is called an ordinary point of 4.17 if b: x: ≠ 0. A singular point of the linear equation is
any point x = x for which b: x = 0. Any point that is not a singular point is an ordinary point.
1 − x y aa − 6xy a − 4y = 0
has x = 1 and x = −1 as its only singular points in the finite complex plane. The equation
y aa + 2xy a + y = 0
xy aa + y a + xy = 0
has the origin, x = 0, as the only singular point in the finite plane.
The solution of linear equations with constant coefficients can be accomplished by methods developed
earlier. For linear, ordinary differential equations with variable coefficients, and of order greater than one,
probably the most generally effective method of attack is that based upon the use of power series.
In solving differential equations near an ordinary point, we assume a solution of a power series of the form
.
y = ? c x (4.18)
5:
and substitute this to the given differential equation. The task then is to find the general term of the series
c . The method described is demonstrated in the following examples.
Example 4.11
1. Solve the equation y aa + 4y = 0 near the ordinary point x = 0.
Answers:
1. y = c: cos 2x + c sin 2x.
s sS u[)[m v
2. y = )[tK K + )[K K
2. 1 − 4x y aa + 8y = 0
3. 1 + x y aa + 10xy a + 20y = 0
4. x − 9y aa + 3xy a − 3y = 0
5. x + 4y aa + 6xy a + 4y = 0
6. y aa + x y = 0
7. 1 + 2x y aa + 3xy a − 3y = 0
8. y aa + xy a + 3y = x
9. y aa + 3xy a + 7y = 0
10. x + 4y aa + xy a − 9y = 0