Differential Equations Guide
Differential Equations Guide
y ! x 3 " 2!x 3. Give an interval I of definition of E across the heart satisfies the linear differential equation
each of these solutions. Graph the solution curves. Is dE 1
there an initial-value problem whose solution is !# E.
dt RC
defined on (#$, $)?
(c) Is each IVP found in part (b) unique? That is, can Solve the DE subject to E(4) ! E 0.
there be more than one IVP for which, say,
y ! x 3 # 1!x 3, x in some interval I, is the solution?
44. In determining the integrating factor (5), we did not use Computer Lab Assignments
a constant of integration in the evaluation of "P(x) dx. 48. (a) Express the solution of the initial-value problem
Explain why using "P(x) dx " c has no effect on the y% # 2xy ! #1, y(0) ! 1( !2, in terms of erfc(x).
solution of (2).
(b) Use tables or a CAS to find the value of y(2). Use a
45. Suppose P(x) is continuous on some interval I and a is a CAS to graph the solution curve for the IVP on
number in I. What can be said about the solution of the (#$, $).
initial-value problem y% " P(x)y ! 0, y(a) ! 0?
49. (a) The sine integral function is defined by
Si(x) ! "x0 (sin t>t) dt, where the integrand is
Mathematical Models defined to be 1 at t ! 0. Express the solution y(x) of
the initial-value problem x 3 y% " 2x 2 y ! 10sin x,
46. Radioactive Decay Series The following system
y(1) ! 0 in terms of Si(x).
of differential equations is encountered in the study of the
decay of a special type of radioactive series of elements: (b) Use a CAS to graph the solution curve for the IVP
for x ' 0.
dx
! #&1x (c) Use a CAS to find the value of the absolute maxi-
dt mum of the solution y(x) for x ' 0.
dy
! &1x # &2 y, 50. (a) The Fresnel sine integral is defined by
dt x
S(x) ! " 0 sin(pt2>2) dt. Express the solution y(x)
where !1 and !2 are constants. Discuss how to solve this of the initial-value problem y%# (sin x 2 )y ! 0,
system subject to x(0) ! x 0 , y(0) ! y 0. Carry out your y(0) ! 5, in terms of S(x).
ideas. (b) Use a CAS to graph the solution curve for the IVP
on (#$, $).
47. Heart Pacemaker A heart pacemaker consists of a
switch, a battery of constant voltage E 0, a capacitor with (c) It is known that S(x) : 12 as x : $ and S(x) : #12
constant capacitance C, and the heart as a resistor with as x : #$ . What does the solution y(x) approach
constant resistance R. When the switch is closed, the as x : $? As x : #$?
capacitor charges; when the switch is open, the capacitor (d) Use a CAS to find the values of the absolute
discharges, sending an electrical stimulus to the heart. maximum and the absolute minimum of the
During the time the heart is being stimulated, the voltage solution y(x).
$f $f
dz ! dx # dy. (1)
$x $y
In the special case when f (x, y) ! c, where c is a constant, then (1) implies
$f $f
dx # dy ! 0. (2)
$x $y
Let M(x, y) and N(x, y) be continuous and have continuous first partial
derivatives in a rectangular region R defined by a % x % b, c % y % d. Then a
necessary and sufficient condition that M(x, y) dx # N(x, y) dy be an exact
differential is
$M $N
! . (4)
$y $x
64 ● CHAPTER 2 FIRST-ORDER DIFFERENTIAL EQUATIONS
PROOF OF THE NECESSITY For simplicity let us assume that M(x, y) and
N(x, y) have continuous first partial derivatives for all (x, y). Now if the expression
M(x, y) dx ! N(x, y) dy is exact, there exists some function f such that for all x in R,
"f "f
M(x, y) dx ! N(x, y) dy # dx ! dy.
"x "y
"f "f
Therefore M(x, y) # , N(x, y) # ,
"x "y
and
"M
"y
#
" "f
"y "x
#& '
"2 f
#
" "f
"y "x "x "y
"N
# .
"x & '
The equality of the mixed partials is a consequence of the continuity of the first par-
tial derivatives of M(x, y) and N(x, y).
The sufficiency part of Theorem 2.4.1 consists of showing that there exists a
function f for which "f !"x # M(x, y) and "f !"y # N(x, y) whenever (4) holds. The
construction of the function f actually reflects a basic procedure for solving exact
equations.
where the arbitrary function g(y) is the “constant” of integration. Now differentiate
(5) with respect to y and assume that "f !"y # N(x, y):
"f
#
"y "y
"
# M(x, y) dx ! g%(y) # N(x, y).
Finally, integrate (6) with respect to y and substitute the result in (5). The implicit
solution of the equation is f (x, y) # c.
Some observations are in order. First, it is important to realize that the expres-
sion N(x, y) $ ("!"y) " M(x, y) dx in (6) is independent of x, because
"
"x $
N(x, y) $
"
"y
# M(x, y) dx % # "N"x $ "y" &"x" # M(x, y) dx' # "N"x $ "M"y # 0.
Second, we could just as well start the foregoing procedure with the assumption that
"f !"y # N(x, y). After integrating N with respect to y and then differentiating that
result, we would find the analogues of (5) and (6) to be, respectively,
%M %N
# 2x # .
%y %x
Thus the equation is exact, and so by Theorem 2.4.1 there exists a function f (x, y)
such that
%f %f
# 2xy and # x2 " 1.
%x %y
f (x, y) # x2 y ! g(y).
Taking the partial derivative of the last expression with respect to y and setting the
result equal to N(x, y) gives
%f
# x2 ! g$(y) # x2 " 1. ; N(x, y)
%y
It follows that g$(y) # "1 and g(y) # "y. Hence f (x, y) # x 2 y " y, so the solution
of the equation in implicit form is x 2 y " y # c. The explicit form of the solution is
easily seen to be y # c!(1 " x 2) and is defined on any interval not containing either
x # 1 or x # "1.
%M %N
# 2e 2y ! xy sin xy " cos xy # .
%y %x
%f %f
M(x, y) # and N(x, y) # .
%x %y
Now for variety we shall start with the assumption that %f !%y # N(x, y); that is,
%f
# 2xe2y " x cos xy ! 2y
%y
Remember, the reason x can come out in front of the symbol ! is that in the integra-
tion with respect to y, x is treated as an ordinary constant. It follows that
%M %N
" #2xy " .
%y %x
%f
Now " y(1 # x2)
%y
y2
f(x, y) " (1 # x 2 ) $ h(x)
2
%f
" #xy2 $ h!(x) " cos x sin x # xy2.
%x
The last equation implies that h!(x) " cos x sin x. Integrating gives
# 1
h(x) " # (cos x)(#sin x dx) " # cos2 x.
2
y2 1
y Thus (1 # x2) # cos2 x " c1 or y2 (1 # x2) # cos2 x " c, (7)
2 2
where 2c1 has been replaced by c. The initial condition y " 2 when x " 0 demands
x
that 4(1) # cos 2 (0) " c, and so c " 3. An implicit solution of the problem is then
y 2 (1 # x 2 ) # cos 2 x " 3.
The solution curve of the IVP is the curve drawn in dark blue in Figure 2.4.1;
it is part of an interesting family of curves. The graphs of the members of the one-
parameter family of solutions given in (7) can be obtained in several ways, two of
FIGURE 2.4.1 Some graphs which are using software to graph level curves (as discussed in Section 2.2) and
of members of the family using a graphing utility to carefully graph the explicit functions obtained for var-
y 2(1 # x 2) # cos 2x " c ious values of c by solving y 2 " (c $ cos 2 x)"(1 # x 2) for y.
INTEGRATING FACTORS Recall from Section 2.3 that the left-hand side of
the linear equation y! $ P(x)y " f (x) can be transformed into a derivative when
we multiply the equation by an integrating factor. The same basic idea sometimes
works for a nonexact differential equation M(x, y) dx $ N(x, y) dy " 0. That is, it is
2.4 EXACT EQUATIONS ● 67
sometimes possible to find an integrating factor !(x, y) so that after multiplying, the
left-hand side of
!(x, y)M(x, y) dx ! !(x, y)N(x, y) dy " 0 (8)
is an exact differential. In an attempt to find !, we turn to the criterion (4) for exact-
ness. Equation (8) is exact if and only if (!M)y " (!N)x , where the subscripts
denote partial derivatives. By the Product Rule of differentiation the last equation is
the same as ! My ! ! y M " !Nx ! ! x N or
! x N # ! y M " (My # Nx)!. (9)
Although M, N, My , and Nx are known functions of x and y, the difficulty here in
determining the unknown !(x, y) from (9) is that we must solve a partial differential
equation. Since we are not prepared to do that, we make a simplifying assumption.
Suppose ! is a function of one variable; for example, say that ! depends only on x. In
this case, ! x " d!!dx and ! y " 0, so (9) can be written as
d$ My # Nx
" $. (10)
dx N
We are still at an impasse if the quotient (My # Nx )!N depends on both x and y.
However, if after all obvious algebraic simplifications are made, the quotient
(My # Nx )!N turns out to depend solely on the variable x, then (10) is a first-order
ordinary differential equation. We can finally determine ! because (10) is separa-
ble as well as linear. It follows from either Section 2.2 or Section 2.3 that
!(x) " e "((My#Nx)/N )dx. In like manner, it follows from (9) that if ! depends only on
the variable y, then
d$ Nx # My
" $. (11)
dy M
In this case, if (N x # My)!M is a function of y only, then we can solve (11) for !.
We summarize the results for the differential equation
M(x, y) dx ! N(x, y) dy " 0. (12)
• If (My # Nx)!N is a function of x alone, then an integrating factor for (12) is
$(x) " e
#
My#Nx
N
dx
. (13)
• If (Nx # My)!M is a function of y alone, then an integrating factor for (12) is
$(y) " e
# Nx#My
M
dy
. (14)
The integrating factor is then e !3dy/y ! e 3lny ! e lny ! y 3. After we multiply the given
3
REMARKS
(i) When testing an equation for exactness, make sure it is of the precise
form M(x, y) dx " N(x, y) dy ! 0. Sometimes a differential equation
is written G(x, y) dx ! H(x, y) dy. In this case, first rewrite it as
G(x, y) dx # H(x, y) dy ! 0 and then identify M(x, y) ! G(x, y) and
N(x, y) ! #H(x, y) before using (4).
(ii) In some texts on differential equations the study of exact equations
precedes that of linear DEs. Then the method for finding integrating factors
just discussed can be used to derive an integrating factor for
y$ " P(x)y ! f (x). By rewriting the last equation in the differential form
(P(x)y # f (x)) dx " dy ! 0, we see that
M y # Nx
! P(x).
N
From (13) we arrive at the already familiar integrating factor e !P(x)dx, used in
Section 2.3.
6. " 1
2y # " cos 3x
x
dy y
#
" 2 # 4x3 " 3y sin 3x ! 0
dx x 16. (5y # 2x)y$ # 2y ! 0
2 2 2
7. (x # y ) dx " (x # 2xy) dy ! 0 17. (tan x # sin x sin y) dx " cos x cos y dy ! 0
2
! (x # sin2 x # 4xye xy ) dy
9. (x # y 3 " y 2 sin x) dx ! (3xy 2 " 2y cos x) dy
10. (x 3 " y 3) dx " 3xy 2 dy ! 0 19. (4t 3 y # 15t 2 # y) dt " (t 4 " 3y 2 # t) dy ! 0
! "
2 dy#dx # g(x)h(y) is exact.
34. cos x dx ! 1 ! sin x dy # 0
y
dy %g dx %g du dy du
" ! gives " gx (x, u) ! gu(x, u) .
dx %x dx %u dx dx dx
If we replace dy#dx by the foregoing derivative and replace y in f (x, y) by g (x, u), then
du
the DE dy#dx " f (x, y) becomes g x (x, u) ! gu (x, u) " f (x, g(x, u)), which, solved
dx
du
for du#dx, has the form " F(x, u). If we can determine a solution u " !(x) of this
dx
last equation, then a solution of the original differential equation is y " g(x, !(x)).
In the discussion that follows we examine three different kinds of first-order
differential equations that are solvable by means of a substitution.