Reduction of Order
Reduction of Order
Reduction of Order
We shall take a brief break from developing the general theory for linear differential equations
to discuss one method (the reduction of order method) for finding the general solution to any
linear differential equation. In some ways, this method may remind you of the material in chapter
11. Indeed, part of the method involves solving a higher-order equation via first-order methods
as discussed in chapter 11. The general theory developed in chapter 12 will not, however, be
used to any great extent. Instead, the material developed here will help us finish that general
theory (at least partially confirming the suspicions raised at the end of the chapter), and will help
lead us to the complete result on constructing general solutions from particular solutions.
But why worry about completing that general theory if any linear differential equation can
be completely solved by this reduction of order method? Because this method requires that
one solution to the differential equation already be known. This limits the methods applicability.
Also, serious practical difficulties arise when the differential equation to be solved is of order
three or more. Still, there are situations where the method is of practical value, and it will help
us confirm suspicions we already have about general solutions.
Oh yes, there is another reason to develop this method: A rather powerful method for solving
nonhomogeneous equations, the variation of parameters method described in chapter 23, is
simply a clever refinement of the reduction of order method.
13.1
The reduction of order method is a method for converting any linear differential equation to
another linear differential equation of lower order, and then constructing the general solution
to the original differential equation using the general solution to the lower-order equation. In
general, to use this method with an N th -order linear differential equation
a0 y (N ) + a1 y (N 1) + + a N 2 y + a N 1 y + a N y = g
we need one known nontrivial solution y1 = y1 (x) to the corresponding homogeneous differential equation. We then try a substitution of the form
y = y1 u
where u = u(x) is a yet unknown function (and y1 = y1 (x) is the aforementioned known
solution). Plugging this substitution into the differential equation then leads to a linear differential
283
284
Reduction of Order
equation for u . As we will see, because y1 satisfies the homogeneous equation, the differential
equation for u ends up being of the form
A0 u (N ) + A1 u (N 1) + + A N 2 u + A N 1 u = g
remarkably, there is no A N u term. This means we can use the substitution
v = u
as discussed in chapter 11, to rewrite the differential equation for u as a (N 1)th -order
differential equation for v ,
A0 v (N 1) + A1 v (N 2) + + A N 2 v + A N 1 v = g
So we have reduced the order of the equation to be solved. If a general solution v = v(x)
for this equation can be found, then the most general formula for u can be obtained from v
by integration (since u = v ). Finally then, by going back to the original substitution formula
y = y1 u , we can obtain a general solution to the original differential equation.
This method is especially useful for solving second-order homogeneous linear differential
equations since (as we will see) it reduces the problem to one of solving relatively simple firstorder differential equations. Accordingly, we will first concentrate on its use in finding general
solutions to second-order, homogeneous linear differential equations. Then we will briefly
discuss using reduction of order with linear homogeneous equations of higher order, and with
nonhomogeneous linear equations.
13.2
The Method
Here we lay out the details of the reduction of order method for second-order homogeneous
linear differential equations. To illustrate the method, well use the differential equation
x 2 y 3x y + 4y = 0
Note that the first coefficient, x 2 , vanishes when x = 0 . From comments made in chapter 12
(see page 261), we should suspect that x = 0 ought not be in any interval of interest for this
equation. So we will be solving over the intervals (0, ) and (, 0) .
Before starting the reduction of order method, we need one nontrivial solution y1 to our
differential equation. Ways for finding that first solution will be discussed in later chapters. For
now let us just observe that if
y1 (x) = x 2 ,
then
x 2 y1 3x y1 + 4y1 = x 2
d 2
d 2 2
x 3x
x + 4 x2
2
dx
dx
= x 2 [2 1] 3x[2x] + 4x 2
= x 2 [2 (3 2) + 4] = 0
|
{z
}
0!
285
For our example (as already noted), we will seek a general solution to
x 2 y 3x y + 4y = 0
(13.1)
Let
y = y1 u
where u = u(x) is a function yet to be determined. To simplify the plugging into
the differential equation, go ahead and compute the corresponding formulas for the
derivatives y and y using the product rule:
y = (y1 u) = y1 u + y1 u
and
y =
=
=
=
y1 u + y1 u
y1 u + y1 u
y1 u + y1 u + y1 u + y1 u
= y1 u + 2y1 u + y1 u
y = (y ) = 2xu + x 2 u
= (2xu) + x 2 u
= 2u + 2xu + 2xu + x 2 u
= 2u + 4xu + x 2 u
2.
Plug the formulas just computed for y , y and y into the differential equation, group
together the coefficients for u and each of its derivatives, and simplify as far as possible.
(Well do this with the example first and then look at the general case.)
286
Reduction of Order
Plugging the formulas just computed above for y , y and y into equation
(13.1), we get
0 = x 2 y 3x y + 4y
= x 2 2u + 4xu + x 2 u 3x 2xu + x 2 u + 4 x 2 u
= 2x 2 u + 4x 3 u + x 4 u 6x 2 u 3x 3 u + 4x 2 u
= x 4 u + 4x 3 3x 3 u + 2x 2 6x 2 + 4x 2 u
= x 4 u + x 3 u + 0 u
Notice that the u term drops out! So the resulting differential equation for u
is simply
x 4 u + x 3 u = 0 ,
which we can further simplify by dividing out x 4 ,
u +
1
u = 0
x
In general, plugging in the formulas for y and its derivatives into the given differential
equation yields
0 = ay + by + cy
= a y1 u + 2y1 u + y1 u + b y1 u + y1 u + c y1 u
= a y1 u + 2a y1 u + ay1 u + by1 u + by1 u + cy1 u
= ay1 u + 2a y1 + by1 u + a y1 + by1 + cy1 u .
Au + Bu + Cu = 0
where
A = ay1
B = 2a y1 + by1
C = a y1 + by1 + cy1
and
Consequently,
C = a y1 + by1 + cy1 = 0
Now find the general solution to the second-order differential equation just obtained for
u,
Au + Bu = 0 ,
via the substitution method discussed in section 11.1:
(a)
287
dv
+ Bv = 0
dx
(It is worth noting that this first-order differential equation will be both linear and
separable.)
(b)
Find the general solution v(x) to this first-order equation. (Since it is both linear
and separable, you can solve it using either the procedure developed for first-order
linear equations or the approach developed for first-order separable equations.)
(c)
Using the formula just found for v , integrate the substitution formula u = v to
obtain the formula for u ,
Z
u(x) =
v(x) d x .
Dont forget all the arbitrary constants.
1
u = 0
x
Equivalently,
1
dv
= v
dx
x
1 dv
dx =
v dx
1
dx
x
ln |v| = ln |x| + c0
v = e ln|x| + c0
v = ec0 x 1
A
x
288
4.
Reduction of Order
Finally, plug the formula just obtained for u(x) into the first substitution,
y = y1 u
used to convert the original differential equation for y to a differential equation for u .
The resulting formula for y(x) will be a general solution for that original differential
equation. (Sometimes that formula can be simplified a little. Feel free to do so.)
In our example, the solution we started with was y1 (x) = x 2 . Combined with
the u(x) just found, we have
y = y1 u = x 2 [ A ln |x| + B] .
That is,
y(x) = Ax 2 ln |x| + Bx 2
and
y2 (x) = x 2 ln |x|
Since the A and B in the above formula for y(x) are arbitrary constants, and y2 is given
by that formula for y with A = 1 and B = 0 , it must be that this y2 is another particular
solution to our original homogeneous linear differential equation. Whats more, it is clearly not
a constant multiple of y1 . This should strengthen an earlier suspicion that the general solution
to a homogeneous linear second-order differential equation can be written as just such a linear
combination.
We will later examine more closely the general form for general solution to any homogeneous
linear differential equation. In the meantime, while practicing this method, do observe that the
general solution you obtain for each second-order homogeneous linear differential equation can,
invariably, be written as
y(x) = Ay2 (x) + By1 (x)
where y1 and y2 are two solutions that are not constant multiples of each other. Keep in mind
that this form is the same as the form earlier anticipated, namely,
y(x) = c1 y1 (x) + c2 y2 (x) .
Weve just renamed the arbitrary constants from c1 and c2 to B and A , respectively.
13.3
289
If you look back over our discussion in section 13.2, you will see that the reduction of order
method applies almost as well in solving a nonhomogeneous equation
ay + by + cy = g
Then, letting y = y1 u in the nonhomogeneous equation and then replacing u with v leads to
an equation of the form
Av + Bv = g
instead of
Av + Bv = 0
So we dont end up with a first-order equation for v which is both separable and linear; it is just
linear. Still, we know how to solve such equations. Solving that first-order linear differential
equation for v and continuing with the method already described finally yields the general
solution to the desired nonhomogeneous differential equation.
We will do one example. Then Ill tell you why the method is rarely used in practice.
!Example 13.1:
equation
(13.2)
= 2xu + x 2 u
and
y = (y ) = 2xu + x 2 u
= 2u + 2xu + 2xu + x 2 u
= 2u + 4xu + x 2 u
290
Reduction of Order
x = x 2 y 3x y + 4y
= x 2 2u + 4xu + x 2 u 3x 2xu + x 2 u + 4 x 2 u
= 2x 2 u + 4x 3 u + x 4 u 6x 2 u 3x 3 u + 4x 2 u
= x 4 u + 4x 3 3x 3 u + 2x 2 6x 2 + 4x 2 u
= x 4 u + x 3 u + 0 u
As before, the u term drops out. In this case, we are left with
x 4 u + x 3 u = x .
That is,
x 4v + x 3v = x
1/
2
with v = u
This is a relatively simple first-order linear equation. To help find the integrating factor, we
now divide through by x 4 , obtaining
dv
1
7
+ v = x /2
dx
x
1
x
dx
= eln|x| = |x|
Since we are just attempting to solve over the interval (0, ) , we really just have
.
= x
Multiplying the last differential equation for v and proceeding as usual when solving firstorder linear differential equations:
h 7 i
1
dv
x
+ v = x x /2
dx
dv
5
+ v = x /2
dx
d
5
xv = x /2
dx
d
xv d x =
dx
x /2 d x
2
3
2
3
xv = x /2 + c1
v = x /2 +
c1
x
Thus,
u =
du
dx =
dx
291
Z
2 5/2
c1
x
+
dx
3
2
4 3/2
x
+ c1 ln(x) + c2
9
2
3
x /2 + c1 ln(x) + c2
4 1/2
x + c1 x 2 ln(x) + c2 x 2
9
For no obvious reason at this point, lets observe that we can write this solution as
y(x) = c1 x 2 ln(x) + c2 x 2 +
4
x
9
(13.3)
It should be observed that, in the above example, we only used one particular solution,
y1 (x) = x 2 , to the homogeneous differential equation
x 2 y 3x y + 4y = 0
even though we had already found the general solution
Ax 2 ln |x| + Bx 2
Later, in chapter 23, we will develop a refinement of the reduction of order method for solving second-order nonhomogeneous linear differential equations which makes use of the entire
general solution to the corresponding homogeneous equation. This refinement (the variation
of parameters method) has two distinct advantages over the reduction of order method when
solving nonhomogeneous differential equations:
1.
The computations required for the refined procedure tend to be simpler and more easily
carried out.
2.
With a few straightforward modifications, the refined procedure readily extends to being
a useful method for dealing with nonhomogeneous linear differential equations of any
order. For the reasons discussed in the next section, the same cannot be said about the
basic reduction of order method.
That is why, in practice, the basic reduction of order method is rarely used with nonhomogeneous
equations.
292
13.4
Reduction of Order
In theory, reduction of order can be applied to any linear equation of any order, homogeneous or
not. Whether its application is useful is another issue.
!Example 13.2:
(13.4)
and think about what happens when you differentiate exponentials, you will realize that
y1 (x) = e2x
= 2e2x u + e2x u
and
y = e2x u
8e2x u + 12e2x u + 6e2x u + e2x u 8 e2x u = 0
e2x u + 6e2x u + 12e2x u + 8e2x 8e2x u = 0
Letting v = u and dividing out the exponential, this becomes the second-order differential
equation
v + 6v + 12v = 0 .
(13.5)
Additional Exercises
293
Thus we have changed the problem of solving the third-order differential equation to one
of solving a second-order differential equation. If we can now correctly guess a particular
solution v1 to that second-order differential equation, we could again use reduction of order
to get the general solution v(x) to that second-order equation, and then use that and the fact
that y = e2x u with v = u to obtain the general solution to our original differential equation.
Unfortunately, even though the order is less, guessing a solution to equation (13.5) is a good
deal more difficult than was guessing a particular solution to the original differential equation,
equation (13.4).
As the example illustrates, even if we can, somehow, obtain one particular solution to a given
N -order linear homogeneous linear differential equation, and then use it to reduce the problem
to solving an (N 1)th -order differential equation, that lower order differential equation may be
just as hard to solve as the original differential equation (unless N = 2 ). In fact, we will learn
how to solve differential equations such as equation (13.5), but those methods can also be used
to find the general solution to the original differential equation, equation (13.4), as well.
Still it does no harm to know that the problem of solving an N th -order linear homogeneous
linear differential equation can reduced to that of solving an (N 1)th -order differential equation,
especially since we may refer to this fact in the next chapter. For the record, here is a theorem to
that effect:
th
(13.6)
where g and the ak s are known functions on some interval of interest I , and let y1 be a
nontrivial particular solution to the corresponding homogeneous equation
a0 y (N ) + a1 y (N 1) + + a N 2 y + a N 1 y + a N y = 0
Set
u =
y
y1
(so that y = y1 u ) .
where the Ak s are functions on the interval I that can be determined from the ak s along with
y1 and its derivatives.
The proof is relatively straightforward: You see what happens when you repeatedly use the
product rule with y = y1 u , and plug the results into the equation (13.6). Ill leave the details
to you (see exercise 13.3).
294
Reduction of Order
Additional Exercises
13.1. For each of the following, first verify that the given y1 is a solution to the given
differential equation, and then find the general solution to the differential equation
using the given y1 with the method of reduction of order.
y1 (x) = e2x
a. y 5y + 6y = 0 ,
y1 (x) = e5x
b. y 10y + 25y = 0 ,
c. x 2 y 6x y + 12y = 0 ,
y1 (x) = x 3
d. 4x 2 y + y = 0
on x > 0 , y1 (x) = x
2
4
e. y 4 +
y + 4+
y = 0 , y1 (x) = e2x
x
f. y
1
y 4x 2 y = 0
x
g. y + y = 0 ,
y1 (x) = ex
y1 (x) = sin(x)
h. x y + (2 + 2x)y + 2y = 0 ,
y1 (x) = x 1
i. sin2 (x)y 2 cos(x) sin(x)y + 1 + cos2 (x) y = 0 ,
j. x 2 y 2x y + x 2 + 2 y = 0 , y1 (x) = x sin(x)
y1 (x) = sin(x)
13.2. Several nonhomogeneous differential equations are given below. For each, first verify
that the given y1 is a solution to the corresponding homogeneous differential equation,
and then find the general solution to the given nonhomogeneous differential equation
using reduction of order with the given y1 .
a. y 4y + 3y = 9e2x
y1 (x) = e3x
c. x 2 y + x y y = x , y1 (x) = x
d. x 2 y 20y = 27x 5
y1 (x) = x 5
e. x y + (2 + 2x)y + 2y = 8e2x
f. (x + 1)y + x y y = (x + 1)2
y1 = x 1
,
,
y1 = ex
b. N = 4
13.4. Each of the following is one of the relatively few third- and fourth-order differential
equations that can be easily solved via reduction of order. For each, first verify that the
Additional Exercises
295
given y1 is a solution to the given differential equation or to the corresponding homogeneous equation (as appropriate), and then find the general solution to the differential
equation using the given y1 with the method of reduction of order.
a. y 9y + 27y 27y = 0 ,
y1 = e3x
y1 = e3x
y1 = e2x
y1 = x 2