MATH2352 Differential Equations and Applications Tutorial Notes 3
MATH2352 Differential Equations and Applications Tutorial Notes 3
MATH2352 Differential Equations and Applications Tutorial Notes 3
Tutorial Notes 3
Samsung, SUM Sung Fung
6 July, 2012
10 Reduction of Order
For Linear Homogeneous Equation ONLY!!
Suppose we know y1 (t) is one of the two linear independent solutions of the equation
y (t) + p(t)y (t) + q(t)y(t) = 0
where p(t) and q(t) are some functions in t, we can use the Reduction of Order Method to find out
another linear independent solution.
Procedures:
1. Consider the test solution y = v(t)y1 (t), then
y (t) = v (t)y1 (t) + v(t)y1 (t)
y (t) = v (t)y1 (t) + 2v (t)y1 (t) + v(t)y1 (t)
and substitute these into the differential equation to obtain
y1 v + (2y1 + py1 )v + (y1 + py1 + qy1 )v = 0
2. Since y1 is a solution to the given DE, we have y1 + py1 + qy1 = 0, so the equation of v(t) can be
reduced into
y1 v + (2y1 + py1 )v = 0
which is a first order equation for the function v .
3. Perform integration to obtain v .
4. Integrating v once more to obtain v(t), and hence getting y2 (t) = v(t)y1 (t).
5. The general solution becomes y(t) = A1 y1 (t) + A2 y2 (t).
Example:
Given y1 (t) = et , to find the general solution of the equation
(t 1)y ty + y = 0 (t > 1)
t
Let y = v(t)e and substitute this into the DE, upon simplifications we have
(t 1)v + (t 2)v = 0
Solving this equation in v yields
v = A(t 1)et
Integrate once more, we get
v(t) = Atet + B y(t) = At + Bet
Note that we can set A = 1 and B = 0 because Bet and y1 are LD. So, y2 (t) = t and the
general solution becomes
y(t) = A1 et + A2 t
1
11 Inhomogeneous Linear Equations
Given
N
X
L[y(t)] = n (t)y (n) (t) = g(t)
n=0
is a linear inhomogeneous differential equation, then the general solution of L[y] = g can be written in
the form
y(t) = Yh (t) + Yp (t)
where Yh is the general homogeneous solution satisfying L[Yh (t)] = 0 and Yp is a particular solution that
satisfies L[Yp (t)] = g(t).
Thus, in order to solve L[y] = g, we have to solve the homogeneous equation L[Yh ] = 0 and use a few
tricks to create a Yp .
g(t) Yp (t)
Pn (t) = A0 + A1 t + + An tn ts (B0 + B1 t + + Bn tn )
Here s is the smallest nonnegative integer (s = 0, 1, 2) that Yp does NOT contain any homogeneous
solution.
Example:
Find a particular solution to the equation y + y = 2tet .
The general homogeneous solution to the given DE is
Consider a trial solution of the form Y = (B0 + B1 t)et , plugging Y into the given DE gives
2
Example:
Solve the following IVP
y (x) 2y (x) + y(x) = 24x2 ex
y(0) = 3
y (0) = 0
Yh (x) = A1 ex + A2 xex
Y = (B1 + B2 x + B3 x2 )ex
Y 2Y + Y = 2B3 ex 6= 24x2 ex
The guess Y is NOT a good guess because the first two terms constitute the homogeneous
solution.
To obtain a valid guess, we have to multiply xn (with some n N) to the previous guess so
that there is no homogeneous part left, in this case, n = 2 should be used:
Y = (B1 x2 + B2 x3 + B3 x4 )ex
Substituting the general solution into the ICs, we get A1 = 3 and A2 = 3, so the solution
of the IVP is
y(x) = 3ex 3xex + 2x4 ex
3
13 Method of Variation of Parameters
Another useful trick for finding Yp for y (t) + p(t)y (t) + q(t)y(t) = g(t) is to guess the solution using
the Method of Variation of Parameters.
Procedures:
1. Consider the test solution y = u1 (t)y1 (t) + u2 (t)y2 (t), then
y (t) = u1 (t)y1 (t) + u2 (t)y2 (t) +u1 (t)y1 (t) + u2 (t)y2 (t)
| {z }
set to 0
We can set the sum of the first two terms to zero to simplify the procedures later on. Also, we
have
y (t) = u1 (t)y1 (t) + u2 (t)y2 (t) + u1 (t)y1 (t) + u2 (t)y2 (t)
Substitute these into the differential equation to obtain
u1 y1 + u2 y2 + u1 (y1 + py1 + qy1 ) + u2 (y2 + py2 + qy2 ) = g(t)
2. Since y1 and y2 are homogeneous solutions to the given DE, we have yi + pyi + qyi = 0 for i = 1, 2,
so the above equation can be reduced into
u1 y1 + u2 y2 = g
Combining with the equation u1 y1 + u2 y2 = 0, we have the system of equations
y1 y2 u1 0
=
y1 y2 u2 g
Example:
To find a particular solution to the equation y + y = sec t for t 2 , 2 , of which the
general homogeneous solution is Yh (t) = A1 cos t + A2 sin t, denote y1 = cos t and y2 = sin t.
Suppose Y = u1 (t)y1 (t) + u2 (t)y2 (t) is a particular solution in the method of variation of
parameters, then the corresponding matrix equation for u1 and u2 is
cos t sin t u1 0
=
sin t cos t u2 sec t
which upon solving gives u1 = tan t and u2 = 1.
Integrating u1 and u2 once more gives
u1 = ln(cos t), u2 = t
where the constants were dropped as they will eventually produce homogeneous solutions.
So a particular solution to the given equation is
Yp (t) = ln(cos t) cos t + t sin t