Math Methods 1
Lia Vas
Series Solutions of Ordinary Differential Equations
Let us consider a linear differential equation of the second order
a(x)y 00 + b(x)y 0 + c(x)y = g(x)
where a, b, c, and g are some real-valued functions.
Recall that the equation a(x)y 00 + b(x)y 0 + c(x)y = 0 is the homogeneous part of the above
relation. If yh is the solution of the homogeneous part of the non-homogeneous equation a(x)y 00 +
b(x)y 0 + c(x)y = g(x), the particular solution yp can be found using the variation of parameters
method and the solution of the homogeneous part. In this case, the general solution of the equation
has the form y = yh + yp .
Knowing the variation of parameters method, the problem boils down to finding the solution
of the homogeneous part. If the functions a(x), b(x) and c(x) are constant, methods of Differential
Equations course can be used. However, if these functions are not constant, another method is
needed.
So, we shall concentrate on methods for finding the solution of the homogeneous equation a(x)y 00 +
b(x)y 0 + c(x)y = 0 if the functions a(x), b(x) or c(x) are not constant.
There are certain methods for finding this solution in some particular cases. For example,
• Legendre’s linear equation a(px + q)2 y 00 + b(px + q)y 0 + cy = 0 where a, b, c, p and q are
constants. This equation can be reduced to an equation with constant coefficients using the
substitution px + q = et . If p = 1 and q = 0, the equation ax2 y 00 + bxy 0 + cy = 0 is called
Euler’s equation.
• Exact equations. If the left hand side of the equation a(x)y 00 + b(x)y 0 + c(x)y = 0 is a
derivative of another equation, it is said to be exact. It can be shown that an equation is
exact if c(x) − b0 (x) + a00 (x) = 0. Note that if a(x) = 0, we can write the equation in the form
b(x)dy = −c(x)ydx and denote Q = b and P = −cy, then the condition c(x) − b0 (x) = 0 is
equivalent to Py = Qx which is the criterion of exactness for the first order differential equation.
• Partially known complementary functions. If one solution y1 of the homogeneous equation a(x)y 00 +b(x)y 0 +c(x)y = 0 is known, to find the second solution y2 and the general solution
y = c1 y1 + c2 y2 as well, assume that it has the form y2 = vy1 where v is an unknown function.
Substituting the derivatives of y2 in the equation yields a separable differential equation in v 0 .
Solving for v 0 first and then finding v produces the second complementary function y2 of the
general solution.
However, none of the above methods are general. The goal of this section is to present a general
method that produces the solution of the equation a(x)y 00 + b(x)y 0 + c(x)y = 0 in a form of a power
series. Having a solution in the form of a power series should not be considered as a downside –
1
most commonly occurring functions such as polynomials, exponential and trigonometric functions
have power series extension.
Let us start by assuming that a(x) 6= 0 (so that the equation is really of the second order) and
by dividing by a(x) to obtain the form
y 00 + p(x)y 0 + q(x)y = 0
Although we shall denote the independent variable as x suggesting it is a real variable, our
following discussion and the methods presented are valid for functions of complex variables as
well.
Regular, Singular and Regular-Singular Points
Let us consider the equation y 00 + p(x)y 0 + q(x)y = 0 and let us assume that the function p(x)
and q(x) are analytic that is, that they have power series expansions
p(x) =
∞
X
n
pn (x − x0 ) and q(x) =
n=0
∞
X
qn (x − x0 )n
n=0
centered at a point x = x0 which are convergent on some interval containing x0 . In this case, the
point x = x0 is said to be a regular (or ordinary) point.
If the functions p(x) or q(x) are not analytic at x = x0 , then the point x = x0 is said to be a
singular point. In this case when it is possible to write the equation in the form
(x − x0 )2 y 00 + p(x)(x − x0 )y 0 + q(x)y = 0
where the functions p(x) and q(x) are analytic at x = x0 , the point x = x0 is said to be a regularsingular point. If it is not possible to represent the equation in this way, the point x = x0 is said
to be irregular or essential singularity.
The following examples showcase some equations that appear in applications in physics, and the
classifications of their relevant singularities.
Equation
Legendre
Chebyshev
Bessel
Laguerre
Hermite
singular points
type
(1 − x2 )y 00 − 2xy 0 + l(l + 1)y = 0
±1
regular-singular
(1 − x2 )y 00 − xy 0 + p2 y = 0
±1
regular-singular
x2 y 00 + xy 0 + (x2 − a2 )y = 0
0
regular-singular
xy 00 + (1 − x)y 0 + ay = 0
0
regular-singular
none
all pts regular
none
all pts regular
y 00 − 2xy 0 + 2ay = 0
Simple harmonic oscillator
y 00 + a2 y = 0
2
Solutions about an ordinary point
If x = x0 is a regular point of y 00 + p(x)y 0 + q(x)y = 0, then every solution of this equation can
be represented as a power series centered at x0 .
y=
∞
X
an (x − x0 )n
n=0
Note that in this case
y0 =
∞
X
nan (x − x0 )n−1 =
n=0
00
y =
∞
X
∞
X
(n + 1)an+1 (x − x0 )n and
n=0
n−2
n(n − 1)an (x − x0 )
=
n=0
∞
X
(n + 2)(n + 1)an+2 (x − x0 )n
n=0
The solution y can be obtained following the steps below.
1. Substitute y, y 0 and y 00 in the equation.
2. Write the left hand side of the equation as a single power series.
3. Equate the coefficients of the series you obtained to zero. This produces a recursive equation
that computes the terms of the sequence an .
4. If possible, obtain the sequence an as a function of n. Then obtain the solution y as a power
series. If you can express this series as an elementary function, then it is said that the solution
is in closed form. The close form will not always be possible to obtain.
When finding the closed form of the solutions, the power series expansions of the following
elementary functions may be useful:
∞
X
1 n
e =
x
n!
n=0
∞
X
1
=
xn
1 − x n=0
x
sin x =
∞
X
n=0
∞
X
(−1)n 2n
cos x =
x
(2n)!
n=0
(−1)n
x2n+1
(2n + 1)!
5. The radius of convergence of this power series is the distance from the center x = x0 to the
next nearest singularity of p and q.
1
To illustrate this last point, consider the equation y 00 − x−1
y = 0 about the regular point x = 0.
−1
1
For this equation p = 0 and q = x−1 = 1−x . Both of these functions are analytic: first because it
is a polynomial and the second since it is a rational function so both have continuous derivatives
of any order, p at any point and q at any point different than 1. Thus, the point x = 0 is regular.
Alternatively, the functions p and q are analytic since they have convergent power series expansions
1
at x = 0. The expansion of p is 0 + 0x + 0x2 + . . . . The expansion of q is 1−x
= 1 + x + x2 + . . . .
3
This series is convergent on (−1, 1) and so the power series solution of this equation is convergent
on (−1, 1) too. Note that the radius of convergence is 1 corresponds exactly the distance from the
center x = 0 to the first singularity x = 1.
To illustrate the method, we will start by an example that can be solved using much simpler
approach. Thus, we stress that this example is here just to illustrate how the method works, not to
illustrate the usefulness of the method.
Example 1 – “Fake” example. Consider the equation y 00 − 2y 0 + y = 0. Show that x = 0 is
a regular point of this equation. Then find the series solution at x = 0 and express your answers in
the closed form.
Solution. In this case p = −2 and q = 1. The functions p and q are analytic since the derivatives
of these functions (the constant functions 0) are continuous. Alternatively, you can say that p
and q are analytic since their power series expansions at x = 0 are p = 2 + 0x + 0x2 + . . . and
q = 1 + 0x + 0x2 + . . . which are convergent for every value of x. Thus, the point x = 0 is regular.
The solution of the equation can be found in the form
y=
∞
X
an x n .
n=0
Find the derivatives to be
0
y =
∞
X
nan x
n−1
=
y =
∞
X
n(n − 1)an x
(n + 1)an+1 xn and
n=0
n=0
00
∞
X
n−2
=
n=0
∞
X
(n + 2)(n + 1)an+2 xn
n=0
P∞
and substitute them into the equation to get n=0 (n + 2)(n + 1)an+2 xn − 2
P
∞
n
n=0 an x = 0 ⇒
∞
X
[(n + 2)(n + 1)an+2 − 2(n + 1)an+1 + an ]xn = 0
P∞
n=0 (n
+ 1)an+1 xn +
n=0
Since the series is identically equal to zero, all its terms have to be equal to zero. Thus,
(n + 2)(n + 1)an+2 − 2(n + 1)an+1 + an = 0 ⇒ an+2 =
2(n + 1)an+1 − an
for all n = 0, 1, . . .
(n + 1)(n + 2)
The above formula represents the recursive equation that computes the coefficients of the sequence
an . Note that this depends on two initial conditions a0 and a1 . Now you are looking for convenient
1 −a0
choices of a0 and a1 to obtain two solutions. In this case, the next term is a2 = 2a1·2
and with
a0
2a0 −a0
a0 = a1 this simplifies to a2 = 2 . Continuing computing the terms, you obtain a3 = 3·2 = a3!0 ,
a0
a4 = 2(3)a4·33 −a2 = 4·3·2
= a4!0 . . . an = an!0 . To simplify further assume that a0 = 1 and obtain the first
solution as
∞
X
1 2 1 3
1 n
y1 = 1 + x + x + x + . . . =
x = ex .
2!
3!
n!
n=0
Since you are looking for linearly independent second solution, make a fundamentally different
choice of a0 and a1 . For example, take a0 = 0. In this case a2 = 2a21 = a1 , a3 = 2(2)a3·21 −a1 = a21 ,
4
8a1
a
− 1
a1
1 −a1
4 −a3
a4 = 2(3)a4·33 −a2 = 3a4·3
= 2·3
= a3!1 , a5 = 8a5·4
= 6 5·4 2 =
Assuming that a1 = 1, you obtain the second solution as
5a1
5!
=
a1
4!
. . . an =
a1
(n−1)!
or an+1 =
a1
n!
∞
∞
X
X
1 3 1 4
1 n+1
1 n
y1 = x + x + x + x + . . . =
x
=x
x = xex .
2!
3!
n!
n!
n=0
n=0
2
Note that the radius of the convergence is infinity. We can make this conclusion either by noting
that the radius of convergence of both series representing two solutions is infinity or, easier, by noting
that the equation has no singular points.
Thus, the general solution is y = c1 ex + c2 xex .
As we pointed out, the equation y 00 − 2y 0 + y = 0 has constant coefficients so it could be solved by
considering the characteristic equation r2 − 2r + 1 = 0 ⇒ (r − 1)(r − 1) = 0 and finding the solution
in the form y = c1 ex + c2 xex . So, we turn to more relevant examples – examples of equations that
do not have constant coefficients.
Example 2 – “Real” Example. Consider the equation (1 − x)2 y 00 − 2y = 0. Show that x = 0
is a regular point of this equation. Then find the series solution at x = 0. Write your solution in the
closed form and determine the radius of the convergence of the solution.
−2
2
Solution. Dividing by (1 − x)2 we obtain the form y 00 − (1−x)
2 y = 0 so that p = 0 and q = (1−x)2 .
Both of these functions are analytic: first because it is a polynomial and the second since it is a
rational function so both have continuous derivatives of any order, p at any point and q at any
point different than 1. Thus, the point x = 0 is regular. Since 1 is a singularity of q, the radius of
convergence of the solutions is 1 (which is the distance from the center 0 to the singularity 1). So,
the interval of convergence of the solutions is (−1, 1).
Alternatively, you can argue that p and q are analytic since they have convergent power series
expansions at x = 0. The expansion of p is 0 + 0x + 0x2 + . . . . The expansion of q can be obtained
1
= −2(1+x+x2 +. . .). This series is convergent on (−1, 1) and so the expansion
differentiating −2 1−x
of q is convergent on (−1, 1). Thus, the solutions
on (−1, 1), too.
P are convergent
n
a
x
and
the
derivatives are
The solution can be found in the form y = ∞
n=0 n
0
y =
∞
X
nan x
∞
∞
∞
X
X
X
n
00
n−2
=
(n + 1)an+1 x and y =
n(n − 1)an x
=
(n + 2)(n + 1)an+2 xn
n−1
n=0
n=0
n=0
n=0
Substitute the function and its derivatives into the equation.
(1 − 2x + x2 )
∞
∞
X
X
(n + 2)(n + 1)an+2 xn − 2
an xn = 0 ⇒
n=0
∞
X
n=0
∞
∞
∞
X
X
X
n+1
n+2
(n + 2)(n + 1)an+2 x − 2
(n + 2)(n + 1)an+2 x +
(n + 2)(n + 1)an+2 x − 2
an x n = 0 ⇒
n
n=0
n=0
∞
X
(n + 2)(n + 1)an+2 xn − 2
n=0
n=0
∞
X
(n + 1)nan+1 xn +
n=0
∞
X
n=0
∞
X
n=0
n(n − 1)an xn − 2
∞
X
an x n = 0 ⇒
n=0
[(n + 2)(n + 1)an+2 − 2(n + 1)nan+1 + n(n − 1)an − 2an ]xn = 0.
n=0
5
Combine the terms with an and factor out n + 1 to get
∞
X
[(n + 2)(n + 1)an+2 − 2(n + 1)nan+1 + (n2 − n − 2)an ]xn = 0 ⇒
n=0
∞
X
[(n + 2)(n + 1)an+2 − 2(n + 1)nan+1 + (n + 1)(n − 2)an ]xn = 0 ⇒
n=0
∞
X
(n + 1)[(n + 2)an+2 − 2nan+1 + (n − 2)an ]xn = 0
n=0
Since the series is identically equal to zero, all its terms have to be equal to zero. Thus,
(n + 2)an+2 − 2nan+1 + (n − 2)an = 0 ⇒ an+2 =
2nan+1 − (n − 2)an
for all n = 0, 1, . . .
(n + 2)
3
Note that a2 = 2a20 = a0 , a3 = 2a03+a1 , a4 = 4a43 = a3 , a5 = 6a35−a3 = a3 , a6 = 8a3 −2a
= a3 . . .
6
Thus, choosing 2a0 + a1 = 0 would produce 0 = a3 = a4 = a5 = . . . . So the only non-zero coefficients
are a0 , a1 = −2a0 and a2 = a0 . This yields a polynomial solution. Choosing a0 = 1 we obtain the
first solution
y1 = 1 − 2x + x2 = (1 − x)2 .
Choosing a1 = a0 we obtain a2 = a0 , a3 = 2a03+a0 = a0 , a4 = a3 = a0 , . . . an = a0 . Choose a0 = 1
for simplicity again and obtain the second solution
2
3
4
y2 = 1 + x + x + x + x + . . . =
∞
X
n=0
xn =
1
.
1−x
Note that the radius of convergence of this series is 1 so it converges for −1 < x < 1 just as we
projected at the beginning.
The two solutions y1 and y2 are obviously independent (since one is a polynomial and the other
is not) so the general solution is
1
.
y = c1 (1 − x)2 + c2
1−x
Practice Problems.
1. Consider the differential equation y 00 − y = 0. Show that x = 0 is a regular point of this
equation. Find the series solutions of the given equation about x = 0. Express your answers in
closed form.
2. Consider Hermite equation y 00 −2xy 0 +4y = 0. Show that x = 0 is a regular point of this equation.
Find the series solutions of the given equation about x = 0. Find the closed form of one solution
and list first few terms of the second solution. Determine the interval of convergence.
Solutions.
6
1. For this equation p = 0 and q = −1 so both of these functions are analytic (all derivatives are
constant or, alternatively, the convergent power series expansions are p = 0 + 0x + 0x2 + . . .
and q = −1 +P
0x + 0x2 + . . .). Thus, x = 0 is a regular point and the solution can be found in
n
the form y = ∞
n=0 an x . Find the derivatives to be
0
y =
∞
X
nan x
n−1
∞
∞
∞
X
X
X
n
00
n−2
=
(n + 1)an+1 x and y =
n(n − 1)an x
=
(n + 2)(n + 1)an+2 xn
n=0
n=0
n=0
n=0
and substitute them into the equation.
∞
X
n
(n + 2)(n + 1)an+2 x −
∞
X
n=0
n
an x = 0 ⇒
n=0
∞
X
[(n + 2)(n + 1)an+2 − an ]xn = 0
n=0
Since the series is identically equal to zero, all its terms have to be equal to zero. Thus,
(n + 2)(n + 1)an+2 − an = 0 for all n = 0, 1, . . . ⇒ an+2 =
an
for all n = 0, 1, . . .
(n + 1)(n + 2)
The above formula represents the recursive equation that computes the coefficients of the
sequence an . In order to find the solution of the recursive equation, consider the first couple of
terms.
a2
a0
a0
a0
, a4 =
=
⇒ a2n =
a2 =
1·2
3·4
1·2·3·4
(2n)!
and
a3 =
a3
a1
a1
a1
, a5 =
=
⇒ a2n+1 =
2·3
4·5
2·3·4·5
(2n + 1)!
Recall that we are looking for two linearly independent solutions y1 and y2 to be produced by
these formulas. For example,
P
1 n
x
if we take a0 = a1 = 1, we obtain an = n!1 ⇒ y1 = ∞
n=0 n! x = e and
P
P
n
n
∞
(−x)n
If we take a0 = −a1 = 1, we obtain an = (−1)
⇒ y2 = n=0 (−1)
xn = ∞
= e−x .
n=0 n!
n!
n!
Alternatively, by taking a0 = 1 and a1 = 0, you obtain cosh(x) and by taking a0 = 0 and
a1 = 1, you obtain sinh(x).
Note that the radius of the convergence is infinity because the equation has no singular points.
Thus, the general solution is y = c1 ex + c2 e−x (alternatively y = c1 sinh(x) + c2 cosh(x)).
2. Here p = −2x and q = 4 are analytic since their derivatives are constant functions. Alternatively, you can deduce that p and q are analytic since their power series expansions are
p = 0 − 2x + 0x2 + . . . and q = 4 + 0x + 0x2 + . . . . So the point x = 0 is regular. Both
expansions for p and q are convergent for every x and so the series solution is convergent for
every x too.
The solution can be found in the form
y=
∞
X
n=0
n
0
an x ⇒ y =
∞
X
nan x
n−1
00
and y =
n=0
∞
X
n=0
7
n(n − 1)an x
n−2
=
∞
X
n=0
(n + 2)(n + 1)an+2 xn
Substitute these into the equation and get
∞
X
n
(n + 2)(n + 1)an+2 x − 2
n=0
∞
X
n
nan x + 4
n=0
∞
X
∞
X
an x n = 0 ⇒
n=0
[(n + 2)(n + 1)an+2 − 2nan + 4an ] xn = 0.
n=0
Since the series is identically equal to zero, all its terms have to be equal to zero. Thus,
(n + 2)(n + 1)an+2 − 2nan + 4an for all n = 0, 1, . . . ⇒ an+2 =
(2n − 4)an
for all n = 0, 1, . . .
(n + 1)(n + 2)
Thus the even-indexed coefficients depends on a0 and the odd-indexed coefficients on a1 . We
can obtain two linearly independent solutions by taking a0 = 1 and a1 = 0 for y1 and a0 = 0
and a1 = 1 for y2 .
In the first case, note that a2 =
polynomial y1 = 1 − 2x2 .
−4
2
= −2 and a4 = 0. Thus a6 = a8 = a10 = . . . = 0. So y1 is a
−2·2
In the second case, a3 = −2
, a5 = 2·3·4·5
= −2·2
, a7 = −2·2·4
. . . Thus y2 = x −
2·3
5!
7!
·2·4 7
2 7! x . . . . The general solution y = c1 y1 + c2 y2 is convergent on (−∞, ∞).
8
2 3
x
3!
−
2·2 5
x
5!
−
Solutions about a regular-singular point
Recall that x = x0 is a regular-singular point of the equation
(x − x0 )2 y 00 + p(x)(x − x0 )y 0 + q(x)y = 0
if the functions p(x) and q(x) are analytic at an interval containing x = x0 .
For simplicity, we can assume that x0 = 0. If that is not the case, the substitution t = x − x0 can
convert the equation with a regular-singular point x = x0 to an equation with regular-singular point
t = 0.
With the assumption that x0 = 0, we search for a solution in the form
y=x
r
∞
X
n
an x =
n=0
∞
X
an xn+r
n=0
for some number r. This series is called Frobenius series. We can assume that a0 6= 0 since
otherwise wePcan redefine a1 or some higher coefficient as a0 .
n+r
, we have that
For y = ∞
n=0 an x
y0 =
∞
∞
X
X
(n + r)(n + r − 1)an xn+r−2 .
(n + r)an xn+r−1 and y 00 =
n=0
n=0
The solution y can be obtained following the steps below.
1. Substitute y, y 0 and y 00 in the equation.
2. Write the left hand side of the equation as a single power series.
3. Consider the coefficient with the lowest power of x (equivalently, cancel the smallest power of x
and then plug x = 0 to obtain this coefficient). Equate this coefficient to zero. The quadratic
equation in r that you obtained in this way is called an indical equation. Let r1 and r2 denote
the two solutions of the indical equation, called the indices, and let r1 ≥ r2 . We distinguish
three relevant cases.
(a) The difference r1 −r2 is not an integer. In this case, the two linearly independent solutions
y1 and y2 are given by
y1 = x
r1
∞
X
an x
n
n=0
and y2 = x
r2
∞
X
bn x n .
n=0
(b) The difference r1 − r2 is a nonzero integer. If r1 is the larger root, the two linearly
independent solutions y1 and y2 are given by
y1 = x
r1
∞
X
an x
n
and y2 = cy1 ln x + x
n=0
r2
∞
X
n=0
9
bn x n .
(c) The difference r1 − r2 is zero. In this case, the two linearly independent solutions y1 and
y2 are given by
y1 = x
r1
∞
X
an x
n
and y2 = y1 ln x + x
r1 +1
∞
X
bn x n .
n=0
n=0
The coefficients an and bn can be determined from two recursive equations obtained
same
P∞ in the
r
n
way as in the case of a regular point – by substituting the solution y = x
n=0 an x in the
equation for two values r1 and r2 .
4. Just in the case of a regular point, if you can express these series as elementary functions, then
the solution is said to be in the closed form. The close form will not always be possible to
obtain. Also, the radius of convergence of this power series is the distance from the center
x = x0 to the next nearest singularity of p and q.
In the second or third case, the form of the second solution listed above can be obtained using a
method known as the derivative method. One of your project topics asks you to look into this
method in more details.
The following three examples illustrate these three cases.
Example 1 – Case 1 Example. Consider the equation 4xy 00 + 2y 0 + y = 0. Show that x = 0 is
a regular-singular point. Then find the closed form of the series solutions about x = 0 and determine
the interval of convergence.
1
1
and q = 4x
are not defined at x = 0.
Solution. The point x = 0 is not regular since p = 2x
However, x = 0 is a regular-singular point since the equation can be written as x2 y 00 + x2 y 0 + x4 y = 0
so p = 12 = 12 + 0x + 0x2 + . . . and q = x4 = 0 + 41 x + 0x2 + 0x3 + . . . are analytic at x = 0 and converge
for any point x. Thus,
y=x
r
∞
X
n
an x =
∞
X
an x
n+r
∞
∞
X
X
n+r−1
00
(n + r)(n + r − 1)an xn+r−2 .
(n + r)an x
and y =
⇒y =
0
n=0
n=0
n=0
n=0
Substitute into the equation to get
∞
X
n=0
4(n + r)(n + r − 1)an x
n+r−1
+
∞
X
n=0
2(n + r)an x
n+r−1
+
∞
X
an xn+r = 0.
n=0
The smallest power of x appears in the first two sums for n = 0 and it is xr−1 . The coefficient
with xr−1 is 4r(r − 1)a0 + 2ra0 . This gives you the indical equation
1
4r(r − 1)a0 + 2ra0 = 0 ⇒ 2a0 (2r2 − 2r + r) = 0 ⇒ 2r2 − r = 0 ⇒ r(2r − 1) = 0 ⇒ r = 0 or r = .
2
Thus, this is the first case since the difference r1 − r2 is not an integer.
P
n
Let us first consider the case r = 0 which will yield the first solution y = x0 ∞
n=0 an x from the
equation
10
∞
X
4n(n − 1)an xn +
n=0
∞
X
n[4(n − 1) + 2]an xn +
n=0
∞
X
2nan xn +
n=0
∞
X
∞
X
an xn+1 = 0 ⇒
n=0
an−1 xn = 0 ⇒ n(4(n − 1) + 2)an + an−1 = 0 ⇒ n(4n − 2)an = −an−1
n=1
−an−1
−an−1
= 2n(2n−1)
for n =
n(4n−2)
(−1)n
−1
. . . an = (2n)! . Thus, the
6!
This produces the recursive relation an =
produces a1 =
−1
,
2
a2 =
−a1
4·3
=
1
4·3·2
=
1
,
4!
a3 =
∞
X
(−1)n
y1 =
n=0
(2n)!
1, 2, . . . Choosing a0 = 1
first solution is
xn
Comparing this sum with the series expansion for cos x =
√
cos x.
(−1)n 2n
n=0 (2n)! x
P∞
we conclude that y1 =
P
n
Let us consider now the case r = 12 which yields the second solution y = x1/2 ∞
n=0 an x from the
equation
∞
∞
∞
X
X
X
1
1
1
4(n + )(n + − 1)an xn +
2(n + )an xn +
an xn+1 = 0 ⇒
2
2
2
n=0
n=0
n=0
∞
X
n
[(2n + 1)(2n − 1) + (2n + 1)]an x +
n=0
∞
X
an−1 xn = 0 ⇒ [(2n + 1)(2n − 1) + (2n + 1)]an + an−1 = 0 ⇒
n=1
(2n + 1)2nan = −an−1 ⇒ an =
Choosing a0 = 1 produces a1 =
the second solution is
−1
2·3
=
−1
,
3!
a2 =
−an−1
for n = 1, 2, . . .
2n(2n + 1)
−a1
5·4
=
1
5·4·3·2
=
1
,
5!
a3 =
−1
7!
. . . an =
(−1)n
.
(2n+1)!
Thus,
∞
X
(−1)n n
y2 = x
x
(2n + 1)!
n=0
P
(−1)n 2n+1
we conclude that y1 =
Comparing this sum with the series expansion for sin x = ∞
n=0 (2n+1)! x
P∞ (−1)n n+1/2 P∞ (−1)n 1/2 2n+1
√
= n=0 (2n+1)! (x )
= sin x.
n=0 (2n+1)! x
√
√
Thus, the general solution is y = c1 cos x+c2 sin x. The functions p and q converge on (−∞, ∞)
and so the interval of convergence of the solutions is also (−∞, ∞).
1/2
Example 2 – Case 2 Example. Consider the equation x2 y 00 + xy 0 − y = 0. Show that x = 0 is
a regular-singular point. Then find the closed form of the series solutions about x = 0.
Solution. Note that x = 0 is regular-singular point since p = 1 and q = −1 are analytic (all
derivatives are zero or, alternatively, constant functions have convergent power series expansion).
y = xr
∞
X
n=0
an x n =
∞
X
n=0
an xn+r ⇒ y 0 =
∞
∞
X
X
(n + r)an xn+r−1 and y 00 =
(n + r)(n + r − 1)an xn+r−2 .
n=0
n=0
P∞
n+r
into the
equation to get
+
n=0 (n + r)(n + r − 1)an x
P∞Substitute
P
∞
n+r
n+r
= 0 ⇒ n=0 [(n + r)(n + r − 1) + (n + r) − 1]an x
=0
n=0 an x
11
P∞
n=0 (n
+ r)an xn+r −
Equating the first term with zero produces the indical equation
r(r − 1) + r − 1 = 0 ⇒ r2 − 1 = 0 ⇒ (r − 1)(r + 1) = 0 ⇒ r1 = 1, r2 = −1
so we are dealing with the case when the difference of the roots is a non-negative integer.
Plugging that r = 1 in the equation gives us
x
∞
X
[(n + 1)n + (n + 1) − 1]an xn = 0 ⇒ [n2 + 2n]an = 0 ⇒ n(n + 2)an = 0, n = 0, 1, 2, . . .
n=0
Since the expression n(n + 2) is not zero for any positive value of n, we conclude that an = 0 for
all n = 1, 2, . . . . Thus, all the coefficients an are zero except possibly the first one, a0 . Taking a0 = 1
you obtain the first solution y1 = x(1 + 0) = x.
P
P∞
n
−1
n
The second
has the form y = cy1 ln x + xr2 ∞
n=0 bn x = cx ln x + x
n=0 bn x =
P∞ solution
n−1
cx ln x + n=0 bn x P. The derivatives are
P
n−3
n−2
.
y 0 = c + c ln x + ∞
and y 00 = xc + ∞
n=0 (n − 1)(n − 2)bn x
n=0 (n − 1)bn x
Plug the function and the derivatives into the equation and obtain
cx +
∞
X
(n − 1)(n − 2)bn xn−1 + cx + cx ln x +
n=0
∞
X
(n − 1)bn xn−1 − cx ln x −
n=0
∞
X
bn xn−1 = 0
n=0
P∞
The terms with ln x cancel and the remainingP
terms in 2cx+ n=0 [(n−1)(n−2)+(n−1)−1]bn xn−1
2
n−1
=0⇒
have to be equal to zero. Simplify to get 2cx + ∞
n=0 (n − 3n + 2 + n − 2)bn x
2cx +
∞
X
n(n − 2)bn xn−1 = 0 ⇒ 2cx − b1 + 3b3 x2 + 4(2)b4 x3 + 5(3)b5 x4 + . . . = 0
n=0
Considering the coefficients with each term we obtain that c = 0, b1 = 0 and b3 = b4 = b5 . . . = 0.
Thus, b0 and b2 are the only two possible nonzero coefficients. Thus, y2 = x1 (b0 + b2 x2 ) = bx0 + b2 x.
Since the last term is a constant multiple of the first solution, we can take b2 = 0. Taking b0 = 1 for
simplicity, we obtain the second solution y2 = x1 .
Thus, the general solution is y = c1 x + c2 x1 .
Example 3 – Case 3 Example. Consider the equation x2 y 00 − xy 0 + y = 0. Show that x = 0 is
a regular-singular point. Then find the closed form of the series solutions about x = 0.
Solution. Note that x = 0 is regular-singular point since p = −1 and q = 1 are analytic (all
derivatives are zero or, alternatively, constant functions have convergent power series expansion).
y = xr
∞
X
n=0
an x n =
∞
X
an xn+r ⇒ y 0 =
n=0
∞
∞
X
X
(n + r)an xn+r−1 and y 00 =
(n + r)(n + r − 1)an xn+r−2 .
n=0
n=0
P∞
P
n+r
n+r
into the equation to get
− ∞
+
n=0 (n + r)(n + r − 1)an x
n=0 (n + r)an x
P∞Substitute
n+r
a
x
=
0
⇒
n=0 n
∞
X
[(n + r)(n + r − 1) − (n + r) + 1]an xn+r = 0
n=0
12
Equate the coefficient of the smallest power of x with zero to produce the indical equation
r(r − 1) − r + 1 = 0 ⇒ (r2 − 2r + 1) = 0 ⇒ (r − 1)(r − 1) = 0 ⇒ r1 = r2 = 1
so we are dealing with the case when the difference of roots is zero.
Plugging that r = 1 in the equation gives us
x
∞
X
[(n + 1)n − (n + 1) + 1]an xn = 0 ⇒ [n2 + n − n − 1 + 1]an = 0 ⇒ n2 an = 0
n=0
Thus, all the coefficients an are zero except the first one, a0 when n = 0. Take a0 = 1 and obtain
the first solution y1 = x(1 + 0) = x.
P
P
n
n
=
= x ln x + x2 ∞
The P
second solution has the form y = y1 ln x + xr1 +1 ∞
n=0 bn x
n=0 bn x
∞
n+2
x ln x + n=0 bn x P. The derivatives are
P∞
1
n+1
00
n
y 0 = 1 + ln x + ∞
(n
+
2)b
x
and
y
=
+
n
n=0
n=0 (n + 2)(n + 1)bn x .
x
Plug the function and the derivatives into the equation and obtain
x+
∞
X
(n + 2)(n + 1)bn x
n+2
− x − x ln x −
∞
X
(n + 2)bn x
n+2
bn xn+2 = 0
n=0
n=0
n=0
+ x ln x +
∞
X
P∞
All the “non-series” terms cancel and the remaining series n=0 [(n+2)(n+1)bn −(n+2)bn +bn ]xn+2
has to have zero terms. Thus (n + 2)(n + 1)bn − (n + 2)bn + bn = 0 ⇒
(n2 + 3n + 2 − n − 2 + 1)bn = 0 ⇒ (n2 + 2n + 1)bn = 0 ⇒ (n + 1)2 bn = 0 ⇒ bn = 0
for all n. Thus, the second solution is y2 = x ln x and the general solution is y = c1 x + c2 x ln x.
The next example illustrate that in some special instances of the second case, the second solution
can beP
obtained from the first one. In other words, both solutions can be obtained in the form
∞
n
r
y=x
n=0 an x .
Example 4. – “Lucky instance of the second case” Example. Consider the equation
x y − 2xy 0 + 2y = 0. Show that x = 0 is a regular-singular point. Then find the closed form of the
series solutions about x = 0.
2 00
Solutions. For this equation, p = −2 and q = 2 and these two functions are analytic (all
derivatives are zero, alternatively, constant functions have convergent
power series expansion). So
P∞
x = 0 is a regular-singular point. Plugging the solution y = n=0 an xn+r and its derivatives yield
the equation
∞
X
[(n + r)(n + r − 1) − 2(n + r) + 2]an xn+r = 0
n=0
Equating the first term with zero produces the indical equation r(r − 1) − 2r + 2 = 0 ⇒ r2 − 3r + 2 =
0 ⇒ (r − 1)(r − 2) = 0. So, the difference r1 − r2 is an integer. The first solution can be obtained
by taking the larger of the two r-values, r = 2. In this case [(n + 2)(n + 1) − 2(n + 2) + 2]an = 0 ⇒
(n2 + 3n + 2 − 2n − 4 + 2)an = 0 ⇒ (n2 + n)an = 0 ⇒ n(n + 1)an = 0. In this case all coefficients an
are zero except possibly a0 . By taking a0 = 1, we obtain the solution y1 = x2 (1 + 0) = x2 .
Consider now r = 1. In this case (n + 1)nan − 2(n + 1)an + 2an = 0 ⇒ (n2 + n − 2n − 2 + 2)an = 0
⇒ (n2 − n)an = 0 ⇒ n(n − 1)an = 0. In this case all coefficients an are zero except possibly a0 (when
13
n = 0) and a1 (when n − 1 = 0). Thus, y2 = x(a0 + a1 x) = a0 x + a1 x2 . Since the second part has the
form of the first solution, you can take a1 = 0. Taking a0 = 1 for simplicity, you obtain the solution
y2 = x that is linearly independent from y1 = x2 .
Thus, the general solution is y = c1 x2 + c2 x.
Practice Problems.
1. Consider the equation 3x2 y 00 − 4xy 0 + 2y = 0. Show that x = 0 is a regular-singular point.
Then find the closed form of the series solutions about x = 0 and determine the interval of
convergence.
2. Consider the equation xy 00 − xy 0 + y = 0. Show that x = 0 is a regular-singular point. Then
find the series solutions about x = 0. Find the closed form of one solution and list first few
terms of the second solution.
3. Consider the equation x2 y 00 + 2xy 0 − x2 y = 0. Show that x = 0 is a regular-singular point. Then
find the series solutions about x = 0. Find the series form of one solution and write down the
form of the second solution (do not need to solve for coefficients of the second solution).
4. Consider the equation x(1 − x)y 00 + (1 − x)y 0 + y = 0. Show that x = 0 is a regular-singular
point. Then find the series solutions about x = 0. Find the closed form of one solution and
write down the form of the second solution (do not need to solve for coefficients of the second
solution). Determine the interval of convergence.
Solutions.
1. For this equation, p = − 43 and q = 32 . Constant functions have power series expansions which
converges at every point. So x = 0 is a regular-singular point and the series solution is
convergent on (−∞, ∞).
P
n+r
and its derivatives yield the equation
Plugging the solution y = ∞
n=0 an x
∞
X
3(n + r)(n + r − 1)an x
n=0
n+r
−
∞
X
4(n + r)an x
n+r
+
∞
X
2an xn+r = 0
n=0
n=0
r
Equating the coefficient of the smallest power of x, x in this case, with zero produces the
indical equation 3r(r − 1) − 4r + 2 = 0 ⇒ 3r2 − 7r + 2 = 0 ⇒ r1 = 2 and r2 = 13 . So, the
difference r1 − r2 is not an integer.
When r = 2 the equation becomes
∞
X
[3(n + 2)(n + 1)an − 4(n + 2)an + 2an ]xn+2 = 0 ⇒ [3(n + 2)(n + 1) − 4(n + 2) + 2]an = 0 ⇒
n=0
(3n + 5)nan = 0 ⇒ an = 0 for all n > 0. By taking a0 = 1 we obtain the solution y1 =
x2 (1 + 0) = x2 .
When r =
1
3
the equation becomes
∞
X
1
2
1
1
2
1
[3(n + )(n − )an − 4(n + )an + 2an ]xn+2 = 0 ⇒ [3(n + )(n − ) − 4(n + ) + 2]an = 0 ⇒
3
3
3
3
3
3
n=0
14
(3n − 5)nan = 0 ⇒ an = 0 for all n > 0. By taking a0 = 1 we obtain the solution y2 =
x1/3 (1 + 0) = x1/3 .
So, the general solution is y = c1 x2 + c2 x1/3 .
2. For this equation, p = −x and q = x are analytic (power series expansions are p = 0 − 1x +
are convergent for any x). So x = 0 is
0x2 + 0x3 + . . . and q = 0 + 1x + 0x2 + 0x3 + . . . and
P∞
a regular-singular point. Plugging the solution y = n=0 an xn+r and its derivatives yield the
equation
∞
X
(n + r)(n + r − 1)an xn+r−1 −
n=0
∞
X
(n + r)an xn+r +
n=0
∞
X
an xn+r = 0
n=0
The smallest power of x is xr−1 . Its coefficient produces the indical equation r(r − 1) = 0 ⇒
r1 = 1 and r2 = 0. So, the difference r1 − r2 is an integer.
One solution can be obtained by considering r = 1. In this case
∞
X
n
n(n + 1)an x −
n=0
∞
X
(n + 1)an x
n+1
+
n=0
∞
X
∞
X
an xn+1 = 0 ⇒
n=0
[(n + 1)(n + 2)an+1 − (n + 1)an + an ]xn+1 = 0 ⇒
n=0
nan
for n = 1, 2 . . . Thus a1 = 0 so a2 = a3 = . . . = 0. Taking a0 = 1
2a1 = 0 and an+1 = (n+2)(n+1)
we obtain the first solution y1 = x(1 + 0) = x.
P
n
The other solution can be found in the form y = cx ln x + ∞
n=0 bn x . Plugging this and its
derivatives in the equation produces
c+
∞
X
n(n − 1)bn x
n−1
∞
X
− cx ln x − cx −
n
nbn x + cx ln x +
n=0
n=0
c − cx +
∞
X
∞
X
bn x n = 0 ⇒
n=0
[(n + 1)nbn+1 − nbn + bn ]xn = 0 ⇒
n=0
Equate the left-hand side terms with xn for any n with zero.
For n = 0, we have that c + b0 = 0. So, b0 = −c.
For n = 1, we have that −c + 2b2 = 0. So, b2 =
c
2
For n = 2, 3, . . . , we have that [(n + 1)nbn+1 − nbn + bn ⇒ bn+1 =
(n−1)bn
.
n(n+1)
b2
c
c
c
c
c
c
3
4
Thus b3 = 3·2
= 3·2·2
= 2·3!
, b4 = 2b
= 4·3·3·2
= 3·4!
, b5 = 3b
= 5·4·4·3·2
= 4·5!
. . . and
4·3
5·4
c
so bn = (n−1)n! for n = 2, 3, . . . . We can take b1 = 0 and c = 1 for simplicity. Thus y2 =
P
xn
x ln x − 1 + ∞
n=2 (n−1)·n! and the general solution is y = c1 x + c2 y2 .
3. For this equation, p = 2 and q = −x2 . These functions are analytic with convergent power
series expansions are p = 2 + 0x + 0x2 + . . . and q = 0 + 0x − 1x2 + 0x3 + 0x4 . . .. So x = 0 is
15
a regular-singular point. Plugging the solution y =
equation
∞
X
(n + r)(n + r − 1)an xn+r + 2
n=0
∞
X
P∞
n=0
an xn+r and its derivatives yield the
(n + r)an xn+r −
n=0
∞
X
an xn+r+2 = 0
n=0
Equating the coefficient of the smallest power of x (the first term in the first two sums in this
case) with zero produces the indical equation r(r − 1) + 2r = 0 ⇒ r(r + 1) = 0 ⇒ r1 = 0 and
r2 = −1 So, the difference r1 − r2 is an integer.
One solution can be obtained by considering r = 0. In this case
P∞
P∞
P∞
n
n
n+2
=0⇒
n=0 n(n − 1)an x + 2
n=0 nan x −
n=0 an x
P∞
2a1 x + n=0 [(n + 2)(n + 1)an+2 + 2(n + 2)an+2 − an ]xn+2 = 0 ⇒
an
2a1 = 0 and an+2 = (n+2)(n+3)
for n = 0, 1, . . . Thus a1 = 0 and so a3 = a5 = . . . = 0. The even
a0
a0
a0
. Taking a0 = 1 we obtain the first solution
terms are a2 = 3·2 = 3! , a4 = a5!0 , . . . a2n = (2n+1)!
P∞
P
x2n
n
y1 = n=0 (2n+1)! . The second solution has the form y2 = cx ln x + x−1 ∞
n=0 bn x . The general
solution is y = c1 y1 + c2 y2 .
x
4. Multiplying the equation by x and dividing by 1 − x we obtain the form x2 y 00 + xy 0 + 1−x
y = 0.
x
For this form, we can see that p = 1 and q = 1−x . The function p is analytic with the power
1
series 1 + 0x + 0x2 + . . . convergent for any n. The function q = x 1−x
= x(1 + x + x2 + . . .) =
1
x + x2 + x3 + . . . is convergent on interval (−1, 1) since the expansion 1 + x + x2 + . . . of 1−x
is convergent on (−1, 1). So, x = 0 is a regular-singular point and the solutions are convergent
on interval (−1, 1). You can reach the same conclusion regarding the interval of convergence by
noting that x = 1 is a singularity of q so the radius of convergence is 1 (which is the distance
from the center 0 to the singularity 1) and, hence, the interval of convergence is (−1, 1).
P
n+r
and its derivatives into the equation xy 00 − x2 y 00 + y 0 −
Plugging the solution y = ∞
n=0 an x
0
xy + y = 0 produces
∞
X
(n + r)(n + r − 1)an x
n+r−1
n=0
−
∞
X
(n + r)(n + r − 1)an x
n+r
n=0
∞
X
+
∞
X
(n + r)an xn+r−1 −
n=0
(n + r)an x
n+r
n=0
+
∞
X
an xn+r = 0
n=0
Equating the coefficient of the smallest power of x, xr−1 in this case, with zero produces the
indical equation r(r − 1) + r = 0 ⇒ r2 = 0 ⇒ r1 = r2 = 0. So, this is the third case.
P∞
P∞
P∞
n−1
n
n−1
For
r
=
0,
the
equation
becomes
n(n
−
1)a
x
−
n(n
−
1)a
x
+
−
n
n
n=0
n=0
n=0 nan x
P∞
P
P
∞
∞
n
n
n
na
x
+
a
x
=
0
⇒
[(n
+
1)na
−
n(n
−
1)a
+
(n
+
1)a
−
na
+
a
]x
=
n
n+1
n
n+1
n
n
n=0
n=0 n
n=0
(n2 −1)an
(n−1)an
0 ⇒ an+1 = (n+1)2 = n+1 ⇒ a1 = −a0 , a2 = 0 ⇒ a3 = a4 = a5 = . . . 0. With a0 = 1, we
obtain y1 = 1 − x.
P
n
The second solution has the form y2 = (1 − x) ln x + x ∞
n=0 bn x and the general solution is
y = c1 (1 − x) + c2 y2 .
16