7.4: Variation of Parameters for Higher Order Equations
( \newcommand{\kernel}{\mathrm{null}\,}\)
Derivation of the Method
If you have not had Math 410 (linear algebra), then you may want to review Appendix 11.4 again before starting this section.
Although the method of variation of parameters for higher order equations is similar to what we did in Section 5.6 for second order equations, there are differences that we'll need to discuss here. For convenience, we consider linear differential equations written as
P0(x)y(n)+P1(x)y(n−1)+⋯+Pn(x)y=F(x),
which can be rewritten as
y(n)+p1(x)y(n−1)+⋯+pn(x)y=f(x).
on any interval on which P0 has no zeros, with p1=P1/P0, …, pn=Pn/P0 and f=F/P0. For simplicity, throughout this chapter we’ll abbreviate the left side of Equation ??? by Ly; that is,
Ly=y(n)+p1(x)y(n−1)+⋯+pn(x)y.
When we speak of solutions of this equation and its homogeneous equation Ly=0, we mean solutions on (a,b). We’ll show how to use the method of variation of parameters to find a particular solution of Ly=f(x), provided that we know a fundamental set of solutions {y1,y2,…,yn} of Ly=0.
We seek a particular solution of Ly=f(x) in the form
yp=u1y1+u2y2+⋯+unyn
where {y1,y2,…,yn} is a known fundamental set of solutions of the homogeneous equation
y(n)+p1(x)y(n−1)+⋯+pn(x)y=0
and u1, u2, …, un are functions to be determined. We begin, as we did in Section 5.6, by imposing conditions that make our calculations easier. These conditions are similar to what we did in Section 5.6.
We impose the following n−1 conditions on u1,u2,…,un:
u′1y1+u′2y2+⋯+u′nyn=0u′1y′1+u′2y′2+⋯+u′ny′n=0u′1y(n1)+u′2y(n−1)2⋮⋯+u′ny(n−1)n=qu′1y(n−2)1+u′2y(n−2)2+⋯+u′ny(n−2)n=0.
These conditions lead to simple formulas for the first n−1 derivatives of yp:
y(r)p=u1y(r)1+u2y(r)2⋯+uny(r)n, 0≤r≤n−1.
These formulas are easy to remember, since they look as though we obtained them by differentiating Equation ??? n−1 times while treating u1, u2, …, un as constants. To see that Equation ??? implies Equation ???, we first differentiate Equation ??? to obtain
y′p=u1y′1+u2y′2+⋯+uny′n+u′1y1+u′2y2+⋯+u′nyn,
which reduces to
y′p=u1y′1+u2y′2+⋯+uny′n
because of the first equation in Equation ???. Differentiating this yields
y″p=u1y″1+u2y″2+⋯+uny″n+u′1y′1+u′2y′2+⋯+u′ny′n,
which reduces to
y″p=u1y″1+u2y″2+⋯+uny″n
because of the second equation in Equation ???. Continuing in this way yields Equation ???.
The last equation in Equation ??? is
y(n−1)p=u1y(n−1)1+u2y(n−1)2+⋯+uny(n−1)n.
Differentiating this yields
y(n)p=u1y(n)1+u2y(n)2+⋯+uny(n)n+u′1y(n−1)1+u′2y(n−1)2+⋯+u′ny(n−1)n.
Substituting this and Equation ??? into Equation ??? yields
u1Ly1+u2Ly2+⋯+unLyn+(u′1y(n−1)1+u′2y(n−1)2+⋯+u′ny(n−1)n)=f(x).
Since Lyi=0 (1≤i≤n), this reduces to
u′1y(n−1)1+u′2y(n−1)2+⋯+u′ny(n−1)n=f(x).
Combining this equation with Equation ??? shows that
yp=u1y1+u2y2+⋯+unyn
is a solution of Equation ??? if
u′1y1+u′2y2+⋯+u′nyn=0u′1y′1+u′2y′2+⋯+u′ny′n=0u′1y(n1)+u′2y(n−1)2⋮⋯+u′ny(n−1)n=qu′1y(n−2)1+u′2y(n−2)2+⋯+u′ny(n−2)n=0u′1y(n−1)1+u′2y(n−1)2+⋯+u′ny(n−1)n=f(x),
where W=[y1y2⋯yny′1y′2⋯y′n⋮⋮⋱⋮y(n−2)1y(n−2)2⋯y(n−2)ny(n−1)1y(n−1)2⋯y(n−1)n] and b=[00⋮0f(x)].
The determinant of W is the Wronskian of the fundamental set of solutions {y1,y2,…,yn}, which we know has no zeros on (a,b) since the solutions are linearly independent. Solving by Cramer’s rule yields
u′j=WjW,1≤j≤n,
where Wj is the determinant obtained by replacing the j-th column of W with b.
Having obtained u′1, u′2,…,u′n, we can integrate to obtain u1,u2,…,un. As in Section 5.6, we take the constants of integration to be zero, and we drop any linear combination of {y1,y2,…,yn} that may appear in yp.
Third Order Equations
If n=3, then
W=|y1y2y3y′1y′2y′3y″1y″2y″3|.
Therefore
W1=|0y2y30y′2y′3f(x)y″2y″3|,W2=|y10y3y′10y′3y″1f(x)y″3|,W3=|y1y20y′1y′20y″1y″2f(x)|,
and Equation ??? becomes
u′1=W1W,u′2=W2W,u′3=W3W.
Find the general solution of
xy‴−y″−xy′+y=8x2ex,
given that y1=x, y2=ex, and y3=e−x form a fundamental set of solutions of the homogeneous equation. Then find the general solution of Equation ???.
Solution
We seek a particular solution of Equation ??? of the form
yp=u1x+u2ex+u3e−x.
The Wronskian of {y1,y2,y3} is
W(x)=|xexe−x1ex−e−x0exe−x|=2x,
so
W1=|0exe−x0ex−e−x8xexexe−x|=−16xex,W2=|x0e−x10−e−x08xexe−x|=8x(x+1),W3=|xex01ex00ex8xex|=8xe2x(x−1).
Therefore, from Equation ???
u′1=W1W=−8ex,u′2=W2W=4(x+1),u′3=W3W=4e2x(x−1).
Integrating and taking the constants of integration to be zero yields
u1=−8ex,u2=2(x+1)2,u3=e2x(2x−3).
Hence,
yp=u1y1+u2y2+u3y3=(−8ex)x+ex(2(x+1)2)+e−x(e2x(2x−3))=ex(2x2−2x−1).
Since −ex is a solution of the complementary equation, we redefine
yp=2xex(x−1).
Therefore the general solution of Equation ??? is
y=2xex(x−1)+c1x+c2ex+c3e−x.
Fourth Order Equations
If n=4, then
W=|y1y2y3y4y′1y′2y′3y′4y″1y″2y″3y″4y‴1y‴2y‴3y‴4|,
Therefore
W1=|0y2y3y40y′2y′3y′40y″2y″3y″4f(x)y‴2y‴3y‴4|,W2=|y10y3y4y′10y′3y′4y″10y″3y″4y‴1f(x)y‴3y‴4|,
W3=|y1y20y4y′1y′20y′4y″1y″20y″4y‴1y‴2f(x)y‴4|,W4=|y1y2y30y′1y′2y′30y″1y″2y″30y‴1y‴2y‴3f(x)|,
and Equation ??? becomes
u′1=W1W,u′2=W2W,u′3=W3W,u′4=W4W.
Find the general solution of
x4y(4)+6x3y‴+2x2y″−4xy′+4y=12x2,
given that y1=x, y2=x2, y3=1/x and y4=1/x2 form a fundamental set of solutions of the homogeneous equation. Then find the general solution of Equation ??? on (−∞,0) and (0,∞).
Solution
We seek a particular solution of Equation ??? of the form
yp=u1x+u2x2+u3x+u4x2.
The Wronskian of {y1,y2,y3,y4} is
W(x)=|xx21/x−1/x212x−1/x2−2/x3022/x36/x400−6/x4−24/x5|=−72/x6,
so
W1=|0x21/x1/x202x−1/x2−2/x3022/x36/x412/x20−6/x4−24/x5,|=144/x6,W2=|x01/x1/x210−1/x2−2/x3002/x36/x4012/x2−6/x4−24/x5|=−72/x7,
W3=|xx20−1/x212x0−2/x30206/x40012/x2−24/x5|=−144/x4,W4=|xx21/x012x−1/x20022/x3000−6/x412/x2|=72/x3,
Therefore, from Equation ???,
u′1=W1W=−2,u′2=W2W=1/x,u′3=W3W=2x2,u′4=W4W=−x3.
Integrating these and taking the constants of integration to be zero yields
u1=−2x,u2=ln|x|,u3=2x33,u4=−x44.
Hence,
yp=u1y1+u2y2+u3y3+u4y4=(−2x)x+(ln|x|)x2+2x331x+(−x44)1x2=x2ln|x|−19x212.
Since −19x2/12 is a solution of the homogeneous equation, we redefine
yp=x2ln|x|.
Therefore
y=x2ln|x|+c1x+c2x2+c3x+c4x2
is the general solution of Equation ??? on (−∞,0) and (0,∞).