SlideShare a Scribd company logo
Chapter 1
Introduction to Differential Equations
1.1 Definitions and Terminology
1. Second order; linear
2. Third order; nonlinear because of (dy/dx)4
3. Fourth order; linear
4. Second order; nonlinear because of cos(r + u)
5. Second order; nonlinear because of (dy/dx)2 or

1 + (dy/dx)2
6. Second order; nonlinear because of R2
7. Third order; linear
8. Second order; nonlinear because of ẋ2
9. First order; nonlinear because of sin

dy
dx

10. First order; linear in x.
11. Writing the differential equation in the form x(dy/dx) + y2 = 1, we see that it is nonlinear
in y because of y2. However, writing it in the form (y2 − 1)(dx/dy) + x = 0, we see that it is
linear in x.
12. Writing the differential equation in the form u(dv/du) + (1 + u)v = ueu, we see that it is
linear in v. However, writing it in the form (v + uv − ueu)(du/dv) + u = 0, we see that it is
nonlinear in u.
13. From y = e−x/2 we obtain y = −1
2e−x/2. Then 2y + y = −e−x/2 + e−x/2 = 0.
14. From y = 6
5 − 6
5 e−20t we obtain dy/dt = 24e−20t, so that
dy
dt
+ 20y = 24e−20t
+ 20

6
5
−
6
5
e−20t

= 24.
1
Contact me in order to access the whole complete document.
WhatsApp: https://wa.me/message/2H3BV2L5TTSUF1
Email: smtb98@gmail.com
Telegram: https://t.me/solutionmanual
s
m
t
b
9
8
@
g
m
a
i
l
.
c
o
m
s
m
t
b
9
8
@
g
m
a
i
l
.
c
o
m
complete document is available on https://unihelp.xyz/ *** contact me if site not loaded
2 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS
15. From y = e3x cos 2x we obtain y = 3e3x cos 2x−2e3x sin 2x and y = 5e3x cos 2x−12e3x sin 2x,
so that y − 6y + 13y = 0.
16. From y = − cos x ln(sec x + tan x) we obtain y
= −1 + sin x ln(sec x + tan x) and
y
= tan x + cos x ln(sec x + tan x). Then y + y = tan x.
17. The domain of the function, found by solving x+2 ≥ 0, is [−2, ∞). From y = 1+2(x+2)−1/2
we have
(y − x)y
= (y − x)[1 + (2(x + 2)−1/2
]
= y − x + 2(y − x)(x + 2)−1/2
= y − x + 2[x + 4(x + 2)1/2
− x](x + 2)−1/2
= y − x + 8(x + 2)1/2
(x + 2)−1/2
= y − x + 8.
An interval of definition for the solution of the differential equation is (−2, ∞) because y is
not defined at x = −2.
18. Since tan x is not defined for x = π/2 + nπ, n an integer, the domain of y = 5 tan 5x is
{x

 5x = π/2 + nπ} or {x

 x = π/10 + nπ/5}. From y = 25 sec2 5x we have
y
= 25(1 + tan2
5x) = 25 + 25 tan2
5x = 25 + y2
.
An interval of definition for the solution of the differential equation is (−π/10, π/10). Another
interval is (π/10, 3π/10), and so on.
19. The domain of the function is {x

 4 − x2 = 0} or {x

 x = −2 or x = 2}. From y =
2x/(4 − x2)2 we have
y
= 2x

1
4 − x2
2
= 2xy2
.
An interval of definition for the solution of the differential equation is (−2, 2). Other intervals
are (−∞, −2) and (2, ∞).
20. The function is y = 1/
√
1 − sin x , whose domain is obtained from 1 − sin x = 0 or sin x = 1.
Thus, the domain is {x

 x = π/2 + 2nπ}. From y = −1
2(1 − sin x)−3/2(− cos x) we have
2y
= (1 − sin x)−3/2
cos x = [(1 − sin x)−1/2
]3
cos x = y3
cos x.
An interval of definition for the solution of the differential equation is (π/2, 5π/2). Another
one is (5π/2, 9π/2), and so on.
1.1 Definitions and Terminology 3
21. Writing ln(2X − 1) − ln(X − 1) = t and differentiating
implicitly we obtain
2
2X − 1
dX
dt
−
1
X − 1
dX
dt
= 1

2
2X − 1
−
1
X − 1

dX
dt
= 1
2X − 2 − 2X + 1
(2X − 1) (X − 1)
dX
dt
= 1
dX
dt
= −(2X − 1)(X − 1) = (X − 1)(1 − 2X).
– 4 –2 2 4
t
– 4
–2
2
4
x
Exponentiating both sides of the implicit solution we obtain
2X − 1
X − 1
= et
2X − 1 = Xet
− et
(et
− 1) = (et
− 2)X
X =
et − 1
et − 2
.
Solving et − 2 = 0 we get t = ln 2. Thus the solution is defined on (−∞, ln 2) or on (ln 2, ∞).
The graph of the solution defined on (−∞, ln 2) is dashed, and the graph of the solution
defined on (ln 2, ∞) is solid.
22. Implicitly differentiating the solution, we obtain
−2x2 dy
dx
− 4xy + 2y
dy
dx
= 0
−x2
dy − 2xy dx + y dy = 0
2xy dx + (x2
− y)dy = 0.
Using the quadratic formula to solve y2 − 2x2y − 1 = 0
for y, we get y =

2x2 ±
√
4x4 + 4

/2 = x2 ±
√
x4 + 1 .
Thus, two explicit solutions are y1 = x2 +
√
x4 + 1 and
y2 = x2 −
√
x4 + 1 . Both solutions are defined on (−∞, ∞).
The graph of y1(x) is solid and the graph of y2 is dashed.
– 4 –2 2 4
x
– 4
–2
2
4
y
4 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS
23. Differentiating P = c1et/

1 + c1et

we obtain
dP
dt
=

1 + c1et

c1et − c1et · c1et
(1 + c1et)2 =
c1et
1 + c1et

1 + c1et

− c1et
1 + c1et
=
c1et
1 + c1et
1 −
c1et
1 + c1et
= P(1 − P).
24. Differentiating y = 2x2 − 1 + c1e−2x2
we obtain
dy
dx
= 4x − 4xc1e−2x2
, so that
dy
dx
+ 4xy = 4x − 4xc1e−2x2
+ 8x3
− 4x + 4c1xe−x2
= 8x3
25. From y = c1e2x + c2xe2x we obtain
dy
dx
= (2c1 + c2)e2x
+ 2c2xe2x
and
d2y
dx2
= (4c1 + 4c2)e2x
+
4c2xe2x
, so that
d2y
dx2
− 4
dy
dx
+ 4y = (4c1 + 4c2 − 8c1 − 4c2 + 4c1)e2x
+ (4c2 − 8c2 + 4c2)xe2x
= 0.
26. From y = c1x−1 + c2x + c3x ln x + 4x2 we obtain
dy
dx
= −c1x−2
+ c2 + c3 + c3 ln x + 8x,
d2y
dx2
= 2c1x−3
+ c3x−1
+ 8,
and
d3y
dx3
= −6c1x−4
− c3x−2
,
so that
x3 d3y
dx3
+ 2x2 d2y
dx2
− x
dy
dx
+ y = (−6c1 + 4c1 + c1 + c1)x−1
+ (−c3 + 2c3 − c2 − c3 + c2)x
+ (−c3 + c3)x ln x + (16 − 8 + 4)x2
= 12x2
In Problems 27–30, we use the Product Rule and the derivative of an integral ((12) of this section):
d
dx
x
a
g(t) dt = g(x).
27. Differentiating y = e3x
x
1
e−3t
t
dt we obtain
dy
dx
= e3x
x
1
e−3t
t
dt +
e−3t
x
· e3x
or
dy
dx
= e3x
x
1
e−3t
t
dt +
1
x
, so that
x
dy
dx
− 3xy = x

e3x
x
1
e−3t
t
dt +
1
x

− 3x

e3x
x
1
e−3t
t
dt

= xe3x
x
1
e−3t
t
dt + 1 − 3xe3x
x
1
e−3t
t
dt = 1
1.1 Definitions and Terminology 5
28. Differentiating y =
√
x
x
4
cos t
√
t
dt we obtain
dy
dx
=
1
2
√
x
x
4
cos t
√
t
dt +
cos x
√
x
·
√
x or
dy
dx
=
1
2
√
x
x
4
cos t
√
t
dt + cos x, so that
2x
dy
dx
− y = 2x

1
2
√
x
x
4
cos t
√
t
dt + cos x

−
√
x
x
4
cos t
√
t
dt
=
√
x
x
4
cos t
√
t
dt + 2x cos x −
√
x
x
4
cos t
√
t
dt = 2x cos x
29. Differentiating y =
5
x
+
10
x
x
1
sin t
t
dt we obtain
dy
dx
= −
5
x2
−
10
x2
x
1
sin t
t
dt +
sin x
x
·
10
x
or
dy
dx
= −
5
x2
−
10
x2
x
1
sin t
t
dt +
10 sin x
x2
, so that
x2 dy
dx
+ xy = x2

−
5
x2
−
10
x2
x
1
sin t
t
dt +
10 sin x
x2

+ x

5
x
+
10
x
x
1
sin t
t
dt

= −5 − 10
x
1
sin t
t
dt + 10 sin x + 5 + 10
x
1
sin t
t
dt = 10 sin x
30. Differentiating y = e−x2
+e−x2
x
0
et2
dt we obtain
dy
dx
= −2xe−x2
−2xe−x2
x
0
et2
dt+ex2
·e−x2
or
dy
dx
= −2xe−x2
− 2xe−x2
x
0
et2
dt + 1, so that
dy
dx
+ 2xy =

−2xe−x2
− 2xe−x2
x
0
et2
dt + 1

+ 2x

e−x2
+ e−x2
x
0
et2
dt

= −2xe−x2
− 2xe−x2
x
0
et2
dt + 1 + 2xe−x2
+ 2xe−x2
x
0
et2
dt = 1
31. From
y =
−x2, x  0
x2, x ≥ 0
we obtain
y
=
−2x, x  0
2x, x ≥ 0
so that xy − 2y = 0.
32. The function y(x) is not continuous at x = 0 since lim
x→0−
y(x) = 5 and lim
x→0+
y(x) = −5. Thus,
y(x) does not exist at x = 0.
33. Force the function y = emx into the equation y + 2y = 0 to get
(emx
)
+ 2(emx
) = 0
memx
+ 2emx
= 0
emx
(m + 2) = 0
Now since emx  0 for all values of x, we must have m = −2 and so y = e−2x is a solution.
6 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS
34. Force the function y = emx into the equation 3y − 4y = 0 to get
3(emx
)
− 4(emx
) = 0
3memx
− 4emx
= 0
emx
(3m − 4) = 0
Now since emx  0 for all values of x, we must have m = 4/3 and so y = e4x/3 is a solution.
35. Force the function y = emx into the equation y − 5y + 6y = 0 to get
(emx
)
− 5(emx
)
+ 6(emx
) = 0
m2
emx
− 5memx
+ 6emx
= 0
emx
(m2
− 5m + 6) = 0
emx
(m − 2)(m − 3) = 0
Now since emx  0 for all values of x, we must have m = 2 and m = 3 therefore y = e2x and
y = e3x are solutions.
36. Force the function y = emx into the equation 2y + 9y − 5y = 0 to get
2(emx
)
+ 9(emx
)
− 5(emx
) = 0
2m2
emx
+ 9memx
− 5emx
= 0
emx
(2m2
+ 9m − 5) = 0
emx
(m + 5)(2m − 1) = 0
Now since emx  0 for all values of x , we must have m = −5 and m = 1/2 therefore y = e−5x
and y = ex/2 are solutions.
37. Force the function y = xm into the equation xy + 2y = 0 to get
x · (xm
)
+ 2(xm
)
= 0
x · m(m − 1)xm−2
+ 2mxm−1
= 0
(m2
− m)xm−1
+ 2mxm−1
= 0
xm−1
[m2
+ m] = 0
xm−1
[m(m + 1)] = 0
The last line implies that m = 0 and m = −1 therefore y = x0 = 1 and y = x−1 are solutions.
1.1 Definitions and Terminology 7
38. Force the function y = xm into the equation 4x2y + y = 0 to get
4x2
· (xm
)
+ (xm
) = 0
4x2
· m(m − 1)xm−2
+ xm
= 0
4(m2
− m)xm
+ xm
= 0
xm
[4m2
− 4m + 1] = 0
xm
[(2m − 1)(2m − 1)] = 0
The last line implies m = 1/2 therefore y = x1/2 is a solution.
39. Force the function y = xm into the equation x2y − 7xy + 15y = 0 to get
x2
· (xm
)
− 7x · (xm
)
+ 15(xm
) = 0
x2
· m(m − 1)xm−2
− 7x · mxm−1
+ 15xm
= 0
(m2
− m)xm
− 7mxm
+ 15xm
= 0
xm
[m2
− 8m + 15] = 0
xm
[(m − 3)(m − 5)] = 0
The last line implies that m = 3 and m = 5 therefore y = x3 and y = x5 are solutions.
40. Force the function y = xm into the equation x2y − 3xy + 3y = 0 to get
x2
· (xm
)
− 3x · (xm
)
+ 3(xm
)
= 0
x2
· m(m − 1)(m − 2)xm−3
− 3x · m(m − 1)xm−2
+ 3 · mxm−1
= 0
(m3
− 3m2
+ 2m)xm−1
− 3(m2
− m)xm−1
+ 3mxm−1
= 0
xm−1
[m3
− 6m2
+ 8m] = 0
xm
[m(m − 2)(m − 4)] = 0
The last line implies that m = 0, m = 2 and m = 4 therefore y = x0 = 1, y = x2, and y = x4
are solutions.
In Problems 41–44, we substitute y = c into the differential equations and use y = 0 and y = 0.
41. Solving 5c = 10 we see that y = 2 is a constant solution.
42. Solving c2 + 2c − 3 = (c + 3)(c − 1) = 0 we see that y = −3 and y = 1 are constant solutions.
43. Since 1/(c − 1) = 0 has no solutions, the differential equation has no constant solutions.
8 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS
44. Solving 6c = 10 we see that y = 5/3 is a constant solution.
45. From y = (x + c1)2 we obtain

dy
dx
2
= (2 (x + c1))2
= 4 (x + c1)2
.
Then
4y = 4(x + c1)2
.
Inspection of the differential equation reveals that y = 0 is a solution of the differential
equation but is not a member of the one-parameter family y = (x + c1)2.
46. From y = 3 sin (x + c1) we obtain

dy
dx
2
= (3 cos (x + c1))2
= 9 cos2
(x + c1).
Then
9 − y2
= 9 − (3 sin (x + c1))2
= 9 − 9 sin2
(x + c1) = 9

1 − sin2
(x + c1)

= 9 cos2
(x + c1).
Inspection of the differential equation reveals that y = 3 and y = −3 are solutions of the
differential equation but not members of the one-parameter family y = 3 sin (x + c1).
47. From x −

16 − y2 = c1 we obtain
d
dx

x −

16 − y2

=
d
dx
(c1)
1 −
1
2

16 − y2
−1/2

−2y
dy
dx

= 0
1 +
y

16 − y2
dy
dx
= 0
y

16 − y2
dy
dx
= −1
y
dy
dx
= −

16 − y2.
Then
y
dy
dx
+

16 − y2 = −

16 − y2 +

16 − y2 = 0.
Inspection of the differential equation reveals that y = −4 and y = 4 are solutions of the given
differential equation but are not members of the one-parameter family x −

16 − y2 = c1.
1.1 Definitions and Terminology 9
48. From y = x − (x − c1)2 we obtain

dy
dx
2
= (1 − 2 (x − c1))2
= 1 − 4 (x − c1) + 4 (x − c1)2
= 1 − 4x + 4c1 + 4x2
− 8xc1 + 4c2
1
−2
dy
dx
= −2 (1 − 2 (x − c1)) = −2 + 4x − 4c1
4y = 4

x − (x − c1)2

= 4

x − x2
+ 2xc1 − c2
1

= 4x − 4x2
+ 8xc1 − 4c2
1.
Then

dy
dx
2
− 2
dy
dx
+ 4y = 4x − 1.
Inspection of the differential equation reveals that y = x is a solution of the differential
equation but not a member of the one-parameter family y = x − (x − c1)2.
49. From x = e−2t + 3e6t and y = −e−2t + 5e6t we obtain
dx
dt
= −2e−2t
+ 18e6t
and
dy
dt
= 2e−2t
+ 30e6t
.
Then
x + 3y = (e−2t
+ 3e6t
) + 3(−e−2t
+ 5e6t
) = −2e−2t
+ 18e6t
=
dx
dt
and
5x + 3y = 5(e−2t
+ 3e6t
) + 3(−e−2t
+ 5e6t
) = 2e−2t
+ 30e6t
=
dy
dt
.
50. From x = cos 2t + sin 2t + 1
5et and y = − cos 2t − sin 2t − 1
5 et we obtain
dx
dt
= −2 sin 2t + 2 cos 2t +
1
5
et
and
dy
dt
= 2 sin 2t − 2 cos 2t −
1
5
et
and
d2x
dt2
= −4 cos 2t − 4 sin 2t +
1
5
et
and
d2y
dt2
= 4 cos 2t + 4 sin 2t −
1
5
et
.
Then
4y + et
= 4(− cos 2t − sin 2t −
1
5
et
) + et
= −4 cos 2t − 4 sin 2t +
1
5
et
=
d2x
dt2
and
4x − et
= 4(cos 2t + sin 2t +
1
5
et
) − et
= 4 cos 2t + 4 sin 2t −
1
5
et
=
d2y
dt2
.
51. (y)2 + 1 = 0 has no real solutions because (y)2 + 1 is positive for all differentiable functions
y = φ(x).
52. The only solution of (y)2 + y2 = 0 is y = 0, since if y = 0, y2  0 and (y)2 + y2 ≥ y2  0.
10 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS
53. The first derivative of f(x) = ex is ex. The first derivative of f(x) = ekx is kekx. The
differential equations are y = y and y = ky, respectively.
54. Any function of the form y = cex or y = ce−x is its own second derivative. The corresponding
differential equation is y −y = 0. Functions of the form y = c sin x or y = c cos x have second
derivatives that are the negatives of themselves. The differential equation is y + y = 0.
55. We first note that

1 − y2 =

1 − sin2
x =
√
cos2 x = | cos x|. This prompts us to consider
values of x for which cos x  0, such as x = π. In this case
dy
dx





x=π
=
d
dx
(sin x)





x=π
= cos x


x=π
= cos π = −1,
but

1 − y2|x=π =

1 − sin2
π =
√
1 = 1.
Thus, y = sin x will only be a solution of y =

1 − y2 when cos x  0. An interval of
definition is then (−π/2, π/2). Other intervals are (3π/2, 5π/2), (7π/2, 9π/2), and so on.
56. Since the first and second derivatives of sin t and cos t involve sin t and cos t, it is plausible that
a linear combination of these functions, A sin t+B cos t, could be a solution of the differential
equation. Using y = A cos t − B sin t and y = −A sin t − B cos t and substituting into the
differential equation we get
y
+ 2y
+ 4y = −A sin t − B cos t + 2A cos t − 2B sin t + 4A sin t + 4B cos t
= (3A − 2B) sin t + (2A + 3B) cos t = 5 sin t
Thus 3A − 2B = 5 and 2A + 3B = 0. Solving these simultaneous equations we find A = 15
13
and B = −10
13 . A particular solution is y = 15
13 sin t − 10
13 cos t.
57. One solution is given by the upper portion of the graph with domain approximately (0, 2.6).
The other solution is given by the lower portion of the graph, also with domain approximately
(0, 2.6).
58. One solution, with domain approximately (−∞, 1.6) is the portion of the graph in the second
quadrant together with the lower part of the graph in the first quadrant. A second solution,
with domain approximately (0, 1.6) is the upper part of the graph in the first quadrant. The
third solution, with domain (0, ∞), is the part of the graph in the fourth quadrant.
1.1 Definitions and Terminology 11
59. Differentiating (x3 + y3)/xy = 3c we obtain
xy(3x2 + 3y2y) − (x3 + y3)(xy + y)
x2y2
= 0
3x3
y + 3xy3
y
− x4
y
− x3
y − xy3
y
− y4
= 0
(3xy3
− x4
− xy3
)y
= −3x3
y + x3
y + y4
y
=
y4 − 2x3y
2xy3 − x4
=
y(y3 − 2x3)
x(2y3 − x3)
.
60. A tangent line will be vertical where y is undefined, or in this case, where x(2y3 − x3) = 0.
This gives x = 0 and 2y3 = x3. Substituting y3 = x3/2 into x3 + y3 = 3xy we get
x3
+
1
2
x3
= 3x

1
21/3
x

3
2
x3
=
3
21/3
x2
x3
= 22/3
x2
x2
(x − 22/3
) = 0.
Thus, there are vertical tangent lines at x = 0 and x = 22/3, or at (0, 0) and (22/3, 21/3).
Since 22/3 ≈ 1.59, the estimates of the domains in Problem 58 were close.
61. The derivatives of the functions are φ
1(x) = −x/
√
25 − x2 and φ
2(x) = x/
√
25 − x2, neither
of which is defined at x = ±5.
62. To determine if a solution curve passes through (0, 3) we let t = 0 and P = 3 in the equation
P = c1et/(1 + c1et). This gives 3 = c1/(1 + c1) or c1 = −3
2 . Thus, the solution curve
P =
(−3/2)et
1 − (3/2)et
=
−3et
2 − 3et
passes through the point (0, 3). Similarly, letting t = 0 and P = 1 in the equation for the
one-parameter family of solutions gives 1 = c1/(1 + c1) or c1 = 1 + c1. Since this equation
has no solution, no solution curve passes through (0, 1).
63. For the first-order differential equation integrate f(x). For the second-order differential equa-
tion integrate twice. In the latter case we get y =

(

f(x)dx) dx + c1x + c2.
64. Solving for y using the quadratic formula we obtain the two differential equations
y
=
1
x

2 + 2

1 + 3x6

and y
=
1
x

2 − 2

1 + 3x6

,
so the differential equation cannot be put in the form dy/dx = f(x, y).
12 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS
65. The differential equation yy − xy = 0 has normal form dy/dx = x. These are not equivalent
because y = 0 is a solution of the first differential equation but not a solution of the second.
66. Differentiating we get y = c1 + 3c2x2 and y = 6c2x. Then c2 = y/6x and c1 = y − xy/2,
so
y =

y
−
xy
2

x +

y
6x

x3
= xy
−
1
3
x2
y
and the differential equation is x2y − 3xy + 3y = 0.
67. (a) Since e−x2
is positive for all values of x, dy/dx  0 for all x, and a solution, y(x), of the
differential equation must be increasing on any interval.
(b) lim
x→−∞
dy
dx
= lim
x→−∞
e−x2
= 0 and lim
x→∞
dy
dx
= lim
x→∞
e−x2
= 0. Since dy/dx approaches 0 as
x approaches −∞ and ∞, the solution curve has horizontal asymptotes to the left and
to the right.
(c) To test concavity we consider the second derivative
d2y
dx2
=
d
dx

dy
dx

=
d
dx

e−x2

= −2xe−x2
.
Since the second derivative is positive for x  0 and negative for x  0, the solution
curve is concave up on (−∞, 0) and concave down on (0, ∞).
(d)
x
y
68. (a) The derivative of a constant solution y = c is 0, so solving 5 − c = 0 we see that c = 5
and so y = 5 is a constant solution.
(b) A solution is increasing where dy/dx = 5 − y  0 or y  5. A solution is decreasing
where dy/dx = 5 − y  0 or y  5.
69. (a) The derivative of a constant solution is 0, so solving y(a − by) = 0 we see that y = 0 and
y = a/b are constant solutions.
(b) A solution is increasing where dy/dx = y(a − by) = by(a/b − y)  0 or 0  y  a/b. A
solution is decreasing where dy/dx = by(a/b − y)  0 or y  0 or y  a/b.
1.1 Definitions and Terminology 13
(c) Using implicit differentiation we compute
d2y
dx2
= y(−by
) + y
(a − by) = y
(a − 2by).
Solving d2y/dx2 = 0 we obtain y = a/2b. Since d2y/dx2  0 for 0  y  a/2b and
d2y/dx2  0 for a/2b  y  a/b, the graph of y = φ(x) has a point of inflection at
y = a/2b.
(d)
y = a/b
y = 0
x
y
70. (a) If y = c is a constant solution then y = 0, but c2 + 4 is never 0 for any real value of c.
(b) Since y = y2 + 4  0 for all x where a solution y = φ(x) is defined, any solution must
be increasing on any interval on which it is defined. Thus it cannot have any relative
extrema.
(c) Using implicit differentiation we compute d2y/dx2 = 2yy = 2y(y2 + 4). Setting
d2y/dx2 = 0 we see that y = 0 corresponds to the only possible point of inflection.
Since d2y/dx2  0 for y  0 and d2y/dx2  0 for y  0, there is a point of inflection
where y = 0.
(d)
x
y
14 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS
71. In Mathematica use
Clear[y]
y[x ]:= x Exp[5x] Cos[2x]
y[x]
y''''[x] − 20y'''[x] + 158y''[x] − 580y'[x] +841y[x]//Simplify
The output will show y(x) = e5xx cos 2x, which verifies that the correct function was entered,
and 0, which verifies that this function is a solution of the differential equation.
72. In Mathematica use
Clear[y]
y[x ]:= 20Cos[5Log[x]]/x − 3Sin[5Log[x]]/x
y[x]
xˆ3 y'''[x] + 2xˆ2 y''[x] + 20x y'[x] − 78y[x]//Simplify
The output will show y(x) = 20 cos(5 ln x)/x − 3 sin(5 ln x)/x, which verifies that the correct
function was entered, and 0, which verifies that this function is a solution of the differential
equation.
1.2 Initial-Value Problems
1. Solving −1/3 = 1/(1 + c1) we get c1 = −4. The solution is y = 1/(1 − 4e−x
).
2. Solving 2 = 1/(1 + c1e) we get c1 = −(1/2)e−1. The solution is y = 2/(2 − e−(x+1)) .
3. Letting x = 2 and solving 1/3 = 1/(4 + c) we get c = −1. The solution is y = 1/(x2 − 1).
This solution is defined on the interval (1, ∞).
4. Letting x = −2 and solving 1/2 = 1/(4 + c) we get c = −2. The solution is y = 1/(x2 − 2).
This solution is defined on the interval (−∞, −
√
2 ).
5. Letting x = 0 and solving 1 = 1/c we get c = 1. The solution is y = 1/(x2 + 1). This solution
is defined on the interval (−∞, ∞).
6. Letting x = 1/2 and solving −4 = 1/(1/4 + c) we get c = −1/2. The solution is y =
1/(x2 − 1/2) = 2/(2x2 − 1). This solution is defined on the interval (−1/
√
2 , 1/
√
2 ).
In Problems 7–10, we use x = c1 cos t + c2 sin t and x = −c1 sin t + c2 cos t to obtain a system of
two equations in the two unknowns c1 and c2.
1.2 Initial-Value Problems 15
7. From the initial conditions we obtain the system
c1 = −1c2 = 8
The solution of the initial-value problem is x = − cos t + 8 sin t.
8. From the initial conditions we obtain the system
c2 = 0 − c1 = 1
The solution of the initial-value problem is x = − cos t.
9. From the initial conditions we obtain
√
3
2
c1 +
1
2
c2 =
1
2
−
1
2
c2 +
√
3
2
= 0
Solving, we find c1 =
√
3/4 and c2 = 1/4. The solution of the initial-value problem is
x = (
√
3/4) cos t + (1/4) sin t.
10. From the initial conditions we obtain
√
2
2
c1 +
√
2
2
c2 =
√
2
−
√
2
2
c1 +
√
2
2
c2 = 2
√
2.
Solving, we find c1 = −1 and c2 = 3. The solution of the initial-value problem is x =
− cos t + 3 sin t.
In Problems 11–14, we use y = c1ex + c2e−x and y = c1ex − c2e−x to obtain a system of two
equations in the two unknowns c1 and c2.
11. From the initial conditions we obtain
c1 + c2 = 1
c1 − c2 = 2.
Solving, we find c1 = 3
2 and c2 = −1
2 . The solution of the initial-value problem is y =
3
2ex − 1
2e−x.
12. From the initial conditions we obtain
ec1 + e−1
c2 = 0
ec1 − e−1
c2 = e.
Solving, we find c1 = 1
2 and c2 = −1
2e2. The solution of the initial-value problem is
y =
1
2
ex
−
1
2
e2
e−x
=
1
2
ex
−
1
2
e2−x
.
16 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS
13. From the initial conditions we obtain
e−1
c1 + ec2 = 5
e−1
c1 − ec2 = −5.
Solving, we find c1 = 0 and c2 = 5e−1. The solution of the initial-value problem is y =
5e−1e−x = 5e−1−x.
14. From the initial conditions we obtain
c1 + c2 = 0
c1 − c2 = 0.
Solving, we find c1 = c2 = 0. The solution of the initial-value problem is y = 0.
15. Two solutions are y = 0 and y = x3.
16. Two solutions are y = 0 and y = x2. (Also, any constant multiple of x2 is a solution.)
17. For f(x, y) = y2/3 we have
∂f
∂y
=
2
3
y−1/3
. Thus, the differential equation will have a unique
solution in any rectangular region of the plane where y = 0.
18. For f(x, y) =
√
xy we have ∂f/∂y = 1
2

x/y . Thus the differential equation will have a
unique solution in any region where x  0 and y  0 or where x  0 and y  0.
19. For f(x, y) =
y
x
we have
∂f
∂y
=
1
x
. Thus, the differential equation will have a unique solution
in any region where x = 0.
20. For f(x, y) = x + y we have
∂f
∂y
= 1. Thus, the differential equation will have a unique
solution in the entire plane.
21. For f(x, y) = x2/(4 − y2) we have ∂f/∂y = 2x2y/(4 − y2)2. Thus the differential equation
will have a unique solution in any region where y  −2, −2  y  2, or y  2.
22. For f(x, y) =
x2
1 + y3
we have
∂f
∂y
=
−3x2y2
(1 + y3)2 . Thus, the differential equation will have a
unique solution in any region where y = −1.
23. For f(x, y) =
y2
x2 + y2
we have
∂f
∂y
=
2x2y
(x2 + y2)2 . Thus, the differential equation will have a
unique solution in any region not containing (0, 0).
1.2 Initial-Value Problems 17
24. For f(x, y) = (y + x)/(y − x) we have ∂f/∂y = −2x/(y − x)2. Thus the differential equation
will have a unique solution in any region where y  x or where y  x.
In Problems 25–28, we identify f(x, y) =

y2 − 9 and ∂f/∂y = y/

y2 − 9. We see that f and
∂f/∂y are both continuous in the regions of the plane determined by y  −3 and y  3 with no
restrictions on x.
25. Since 4  3, (1, 4) is in the region defined by y  3 and the differential equation has a unique
solution through (1, 4).
26. Since (5, 3) is not in either of the regions defined by y  −3 or y  3, there is no guarantee
of a unique solution through (5, 3).
27. Since (2, −3) is not in either of the regions defined by y  −3 or y  3, there is no guarantee
of a unique solution through (2, −3).
28. Since (−1, 1) is not in either of the regions defined by y  −3 or y  3, there is no guarantee
of a unique solution through (−1, 1).
29. (a) A one-parameter family of solutions is y = cx. Since y = c, xy = xc = y and y(0) =
c · 0 = 0.
(b) Writing the equation in the form y = y/x, we see that R cannot contain any point on the
y-axis. Thus, any rectangular region disjoint from the y-axis and containing (x0, y0) will
determine an interval around x0 and a unique solution through (x0, y0). Since x0 = 0 in
part (a), we are not guaranteed a unique solution through (0, 0).
(c) The piecewise-defined function which satisfies y(0) = 0 is not a solution since it is not
differentiable at x = 0.
30. (a) Since
d
dx
tan (x + c) = sec2
(x + c) = 1+tan2
(x + c), we see that y = tan (x + c) satisfies
the differential equation.
(b) Solving y(0) = tan c = 0 we obtain c = 0 and y = tan x. Since tan x is discontinuous at
x = ±π/2, the solution is not defined on (−2, 2) because it contains ±π/2.
(c) The largest interval on which the solution can exist is (−π/2, π/2).
31. (a) Since
d
dx

−
1
x + c

=
1
(x + c)2
= y2
, we see that y = −
1
x + c
is a solution of the differ-
ential equation.
(b) Solving y(0) = −1/c = 1 we obtain c = −1 and y = 1/(1−x). Solving y(0) = −1/c = −1
we obtain c = 1 and y = −1/(1+x). Being sure to include x = 0, we see that the interval
of existence of y = 1/(1−x) is (−∞, 1), while the interval of existence of y = −1/(1+x)
is (−1, ∞).
18 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS
32. (a) Solving y(0) = −1/c = y0 we obtain c = −1/y0 and
y = −
1
−1/y0 + x
= y01 − y0x , y0 = 0
Since we must have −1/y0 + x = 0, the largest interval of existence (which must contain
0) is either (−∞, 1/y0) when y0  0 or (1/y0, ∞) when y0  0.
(b) By inspection we see that y = 0 is a solution on (−∞, ∞).
33. (a) Differentiating 3x2 − y2 = c we get 6x − 2yy = 0 or yy = 3x.
(b) Solving 3x2 − y2 = 3 for y we get
y = φ1(x) =

3(x2 − 1) , 1  x  ∞,
y = φ2(x) = −

3(x2 − 1) , 1  x  ∞,
y = φ3(x) =

3(x2 − 1) , −∞  x  −1,
y = φ4(x) = −

3(x2 − 1) , −∞  x  −1.
–4 –2 2 4
x
–4
–2
2
4
y
(c) Only y = φ3(x) satisfies y(−2) = 3.
34. (a) Setting x = 2 and y = −4 in 3x2 − y2 = c we get
12 − 16 = −4 = c, so the explicit solution is
y = −

3x2 + 4 , −∞  x  ∞.
–4 –2 2 4
x
–4
–2
2
4
y
(b) Setting c = 0 we have y =
√
3x and y = −
√
3x, both defined on (−∞, ∞).
In Problems 35–38, we consider the points on the graphs with x-coordinates x0 = −1, x0 = 0,
and x0 = 1. The slopes of the tangent lines at these points are compared with the slopes given by
y(x0) in (a) through (f).
35. The graph satisfies the conditions in (b) and (f).
36. The graph satisfies the conditions in (e).
37. The graph satisfies the conditions in (c) and (d).
1.2 Initial-Value Problems 19
38. The graph satisfies the conditions in (a).
39. Using the function y = c1 cos 3x + c2 sin 3x and the first boundary condition we get
y(0) = c1 cos 0 + c2 sin 0 = 0
Therefore c1 = 0. Similarly, for the second boundary condition we get
y(π/6) = c2 sin 3 (π/6) = −1
Therefore c2 = −1. The solution to the boundary value problem is y(x) = − sin 3x.
40. Using the function y = c1 cos 3x + c2 sin 3x and the first boundary condition we get
y(0) = c1 cos 0 + c2 sin 0 = 0
Therefore c1 = 0. Thus y(x) = c2 sin 3x. Similarly for the second boundary condition we get
y(π) = c2 sin 3(π) = 0
But the last line is satisfied for any choice of c2 since sin 3π = 0 therefore there are infinitely
many solutions in this case.
41. The derivative of the function y = c1 cos 3x + c2 sin 3x is y = −3c1 sin 3x + 3c2 cos 3x and
using the first boundary condition we get
y
(0) = 0 + 3c2 = 0
Therefore c2 = 0 So far then we have y = −3c1 sin 3x. Similarly, using the second boundary
condition we get
y
(π/4) = −3c1
√
2/2

= 0
Therefore c1 = 0. The solution is thus the trivial solution y = 0.
42. The derivative of the function y = c1 cos 3x + c2 sin 3x is y = −3c1 sin 3x + 3c2 cos 3x and
using the two boundary conditions we get
y(0) = c1 + 0 = 1
Therefore c1 = 1. In addition
y
(π) = 0 − 3c2 = 5
Therefore c2 = −5/3. The solution to this boundary value problem is y(x) = cos 3x−
5
3
sin 3x.
20 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS
43. Using the function y = c1 cos 3x + c2 sin 3x and the first boundary condition we get
y(0) = c1 cos 0 + c2 sin 0 = 0
Therefore c1 = 0. Thus y(x) = c2 sin 3x. Similarly for the second boundary condition we get
c2 sin 3(π) = 4c2 · 0 = 40 = 4
The last line is obviously a contradiction and so therefore the boundary value problem has
no solution in this case.
44. The derivative of the function y = c1 cos 3x + c2 sin 3x is y = −3c1 sin 3x + 3c2 cos 3x and
using the first boundary condition we get
y
(π/3) = 0 − 3c2 = 1
Therefore c2 = −1/3. Similarly, using the second boundary condition we get
y
(π) = 0 − 3c2 = 0
Therefore c2 = 0, which is a contradiction so the problem has no solution.
45. Integrating y = 8e2x + 6x we obtain
y = (8e2x
+ 6x) dx = 4e2x
+ 3x2
+ c.
Setting x = 0 and y = 9 we have 9 = 4 + c so c = 5 and y = 4e2x + 3x2 + 5.
46. Integrating y = 12x − 2 we obtain
y
= (12x − 2) dx = 6x2
− 2x + c1.
Then, integrating y we obtain
y = (6x2
− 2x + c1) dx = 2x3
− x2
+ c1x + c2.
At x = 1 the y-coordinate of the point of tangency is y = −1 + 5 = 4. This gives the initial
condition y(1) = 4. The slope of the tangent line at x = 1 is y(1) = −1. From the initial
conditions we obtain
2 − 1 + c1 + c2 = 4 or c1 + c2 = 3
and
6 − 2 + c1 = −1 or c1 = −5.
Thus, c1 = −5 and c2 = 8, so y = 2x3 − x2 − 5x + 8.
1.2 Initial-Value Problems 21
47. When x = 0 and y = 1
2 , y = −1, so the only plausible solution curve is the one with negative
slope at (0, 1
2 ), or the blue curve.
48. If y(x) is a solution of the initial-value problem y = x2+y2. y(0) = 1, then y(0) = 02+12 = 1.
Differentiating the differential equation, we have y = 2x+2yy. So y(0) = 2·0+2·1·1 = 2.
49. If the solution is tangent to the x-axis at (x0, 0), then y = 0 when x = x0 and y = 0.
Substituting these values into y + 2y = 3x − 6 we get 0 + 0 = 3x0 − 6 or x0 = 2.
50. When y = 1
16 x4, y = 1
4x3 = x(1
4 x2) = xy1/2, and y(2) = 1
16(16) = 1. When
y =
⎧
⎪
⎨
⎪
⎩
0, x  0
1
16
x4, x ≥ 0
we have
y
=
⎧
⎨
⎩
0, x  0
1
4x3, x ≥ 0
= x
⎧
⎨
⎩
0, x  0
1
4x2, x ≥ 0
= xy1/2
,
and y(2) = 1
16 (16) = 1. The two different solutions are the same on the interval (0, ∞), which
is all that is required by Theorem 1.2.1.
51. We note that the initial condition y(0) = 0,
0 =
y
0
1
√
t3 + 1
dt
is satisfied only when y = 0. For any y  0, necessarily
y
0
1
√
t3 + 1
dt  0
because the integrand is positive on the interval of integration. Then from (12) of Section 1.1
and the Chain Rule we have:
d
dx
x =
d
dx
y
0
1
√
t3 + 1
dt
1 =
1

y3 + 1
dy
dx
and
dy
dx
=

y3 + 1
y
(0) =
dy
dx



x=0
=

(y(0))3
+ 1 =
√
0 + 1 = 1.
Computing the second derivative, we see that:
d2y
dx2
=
d
dx

y3 + 1 =
3y2
2

y3 + 1
dy
dx
=
3y2
2

y3 + 1
·

y3 + 1 =
3
2
y2
d2y
dx2
=
3
2
y2
.
This is equivalent to 2
d2y
dx2
− 3y2
= 0.
22 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS
1.3 Differential Equations as Mathematical Models
1.
dP
dt
= kP + r;
dP
dt
= kP − r
2. Let b be the rate of births and d the rate of deaths. Then b = k1P and d = k2P. Since
dP/dt = b − d, the differential equation is dP/dt = k1P − k2P.
3. Let b be the rate of births and d the rate of deaths. Then b = k1P and d = k2P2. Since
dP/dt = b − d, the differential equation is dP/dt = k1P − k2P2.
4.
dP
dt
= k1P − k2P2
− h, h  0
5. From the graph in the text we estimate T0 = 180◦ and Tm = 75◦. We observe that when
T = 85, dT/dt ≈ −1. From the differential equation we then have
k =
dT/dt
T − Tm
=
−1
85 − 75
= −0.1.
6. By inspecting the graph in the text we take Tm to be Tm(t) = 80 − 30 cos (πt/12). Then the
temperature of the body at time t is determined by the differential equation
dT
dt
= k

T −

80 − 30 cos
 π
12
t

, t  0.
7. The number of students with the flu is x and the number not infected is 1000 − x, so
dx/dt = kx(1000 − x).
8. By analogy, with the differential equation modeling the spread of a disease, we assume that
the rate at which the technological innovation is adopted is proportional to the number of
people who have adopted the innovation and also to the number of people, y(t), who have
not yet adopted it. If one person who has adopted the innovation is introduced into the
population, then x + y = n + 1 and
dx
dt
= kx(n + 1 − x), x(0) = 1.
9. The rate at which salt is leaving the tank is
Rout (3 gal/min) ·

A
300
lb/gal

=
A
100
lb/min.
Thus dA/dt = A/100. The initial amount is A(0) = 50.
1.3 Differential Equations as Mathematical Models 23
10. The rate at which salt is entering the tank is
Rin = (3 gal/min) · (2 lb/gal) = 6 lb/min.
Since the solution is pumped out at a slower rate, it is accumulating at the rate of (3 − 2)
gal/min = 1 gal/min. After t minutes there are 300 + t gallons of brine in the tank. The rate
at which salt is leaving is
Rout = (2 gal/min) ·

A
300 + t
lb/gal

=
2A
300 + t
lb/min.
The differential equation is
dA
dt
= 6 −
2A
300 + t
.
11. The rate at which salt is entering the tank is
Rin = (3 gal/min) · (2 lb/gal) = 6 lb/min.
Since the tank loses liquid at the net rate of
3 gal/min − 3.5 gal/min = −0.5 gal/min,
after t minutes the number of gallons of brine in the tank is 300 − 1
2t gallons. Thus the rate
at which salt is leaving is
Rout =

A
300 − t/2
lb/gal

· (3.5 gal/min) =
3.5A
300 − t/2
lb/min =
7A
600 − t
lb/min.
The differential equation is
dA
dt
= 6 −
7A
600 − t
or
dA
dt
+
7
600 − t
A = 6.
12. The rate at which salt is entering the tank is
Rin = (cin lb/gal) · (rin gal/min) = cinrin lb/min.
Now let A(t) denote the number of pounds of salt and N(t) the number of gallons of brine
in the tank at time t. The concentration of salt in the tank as well as in the outflow is
c(t) = x(t)/N(t). But the number of gallons of brine in the tank remains steady, is increased,
or is decreased depending on whether rin = rout, rin  rout, or rin  rout. In any case, the
number of gallons of brine in the tank at time t is N(t) = N0 + (rin − rout)t. The output rate
of salt is then
Rout =

A
N0 + (rin − rout)t
lb/gal

· (rout gal/min) = rout
A
N0 + (rin − rout)t
lb/min.
The differential equation for the amount of salt, dA/dt = Rin − Rout, is
dA
dt
= cinrin − rout
A
N0 + (rin − rout)t
or
dA
dt
+
rout
N0 + (rin − rout)t
A = cinrin.
24 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS
13. The volume of water in the tank at time t is V = Awh. The differential equation is then
dh
dt
=
1
Aw
dV
dt
=
1
Aw

−cAh

2gh

= −
cAh
Aw

2gh .
Using Ah = π

2
12
2
=
π
36
, Aw = 102 = 100, and g = 32, this becomes
dh
dt
= −
cπ/36
100
√
64h = −
cπ
450
√
h .
14. The volume of water in the tank at time t is V = 1
3 πr2h where r is the radius of the tank
at height h. From the figure in the text we see that r/h = 8/20 so that r = 2
5h and
V = 1
3π
2
5h
2
h = 4
75 πh3. Differentiating with respect to t we have dV/dt = 4
25πh2 dh/dt or
dh
dt
=
25
4πh2
dV
dt
.
From Problem 13 we have dV/dt = −cAh
√
2gh where c = 0.6, Ah = π
 2
12
2
, and g = 32.
Thus dV/dt = −2π
√
h/15 and
dh
dt
=
25
4πh2

−
2π
√
h
15

= −
5
6h3/2
.
15. Since i = dq/dt and L d2q/dt2 + R dq/dt = E(t), we obtain L di/dt + Ri = E(t).
16. By Kirchhoff’s second law we obtain R
dq
dt
+
1
C
q = E(t).
17. From Newton’s second law we obtain m
dv
dt
= −kv2
+ mg.
18. Since the barrel in Figure 1.3.17(b) in the text is submerged an additional y feet below its
equilibrium position, the number of cubic feet in the additional submerged portion is the
volume of the circular cylinder: π × (radius)2
× height or π(s/2)2y. Then we have from
Archimedes’ principle
Upward force of water on barrel = Weight of water displaced
= (62.4) × (Volume of water displaced)
= (62.4)π(s/2)2
y = 15.6πs2
y.
It then follows from Newton’s second law that
w
g
d2y
dt2
= −15.6πs2
y or
d2y
dt2
+
15.6πs2g
w
y = 0,
where g = 32 and w is the weight of the barrel in pounds.
1.3 Differential Equations as Mathematical Models 25
19. The net force acting on the mass is
F = ma = m
d2x
dt2
= −k(s + x) + mg = −kx + mg − ks.
Since the condition of equilibrium is mg = ks, the differential equation is
m
d2x
dt2
= −kx.
20. From Problem 19, without a damping force, the differential equation is m d2x/dt2 = −kx.
With a damping force proportional to velocity, the differential equation becomes
m
d2x
dt2
= −kx − β
dx
dt
or m
d2x
dt2
+ β
dx
dt
+ kx = 0.
21. As the rocket climbs (in the positive direction), it spends its amount of fuel and therefore the
mass of the fuel changes with time. The air resistance acts in the opposite direction of the
motion and the upward thrust R works in the same direction. Using Newton’s second law we
get
d
dt
(mv) = −mg − kv + R
Now because the mass is variable, we must use the product rule to expand the left side of the
equation. Doing so gives us the following:
d
dt
(mv) = −mg − kv + R
v ×
dm
dt
+ m ×
dv
dt
= −mg − kv + R
The last line is the differential equation we wanted to find.
22. (a) Since the mass of the rocket is m(t) = mp + mv + mf (t), take the time rate-of-change
and get, by straight-forward calculation,
d
dt
m(t) =
d
dt
(mp + mv + mf (t)) = 0 + 0 + m
f (t) =
d
dt
mf (t)
Therefore the rate of change of the mass of the rocket is the same as the rate of change
of the mass of the fuel which is what we wanted to show.
(b) The fuel is decreasing at the constant rate of λ and so from part (a) we have
d
dt
m(t) =
d
dt
mf (t) = −λ
m(t) = −λt + c
26 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS
Using the given condition to solve for c, m(0) = 0 + c = m0 and so m(t) = −λt + m0.
The differential equation in Problem 21 now becomes
v
dm
dt
+ m
dv
dt
+ kv = −mg + R
−λv + (−λt + m0)
dv
dt
+ kv = −mg + R
(−λt + m0)
dv
dt
+ (k − λ)v = −mg + R
dv
dt
+
k − λ
−λt + m0
v =
−mg
−λt + m0
+
R
−λt + m0
dv
dt
+
k − λ
−λt + m0
v = −g +
R
−λt + m0
(c) From part (b) we have that d
dtmf (t) = −λ and so by integrating this result we get
mf (t) = −λt+c. Now at time t = 0, mf (0) = 0+c = c therefore mf (t) = −λt+mf (0) .
At some later time tb we then have mf (tb) = −λtb +mf (0) = 0 and solving this equation
for that time we get tb = mf (0)/λ which is what we wanted to show.
23. From g = k/R2 we find k = gR2. Using a = d2r/dt2 and the fact that the positive direction
is upward we get
d2r
dt2
= −a = −
k
r2
= −
gR2
r2
or
d2r
dt2
+
gR2
r2
= 0.
24. The gravitational force on m is F = −kMrm/r2. Since Mr = 4πδr3/3 and M = 4πδR3/3,
we have Mr = r3M/R3 and
F = −k
Mrm
r2
= −k
r3Mm/R3
r2
= −k
mM
R3
r.
Now from F = ma = d2r/dt2 we have
m
d2r
dt2
= −k
mM
R3
r or
d2r
dt2
= −
kM
R3
r.
25. The differential equation is
dA
dt
= k(M − A).
26. The differential equation is
dA
dt
= k1(M − A) − k2A.
27. The differential equation is x(t) = r − kx(t) where k  0.
1.3 Differential Equations as Mathematical Models 27
28. Notice from the diagram that the segment from the waterskier to the boat is tangent to the
curve at the point P(x, y). For the sake of simplicity, let’s label the coordinates of the boat
on the x-axis as (a, 0). The equation of this tangent line is y = y(x)(x − a) from which we
get
y = y
(x)(x − a)
y
y
= x − a
−
y
y
= a − x
By the Pythagorean Theorem we have y2 + (a − x)2 = s2 from which we make a substitution
and then solve for the derivative to get
y2
+ (a − x)2
= s2
y2
+ (y

y
)2
= s2
(y

y
)2
= s2
− y2
y
y
= −

s2 − y2
y
= −
y

s2 − y2
Notice that the sign of the derivative is negative because as the boat proceeds along the
positive x-axis, the y-coordinate decreases.
29. We see from the figure that 2θ + α = π. Thus
y
−x
= tan α = tan(π − 2θ) = − tan 2θ = −
2 tan θ
1 − tan2 θ
.
Since the slope of the tangent line is y = tan θ we have
y/x = 2y[1−(y)2] or y−y(y)2 = 2xy, which is the quadratic
equation y(y)2 + 2xy − y = 0 in y. Using the quadratic
formula, we get
y
=
−2x ±

4x2 + 4y2
2y
=
−x ±

x2 + y2
y
.
(x, y)
x
y
α
α
θ
θ
θ
φ
x
y
Since dy/dx  0, the differential equation is
dy
dx
=
−x +

x2 + y2
y
or y
dy
dx
−

x2 + y2 + x = 0.
28 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS
30. The differential equation is dP/dt = kP, so from Problem 49 in Exercises 1.1, a one-parameter
family of solutions is P = cekt.
31. The differential equation in (3) is dT/dt = k(T − Tm). When the body is cooling, T  Tm,
so T − Tm  0. Since T is decreasing, dT/dt  0 and k  0. When the body is warming,
T  Tm, so T − Tm  0. Since T is increasing, dT/dt  0 and k  0.
32. The differential equation in (8) is dA/dt = 6 − A/100. If A(t) attains a maximum, then
dA/dt = 0 at this time and A = 600. If A(t) continues to increase without reaching a
maximum, then A(t)  0 for t  0 and A cannot exceed 600. In this case, if A(t) approaches
0 as t increases to infinity, we see that A(t) approaches 600 as t increases to infinity.
33. This differential equation could describe a population that undergoes periodic fluctuations.
34. (a) As shown in Figure 1.3.23(a) in the text, the resultant of the reaction force of magnitude
F and the weight of magnitude mg of the particle is the centripetal force of magnitude
mω2x. The centripetal force points to the center of the circle of radius x on which the
particle rotates about the y-axis. Comparing parts of similar triangles gives
F cos θ = mg and F sin θ = mω2
x.
(b) Using the equations in part (a) we find
tan θ =
F sin θ
F cos θ
=
mω2x
mg
=
ω2x
g
or
dy
dx
=
ω2x
g
.
35. From Problem 23, d2r/dt2 = −gR2/r2. Since R is a constant, if r = R + s, then d2r/dt2 =
d2s/dt2 and, using a Taylor series, we get
d2s
dt2
= −g
R2
(R + s)2
= −gR2
(R + s)−2
≈ −gR2
[R−2
− 2sR−3
+ · · · ] = −g +
2gs
R
+ · · · .
Thus, for R much larger than s, the differential equation is approximated by d2s/dt2 = −g.
36. (a) If ρ is the mass density of the raindrop, then m = ρV and
dm
dt
= ρ
dV
dt
= ρ
d
dt
4
3
πr3

= ρ

4πr2 dr
dt

= ρS
dr
dt
.
If dr/dt is a constant, then dm/dt = kS where ρ dr/dt = k or dr/dt = k/ρ. Since the
radius is decreasing, k  0. Solving dr/dt = k/ρ we get r = (k/ρ)t+c0. Since r(0) = r0,
c0 = r0 and r = kt/ρ + r0.
(b) From Newton’s second law,
d
dt
[mv] = mg, where v is the velocity of the raindrop. Then
m
dv
dt
+ v
dm
dt
= mg or ρ
4
3
πr3
dv
dt
+ v(k4πr2
) = ρ
4
3
πr3

g.
Dividing by 4ρπr3/3 we get
dv
dt
+
3k
ρr
v = g or
dv
dt
+
3k/ρ
kt/ρ + r0
v = g, k  0.
1.3 Differential Equations as Mathematical Models 29
37. We assume that the plow clears snow at a constant rate of k cubic miles per hour. Let t be the
time in hours after noon, x(t) the depth in miles of the snow at time t, and y(t) the distance
the plow has moved in t hours. Then dy/dt is the velocity of the plow and the assumption
gives
wx
dy
dt
= k,
where w is the width of the plow. Each side of this equation simply represents the volume
of snow plowed in one hour. Now let t0 be the number of hours before noon when it started
snowing and let s be the constant rate in miles per hour at which x increases. Then for
t  −t0, x = s(t + t0). The differential equation then becomes
dy
dt
=
k
ws
1
t + t0
.
Integrating, we obtain
y =
k
ws
[ ln(t + t0) + c ]
where c is a constant. Now when t = 0, y = 0 so c = − ln t0 and
y =
k
ws
ln

1 +
t
t0

.
Finally, from the fact that when t = 1, y = 2 and when t = 2, y = 3, we obtain

1 +
2
t0
2
=

1 +
1
t0
3
.
Expanding and simplifying gives t2
0 + t0 − 1 = 0. Since t0  0, we find t0 ≈ 0.618 hours ≈ 37
minutes. Thus it started snowing at about 11:23 in the morning.
38. . (1) :
dP
dt
= kP is linear (2) :
dA
dt
= kA is linear
(3) :
dT
dt
= k (T − Tm) is linear (5) :
dx
dt
= kx (n + 1 − x) is nonlinear
(6) :
dX
dt
= k (α − X) (β − X) is nonlinear (8) :
dA
dt
= 6 −
A
100
is linear
(10) :
dh
dt
= −
Ah
Aw

2gh is nonlinear (11) : L
d2q
dt2
+ R
dq
dt
+
1
C
q = E(t) is linear
(12) :
d2s
dt2
= −g is linear (14) : m
dv
dt
= mg − kv is linear
(15) : m
d2s
dt2
+ k
ds
dt
= mg is linear (16) :
d2x
dt2
−
64
L
x = 0 is linear
(17) : linearity or nonlinearity is determined by the manner in which W and T1 involve x.
39. At time t, when the population is 2 million cells, the differential equation P(t) = 0.15P(t)
gives the rate of increase at time t. Thus, when P(t) = 2 (million cells), the rate of increase
is P(t) = 0.15(2) = 0.3 million cells per hour or 300,000 cells per hour.
30 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS
40. Setting A(t) = −0.002 and solving A(t) = −0.0004332A(t) for A(t), we obtain
A(t) =
A(t)
−0.0004332
=
−0.002
−0.0004332
≈ 4.6 grams.
Chapter 1 in Review
1.
d
dx
c1ekx
= k
y
  
c1ekx
;
dy
dx
= ky
2.
d
dx
(5 + c1e−2x
) = −2c1e−2x
= −2(
y
  
5 + c1e−2x
−5);
dy
dx
= −2(y − 5) or
dy
dx
= −2y + 10
3.
d
dx
(c1 cos kx + c2 sin kx) = −kc1 sin kx + kc2 cos kx;
d2
dx2
(c1 cos kx + c2 sin kx) = −k2
c1 cos kx − k2
c2 sin kx = −k2
(
y
  
c1 cos kx + c2 sin kx);
d2y
dx2
= −k2
y or
d2y
dx2
+ k2
y = 0
4.
d
dx
(c1 cosh kx + c2 sinh kx) = kc1 sinh kx + kc2 cosh kx;
d2
dx2
(c1 cosh kx + c2 sinh kx) = k2
c1 cosh kx + k2
c2 sinh kx = k2
(
y
  
c1 cosh kx + c2 sinh kx);
d2y
dx2
= k2
y or
d2y
dx2
− k2
y = 0
5. y = c1ex + c2xex; y = c1ex + c2xex + c2ex; y = c1ex + c2xex + 2c2ex;
y + y = 2(c1ex + c2xex) + 2c2ex = 2(c1ex + c2xex + c2ex) = 2y; y − 2y + y = 0
6. y = −c1ex sin x + c1ex cos x + c2ex cos x + c2ex sin x;
y = −c1ex cos x − c1ex sin x − c1ex sin x + c1ex cos x − c2ex sin x + c2ex cos x + c2ex cos x +
c2ex sin x
= −2c1ex sin x + 2c2ex cos x;
y − 2y = −2c1ex cos x − 2c2ex sin x = −2y; y − 2y + 2y = 0
7. a, d 8. c 9. b 10. a, c 11. b 12. a, b, d
13. A few solutions are y = 0, y = c, and y = ex.
14. Easy solutions to see are y = 0 and y = 3.
15. The slope of the tangent line at (x, y) is y, so the differential equation is y = x2 + y2.
Chapter 1 in Review 31
16. The rate at which the slope changes is dy/dx = y, so the differential equation is y = −y
or y + y = 0.
17. (a) The domain is all real numbers.
(b) Since y = 2/3x1/3, the solution y = x2/3 is undefined at x = 0. This function is a
solution of the differential equation on (−∞, 0) and also on (0, ∞).
18. (a) Differentiating y2 − 2y = x2 − x + c we obtain 2yy − 2y = 2x − 1 or (2y − 2)y = 2x − 1.
(b) Setting x = 0 and y = 1 in the solution we have 1 − 2 = 0 − 0 + c or c = −1. Thus, a
solution of the initial-value problem is y2 − 2y = x2 − x − 1.
(c) Solving the equation y2 − 2y − (x2 − x − 1) = 0 by the quadratic formula we get
y = (2 ±

4 + 4(x2 − x − 1) )/2 = 1 ±
√
x2 − x = 1 ±

x(x − 1) . Since x(x − 1) ≥ 0
for x ≤ 0 or x ≥ 1, we see that neither y = 1 +

x(x − 1) nor y = 1 −

x(x − 1) is
differentiable at x = 0. Thus, both functions are solutions of the differential equation,
but neither is a solution of the initial-value problem.
19. Setting x = x0 and y = 1 in y = −2/x + x, we get
1 = −
2
x0
+ x0 or x2
0 − x0 − 2 = (x0 − 2)(x0 + 1) = 0.
Thus, x0 = 2 or x0 = −1. Since x = 0 in y = −2/x + x, we see that y = −2/x + x is a
solution of the initial-value problem xy + y = 2x, y(−1) = 1, on the interval (−∞, 0) and
y = −2/x + x is a solution of the initial-value problem xy + y = 2x, y(2) = 1, on the interval
(0, ∞).
20. From the differential equation, y(1) = 12+[y(1)]2 = 1+(−1)2 = 2  0, so y(x) is increasing in
some neighborhood of x = 1. From y = 2x+2yy we have y(1) = 2(1)+2(−1)(2) = −2  0,
so y(x) is concave down in some neighborhood of x = 1.
21. (a)
–3 –2 1 2 3
x
–3
–2
1
2
3
y
–1
–1 –3 –2 1 2 3
x
–3
–2
1
y = x2 + c1 y = –x2 + c2
2
3
y
–1
–1
(b) When y = x2 + c1, y = 2x and (y)2 = 4x2. When y = −x2 + c2, y = −2x and
(y)2 = 4x2.
32 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS
(c) Pasting together x2, x ≥ 0, and −x2, x ≤ 0, we get y =
⎧
⎨
⎩
−x2, x ≤ 0
x2, x  0 .
22. The slope of the tangent line is y


(−1,4)
= 6
√
4 + 5(−1)3 = 7.
23. Differentiating y = x sin x + x cos x we get
y
= x cos x + sin x − x sin x + cos x
and
y
= −x sin x + cos x + cos x − x cos x − sin x − sin x
= −x sin x − x cos x + 2 cos x − 2 sin x.
Thus
y
+ y = −x sin x − x cos x + 2 cos x − 2 sin x + x sin x + x cos x = 2 cos x − 2 sin x.
An interval of definition for the solution is (−∞, ∞).
24. Differentiating y = x sin x + (cos x) ln(cos x) we get
y
= x cos x + sin x + cos x

− sin x
cos x

− (sin x) ln (cos x)
= x cos x + sin x − sin x − (sin x) ln (cos x)
= x cos x − (sin x) ln (cos x)
and
y
= −x sin x + cos x − sin x

− sin x
cos x

− (cos x) ln (cos x)
= −x sin x + cos x +
sin2
x
cos x
− (cos x) ln (cos x)
= −x sin x + cos x +
1 − cos2 x
cos x
− (cos x) ln (cos x)
= −x sin x + cos x + sec x − cos x − (cos x) ln (cos x)
= −x sin x + sec x − (cos x) ln (cos x).
Thus
y
+ y = −x sin x + sec x − (cos x) ln(cos x) + x sin x + (cos x) ln (cos x) = sec x.
To obtain an interval of definition we note that the domain of ln x is (0, ∞), so we must have
cos x  0. Thus, an interval of definition is (−π/2, π/2).
Chapter 1 in Review 33
25. Differentiating y = sin (ln x) we obtain y = cos (ln x)/x and y = −[sin (ln x) + cos (ln x)]/x2.
Then
x2
y
+ xy
+ y = x2

−
sin (ln x) + cos (ln x)
x2

+ x
cos (ln x)
x
+ sin (ln x) = 0.
An interval of definition for the solution is (0, ∞).
26. Differentiating y = cos (ln x) ln (cos (ln x)) + (ln x) sin (ln x) we obtain
y
= cos (ln x)
1
cos (ln x)

−
sin (ln x)
x

+ ln (cos (ln x))

−
sin (ln x)
x

+ ln x
cos (ln x)
x
+
sin (ln x)
x
= −
ln (cos (ln x)) sin (ln x)
x
+
(ln x) cos (ln x)
x
and
y
= −x ln (cos (ln x))
cos (ln x)
x
+ sin (ln x)
1
cos (ln x)

−
sin (ln x)
x

1
x2
+ ln (cos (ln x)) sin (ln x)
1
x2
+ x (ln x)

−
sin (ln x)
x

+
cos (ln x)
x
1
x2
− (ln x) cos (ln x)
1
x2
=
1
x2
− ln (cos (ln x)) cos (ln x) +
sin2
(ln x)
cos (ln x)
+ ln (cos (ln x)) sin (ln x)
− (ln x) sin (ln x) + cos (ln x) − (ln x) cos (ln x) .
Then
x2
y
+ xy
+ y = − ln (cos (ln x)) cos (ln x) +
sin2
(ln x)
cos (ln x)
+ ln (cos (ln x)) sin (ln x)
− (ln x) sin (ln x) + cos (ln x) − (ln x) cos (ln x) − ln (cos (ln x)) sin (ln x)
+ (ln x) cos (ln x) + cos (ln x) ln (cos (ln x)) + (ln x) sin (ln x)
=
sin2
(ln x)
cos (ln x)
+ cos (ln x) =
sin2
(ln x) + cos2 (ln x)
cos (ln x)
=
1
cos (ln x)
= sec (ln x).
To obtain an interval of definition, we note that the domain of ln x is (0, ∞), so we must
have cos (ln x)  0. Since cos x  0 when −π/2  x  π/2, we require −π/2  ln x  π/2.
Since ex is an increasing function, this is equivalent to e−π/2  x  eπ/2. Thus, an interval
of definition is (e−π/2, eπ/2). (Much of this problem is more easily done using a computer
algebra system such as Mathematica or Maple.)
34 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS
In Problems 27–30 we use (12) of Section 1.1 and the Product Rule.
27. .
y = ecos x
x
0
te− cos t
dt
dy
dx
= ecos x

xe− cos x

− sin xecos x
x
0
te− cos t
dt
dy
dx
+ (sin x) y = ecos x
xe− cos x
− sin xecos x
x
0
te− cos t
dt + sin x

ecos x
x
0
te− cos t
dt

= x − sin xecos x
x
0
te− cos t
dt + sin xecos x
x
0
te− cos t
dt = x
28. .
y = ex2
x
0
et−t2
dt
dy
dx
= ex2
ex−x2
+ 2xex2
x
0
et−t2
dt
dy
dx
− 2xy = ex2
ex−x2
+ 2xex2
x
0
et−t2
dt − 2x

ex2
x
0
et−t2
dt

= ex
29. .
y = x
x
1
e−t
t
dt
y
= x
e−x
x
+
x
1
e−t
t
dt = e−x
+
x
1
e−t
t
dt
y
= −e−x
+
e−x
x
x2
y
+

x2
− x

y
+ (1 − x) y =

−x2
e−x
+ xe−x

+

x2
e−x
+ x2
x
1
e−t
t
dt − xe−x
− x
x
1
e−t
t
dt

+

x
x
1
e−t
t
dt − x2
x
1
e−t
t
dt

= 0
Chapter 1 in Review 35
30. . y = sin x
x
0
et2
cos t dt − cos x
x
0
et2
sin t dt
y
= sin x

ex2
cos x

+ cos x
x
0
et2
cos t dt − cos x

ex2
sin x

+ sin x
x
0
et2
sin t dt
= cos x
x
0
et2
cos t dt + sin x
x
0
et2
sin t dt
y
= cos x

ex2
cos x

− sin x
x
0
et2
cos t dt + sin x

ex2
sin x

+ cos x
x
0
et2
sin t dt
= ex2 
cos2
x + sin2
x

−
⎛
⎜
⎜
⎝
y
  
sin x
x
0
et2
cos t dt − cos x
x
0
et2
sin t dt
⎞
⎟
⎟
⎠
= ex2
− y
y
+ y = ex2
− y + y = ex2
31. Using implicit differentiation we get
x3
y3
= x3
+ 5
3x2
· y3
+ x3
· 3y2 dy
dx
= 3x2
3x2y3
3x2y2
+
x33y2
3x2y2
dy
dx
=
3x2
3x2y2
y + x
dy
dx
=
1
y2
32. Using implicit differentiation we get
(x − 7)2
+ y2
= 1
2(x − 7) + 2y
dy
dx
= 0
2y
dy
dx
= −2(x − 7)
y
dy
dx
= −(x − 7)

y
dy
dx
2
= (x − 7)2

dy
dx
2
=
(x − 7)2
y2
36 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS
Now from the original equation, isolating the first term leads to (x−7)2 = 1−y2. Continuing
from the last line of our proof we now have

dy
dx
2
=
(x − 7)2
y2
=
1 − y2
y2
=
1
y2
− 1
Adding 1 to both sides leads to the desired result.
33. Using implicit differentiation we get
y3
+ 3y = 2 − 3x
3y2
y
+ 3y
= −3
y2
y
+ y
= −1
(y2
+ 1)y
= −1
y
=
−1
y2 + 1
Differentiating the last line and remembering to use the quotient rule on the right side leads
to
y
=
2yy
(y2 + 1)2
Now since y = −1

y2 + 1

we can write the last equation as
y
=
2y
(y2 + 1)2
y
=
2y
(y2 + 1)2
−1
(y2 + 1)
= 2y

−1
y2 + 1
3
= 2y(y
)3
which is what we wanted to show.
34. Using implicit differentiation we get
y = e−xy
y
= e−xy
(−y − xy
)
ye−xy
+ xe−xy
y
+ y
= 0
ye−xy
+ (xe−xy
+ 1)y
= 0
Now since y = e−xy, substitute this into the last line to get
yy + (xy + 1)y
= 0
or (1 + xy)y + y2 = 0 which is what we wanted to show.
Chapter 1 in Review 37
35. Substituting y = c1 + cos3x, y = −3 sin 3x, and y = −9 cos 3x into the left-hand side of the
differential equation gives
y
+ 9y = −9 cos 3x + 9 (c1 + cos 3x) = −9 cos 3x + 9c1 + 9 cos 3x = 9c1.
Setting the result equal to the right-hand side of the differential equation yields 9c1 = 5 or
c1 =
5
9
. Thus y =
5
9
+ cos 3x.
36. Substituting y = c1 + c2x and y = c2 into the left-hand side of the differential equation gives
y
+ 2y = c2 + 2 (c1 + c2x) = c2 + 2c1 + 2c2x.
Setting the result equal to the right-hand side of the differential equation yields
c2 + 2c1 + 2c2x = 3x
(c2 + 2c1) + 2c2x = 0 + 3x.
Therefore, 2c2 = 3 or c2 =
3
2
, and
3
2
+ 2c1 = 0 or c1 = −
3
2
·
1
2
= −3
4. Thus y = −
3
4
+
3
2
x.
37. Differentiating and isolating c gives
y
= −ce−x
+ 4
ce−x
= 4 − y
c = ex

4 − y

.
Substituting the last equation into y = ce−x + 4x − 6 and simplifying yields
y = ex

4 − y

e−x
+ 4x − 6
y = 4 − y
+ 4x − 6
y
+ y = 4x − 2.
Alternatively, isolating c in y = ce−x + 4x − 6 gives
ce−x
= y − 4x + 6
c = ex
(y − 4x + 6) .
Differentiating the last equation yields
0 = ex
(y − 4x + 6) + ex

y
− 4

0 = ex

y
+ y − 4x + 2

0 = y
+ y − 4x + 2
y
+ y = 4x − 2.
38 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS
38. The slope
dy
dx
= 6 − 2x. Integrating with respect to x we have y = 6x − x2 + c1. Since the
graph must pass through the origin (0, 0), we have
y(0) = 6 (0) − (0)2
+ c1
0 = c1.
Thus y = 6x − x2.
In Problems 39–42, y = c1e−3x + c2ex + 4x is given as a two-parameter family of solutions of the
second-order differential equation y + 2y − 3y = −12x + 8.
39. If y(0) = 0 and y(0) = 0, then
c1 + c2 = 0
−3c1 + c2 = −4
Subtracting the second equation from the first gives us 4c1 = 4 or c1 = 1, and thus c2 = −1.
Therefore y = e−3x − ex + 4x.
40. If y(0) = 5 and y(0) = −11, then
c1 + c2 = 5
−3c1 + c2 = −15
Subtracting the second equation from the first gives us 4c1 = 20 or c1 = 5, and thus c2 = 0.
Therefore y = 5e−3x + 4x.
41. If y(1) = −2 and y(1) = 4, then
c1e−3
+ c2e = −6
−3c1e−3
+ c2e = 0
Subtracting the second equation from the first gives us 4c1 = −6 or c1 = −
3
2
e3, and thus
c2 = −
9
2
e−1. Therefore y = −
3
2
e−3x+3 −
9
2
ex−1 + 4x.
42. If y(−1) = 1 and y(−1) = 1, then
c1e3
+ c2e−1
= 5
−3c1e3
+ c2e−1
= −3
Subtracting the second equation from the first gives us 4c1 = 8 or c1 = 2e−3, and thus c2 = 3e.
Therefore y = 2e−3x−3 + 3ex+1 + 4x.
Chapter 1 in Review 39
43. We are to use the Leibniz’s rule
d
dx
v(x)
u(x)
F(x, t) dt = F(x, v(x))
dv
dx
− F(x, u(x))
du
dx
+
v(x)
u(x)
∂
∂x
F(x, t) dt
Since y(x) =
1
3
x
0
f(t) sin (3x − 3t) dt, take F(x, t) = f(t) sin (3x − 3t), u(x) = 0, and
v(x) = x to get
d
dx
y(x) =
d
dx
1
3
x
0
f(t) sin (3x − 3t) dt
=
1
3
F(x, x) · 1 − F(x, 0) · 0 +
x
0
∂
∂x
f(t) sin (3x − 3t) dt
=
1
3
f(x) sin (3x − 3x) +
x
0
3f(t) cos (3x − 3t) dt
=
x
0
f(t) cos (3x − 3t) dt
Apply the Leibniz’s rule a second time to y
=
x
0
f(t) cos (3x − 3t) dt by taking F(x, t) =
f(t) cos (3x − 3t), u(x) = 0, and v(x) = x to get
d
dx
y
(x) =
d
dx
x
0
f(t) cos (3x − 3t) dt
= F(x, x) · 1 − F(x, 0) · 0 +
x
0
∂
∂x
f(t) cos (3x − 3t) dt
= f(x) cos (3x − 3x) +
x
0
−3f(t) sin (3x − 3t) dt
= f(x) − 3
x
0
f(t) sin (3x − 3t) dt
Therefore y
= f(x)−3
x
0
f(t) sin (3x − 3t) dt. Now substituting y and y into the differential
equation we get
y
+ 9y = f(x) − 3
x
0
f(t) sin (3x − 3t) dt + 9 ·
1
3
x
0
f(t) sin (3x − 3t) dt = f(x)
40 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS
44. We are to use the Leibniz’s rule
d
dx
v(x)
u(x)
F(x, t) dt = F(x, v(x))
dv
dx
− F(x, u(x))
du
dx
+
v(x)
u(x)
∂
∂x
F(x, t) dt
Since y(x) =
π
0
ex cos t
dt, take F(x, t) = ex cos t, u(x) = 0, and v(x) = π to get
d
dx
y(x) =
d
dx
π
0
ex cos t
dt
= F(x, π) · 0 − F(x, 0) · 0 +
π
0
∂
∂x
ex cos t
dt
=
π
0
cos tex cos t
dt
Applying the Leibniz’s rule a second time to y
(x) =
π
0
cos tex cos t
dt by taking F(x, t) =
cos tex cos t, u(x) = 0, and v(x) = π to get
d
dx
y
(x) =
π
0
cos tex cos t
dt
= F(x, π) · 0 − F(x, 0) · 0 +
π
0
∂
∂x
cos tex cos t
dt
=
π
0
cos2
tex cos t
dt
An alternative form of y can be obtained by integrating by parts with respect to t:
y
= sin t · ex cos t



π
0
+
π
0
xecos t
sin2
t dt =
π
0
xex cos t
sin2
t dt = x
π
0
ex cos t
sin2
t dt.
We use this last form of y instead of the first in the differential equation:
xy
+ y
= x
π
0
ex cos t
cos2
t dt + x
π
0
ex cos t
sin2
t dt
= x
π
0
⎛
⎜
⎝
1
  
cos2
t + sin2
t
⎞
⎟
⎠ ex cos t
dt =
y
  
x
π
0
ex cos t
dt = xy
xy
+ y
= xy
Therefore, xy + y − xy = 0.
Chapter 1 in Review 41
45. From the graph we see that estimates for y0 and y1 are y0 = −3 and y1 = 0.
46. The differential equation is
dh
dt
= −
cA0
Aw

2gh .
Using A0 = π(1/24)2 = π/576, Aw = π(2)2 = 4π, and g = 32, this becomes
dh
dt
= −
cπ/576
4π
√
64h = −
c
288
√
h .
47. Let x(t) represent the height of the top of the rope at any time t with the positive direction
upward as indicated. The weight of the portion of the rope off the ground is given by
W = (x ft) · (1 lb/ft) = x. The mass of the rope is m = W/g = x/32. The net force is
F = 5 − W = 5 − x. Now by Newton’s second law we get
F =
d
dt
(mv) =
d
dt
 x
32
· v

= 5 − x
Now expand the derivative and remember that v = dx/dt to get
d
dt
 x
32
· v

= 5 − x
1
32

x
dv
dt
+ v
dx
dt

= 5 − x
x
d
dt

dx
dt

+

dx
dt

dx
dt
= 160 − 32x
x
d2x
dt2
+

dx
dt
2
+ 32x = 160
Ad

More Related Content

What's hot (20)

Aula 17: Separação da equação de Schrödinger em coordenadas cartesianas. Part...
Aula 17: Separação da equação de Schrödinger em coordenadas cartesianas. Part...Aula 17: Separação da equação de Schrödinger em coordenadas cartesianas. Part...
Aula 17: Separação da equação de Schrödinger em coordenadas cartesianas. Part...
Adriano Silva
 
Advanced Engineering Mathematics Solutions Manual.pdf
Advanced Engineering Mathematics Solutions Manual.pdfAdvanced Engineering Mathematics Solutions Manual.pdf
Advanced Engineering Mathematics Solutions Manual.pdf
Whitney Anderson
 
Nilai mutlak intoduction and definition
Nilai mutlak   intoduction and definitionNilai mutlak   intoduction and definition
Nilai mutlak intoduction and definition
AtikaFaradilla
 
System of Linear Equations
System of Linear EquationsSystem of Linear Equations
System of Linear Equations
Gautam Varshney
 
Higher Order Deriavatives
Higher Order DeriavativesHigher Order Deriavatives
Higher Order Deriavatives
Padme Amidala
 
Ecu tuning guide contents
Ecu tuning guide contentsEcu tuning guide contents
Ecu tuning guide contents
OBD bytes
 
07 periodic functions and fourier series
07 periodic functions and fourier series07 periodic functions and fourier series
07 periodic functions and fourier series
Krishna Gali
 
La Derivada
La DerivadaLa Derivada
La Derivada
Kike Prieto
 
Integrales 1
Integrales 1Integrales 1
Integrales 1
Juan Moreno Peña
 
Higher Differential Equation
Higher Differential Equation Higher Differential Equation
Higher Differential Equation
Abdul Hannan
 
LU FACTORIZATION METHOD
 LU FACTORIZATION METHOD LU FACTORIZATION METHOD
LU FACTORIZATION METHOD
reach2arkaELECTRICAL
 
IIT Jam math 2016 solutions BY Trajectoryeducation
IIT Jam math 2016 solutions BY TrajectoryeducationIIT Jam math 2016 solutions BY Trajectoryeducation
IIT Jam math 2016 solutions BY Trajectoryeducation
Dev Singh
 
solucionario de purcell 3
solucionario de purcell 3solucionario de purcell 3
solucionario de purcell 3
José Encalada
 
Matemática V - Facultad de Ingeniería Eléctrica y Electrónica de la Universid...
Matemática V - Facultad de Ingeniería Eléctrica y Electrónica de la Universid...Matemática V - Facultad de Ingeniería Eléctrica y Electrónica de la Universid...
Matemática V - Facultad de Ingeniería Eléctrica y Electrónica de la Universid...
Andy Juan Sarango Veliz
 
Partial differential equations
Partial differential equationsPartial differential equations
Partial differential equations
muhammadabullah
 
Exercice continuité et limites
Exercice continuité et limitesExercice continuité et limites
Exercice continuité et limites
Yessin Abdelhedi
 
application of first order ordinary Differential equations
application of first order ordinary Differential equationsapplication of first order ordinary Differential equations
application of first order ordinary Differential equations
Emdadul Haque Milon
 
The Epsilon-Delta Definition of a Limit
The Epsilon-Delta Definition of a LimitThe Epsilon-Delta Definition of a Limit
The Epsilon-Delta Definition of a Limit
ElizabethMastrangelo1
 
Bernoulli y ricatti
Bernoulli y ricattiBernoulli y ricatti
Bernoulli y ricatti
Tensor
 
Sistema de resorte y masa
Sistema de resorte y masaSistema de resorte y masa
Sistema de resorte y masa
Eduardo Fuentes
 
Aula 17: Separação da equação de Schrödinger em coordenadas cartesianas. Part...
Aula 17: Separação da equação de Schrödinger em coordenadas cartesianas. Part...Aula 17: Separação da equação de Schrödinger em coordenadas cartesianas. Part...
Aula 17: Separação da equação de Schrödinger em coordenadas cartesianas. Part...
Adriano Silva
 
Advanced Engineering Mathematics Solutions Manual.pdf
Advanced Engineering Mathematics Solutions Manual.pdfAdvanced Engineering Mathematics Solutions Manual.pdf
Advanced Engineering Mathematics Solutions Manual.pdf
Whitney Anderson
 
Nilai mutlak intoduction and definition
Nilai mutlak   intoduction and definitionNilai mutlak   intoduction and definition
Nilai mutlak intoduction and definition
AtikaFaradilla
 
System of Linear Equations
System of Linear EquationsSystem of Linear Equations
System of Linear Equations
Gautam Varshney
 
Higher Order Deriavatives
Higher Order DeriavativesHigher Order Deriavatives
Higher Order Deriavatives
Padme Amidala
 
Ecu tuning guide contents
Ecu tuning guide contentsEcu tuning guide contents
Ecu tuning guide contents
OBD bytes
 
07 periodic functions and fourier series
07 periodic functions and fourier series07 periodic functions and fourier series
07 periodic functions and fourier series
Krishna Gali
 
Higher Differential Equation
Higher Differential Equation Higher Differential Equation
Higher Differential Equation
Abdul Hannan
 
IIT Jam math 2016 solutions BY Trajectoryeducation
IIT Jam math 2016 solutions BY TrajectoryeducationIIT Jam math 2016 solutions BY Trajectoryeducation
IIT Jam math 2016 solutions BY Trajectoryeducation
Dev Singh
 
solucionario de purcell 3
solucionario de purcell 3solucionario de purcell 3
solucionario de purcell 3
José Encalada
 
Matemática V - Facultad de Ingeniería Eléctrica y Electrónica de la Universid...
Matemática V - Facultad de Ingeniería Eléctrica y Electrónica de la Universid...Matemática V - Facultad de Ingeniería Eléctrica y Electrónica de la Universid...
Matemática V - Facultad de Ingeniería Eléctrica y Electrónica de la Universid...
Andy Juan Sarango Veliz
 
Partial differential equations
Partial differential equationsPartial differential equations
Partial differential equations
muhammadabullah
 
Exercice continuité et limites
Exercice continuité et limitesExercice continuité et limites
Exercice continuité et limites
Yessin Abdelhedi
 
application of first order ordinary Differential equations
application of first order ordinary Differential equationsapplication of first order ordinary Differential equations
application of first order ordinary Differential equations
Emdadul Haque Milon
 
The Epsilon-Delta Definition of a Limit
The Epsilon-Delta Definition of a LimitThe Epsilon-Delta Definition of a Limit
The Epsilon-Delta Definition of a Limit
ElizabethMastrangelo1
 
Bernoulli y ricatti
Bernoulli y ricattiBernoulli y ricatti
Bernoulli y ricatti
Tensor
 
Sistema de resorte y masa
Sistema de resorte y masaSistema de resorte y masa
Sistema de resorte y masa
Eduardo Fuentes
 

Similar to Answers to Problems for Advanced Engineering Mathematics, 7th Edition by Dennis Zill (20)

Solutions for Problems in "A First Course in Differential Equations" (11th Ed...
Solutions for Problems in "A First Course in Differential Equations" (11th Ed...Solutions for Problems in "A First Course in Differential Equations" (11th Ed...
Solutions for Problems in "A First Course in Differential Equations" (11th Ed...
electricaleng2024
 
Answers to Problems for Advanced Engineering Mathematics, 6th Edition – Denni...
Answers to Problems for Advanced Engineering Mathematics, 6th Edition – Denni...Answers to Problems for Advanced Engineering Mathematics, 6th Edition – Denni...
Answers to Problems for Advanced Engineering Mathematics, 6th Edition – Denni...
financemanag
 
Sect1 1
Sect1 1Sect1 1
Sect1 1
inKFUPM
 
differentiol equation.pptx
differentiol equation.pptxdifferentiol equation.pptx
differentiol equation.pptx
PlanningHCEGC
 
MRS EMMAH.pdf
MRS EMMAH.pdfMRS EMMAH.pdf
MRS EMMAH.pdf
Kasungwa
 
First order differential equations
First order differential equationsFirst order differential equations
First order differential equations
Dr.Summiya Parveen
 
Documents.mx eduv
Documents.mx eduvDocuments.mx eduv
Documents.mx eduv
Osmar Meraz
 
Sect5 1
Sect5 1Sect5 1
Sect5 1
inKFUPM
 
258 lecnot2
258 lecnot2258 lecnot2
258 lecnot2
KELLY SANTOS
 
Ma 104 differential equations
Ma 104 differential equationsMa 104 differential equations
Ma 104 differential equations
arvindpt1
 
Engr 213 midterm 2b sol 2009
Engr 213 midterm 2b sol 2009Engr 213 midterm 2b sol 2009
Engr 213 midterm 2b sol 2009
akabaka12
 
Ordinary Differential Equations: Variable separation method
Ordinary Differential Equations: Variable separation method  Ordinary Differential Equations: Variable separation method
Ordinary Differential Equations: Variable separation method
AMINULISLAM439
 
MFMTQP_MAT_nda question paper for nda class 12
MFMTQP_MAT_nda question paper  for nda class 12MFMTQP_MAT_nda question paper  for nda class 12
MFMTQP_MAT_nda question paper for nda class 12
EMALLIKARJUNAREDDY
 
Calculus 05-6 ccccccccccccccccccccccc.ppsx
Calculus 05-6 ccccccccccccccccccccccc.ppsxCalculus 05-6 ccccccccccccccccccccccc.ppsx
Calculus 05-6 ccccccccccccccccccccccc.ppsx
MohamedAli899919
 
2500-Solved-Problems-in-Differential-Equations.pdf
2500-Solved-Problems-in-Differential-Equations.pdf2500-Solved-Problems-in-Differential-Equations.pdf
2500-Solved-Problems-in-Differential-Equations.pdf
LoreimayAndaya1
 
Chapter 5 assignment
Chapter 5 assignmentChapter 5 assignment
Chapter 5 assignment
KarunaGupta1982
 
Engr 213 midterm 2a sol 2009
Engr 213 midterm 2a sol 2009Engr 213 midterm 2a sol 2009
Engr 213 midterm 2a sol 2009
akabaka12
 
Persamaan Differensial Biasa 2014
Persamaan Differensial Biasa 2014 Persamaan Differensial Biasa 2014
Persamaan Differensial Biasa 2014
Rani Sulvianuri
 
Ordinary differential equations
Ordinary differential equationsOrdinary differential equations
Ordinary differential equations
Ahmed Haider
 
Elementary differential equations with boundary value problems solutions
Elementary differential equations with boundary value problems solutionsElementary differential equations with boundary value problems solutions
Elementary differential equations with boundary value problems solutions
Hon Wa Wong
 
Solutions for Problems in "A First Course in Differential Equations" (11th Ed...
Solutions for Problems in "A First Course in Differential Equations" (11th Ed...Solutions for Problems in "A First Course in Differential Equations" (11th Ed...
Solutions for Problems in "A First Course in Differential Equations" (11th Ed...
electricaleng2024
 
Answers to Problems for Advanced Engineering Mathematics, 6th Edition – Denni...
Answers to Problems for Advanced Engineering Mathematics, 6th Edition – Denni...Answers to Problems for Advanced Engineering Mathematics, 6th Edition – Denni...
Answers to Problems for Advanced Engineering Mathematics, 6th Edition – Denni...
financemanag
 
differentiol equation.pptx
differentiol equation.pptxdifferentiol equation.pptx
differentiol equation.pptx
PlanningHCEGC
 
MRS EMMAH.pdf
MRS EMMAH.pdfMRS EMMAH.pdf
MRS EMMAH.pdf
Kasungwa
 
First order differential equations
First order differential equationsFirst order differential equations
First order differential equations
Dr.Summiya Parveen
 
Documents.mx eduv
Documents.mx eduvDocuments.mx eduv
Documents.mx eduv
Osmar Meraz
 
Ma 104 differential equations
Ma 104 differential equationsMa 104 differential equations
Ma 104 differential equations
arvindpt1
 
Engr 213 midterm 2b sol 2009
Engr 213 midterm 2b sol 2009Engr 213 midterm 2b sol 2009
Engr 213 midterm 2b sol 2009
akabaka12
 
Ordinary Differential Equations: Variable separation method
Ordinary Differential Equations: Variable separation method  Ordinary Differential Equations: Variable separation method
Ordinary Differential Equations: Variable separation method
AMINULISLAM439
 
MFMTQP_MAT_nda question paper for nda class 12
MFMTQP_MAT_nda question paper  for nda class 12MFMTQP_MAT_nda question paper  for nda class 12
MFMTQP_MAT_nda question paper for nda class 12
EMALLIKARJUNAREDDY
 
Calculus 05-6 ccccccccccccccccccccccc.ppsx
Calculus 05-6 ccccccccccccccccccccccc.ppsxCalculus 05-6 ccccccccccccccccccccccc.ppsx
Calculus 05-6 ccccccccccccccccccccccc.ppsx
MohamedAli899919
 
2500-Solved-Problems-in-Differential-Equations.pdf
2500-Solved-Problems-in-Differential-Equations.pdf2500-Solved-Problems-in-Differential-Equations.pdf
2500-Solved-Problems-in-Differential-Equations.pdf
LoreimayAndaya1
 
Engr 213 midterm 2a sol 2009
Engr 213 midterm 2a sol 2009Engr 213 midterm 2a sol 2009
Engr 213 midterm 2a sol 2009
akabaka12
 
Persamaan Differensial Biasa 2014
Persamaan Differensial Biasa 2014 Persamaan Differensial Biasa 2014
Persamaan Differensial Biasa 2014
Rani Sulvianuri
 
Ordinary differential equations
Ordinary differential equationsOrdinary differential equations
Ordinary differential equations
Ahmed Haider
 
Elementary differential equations with boundary value problems solutions
Elementary differential equations with boundary value problems solutionsElementary differential equations with boundary value problems solutions
Elementary differential equations with boundary value problems solutions
Hon Wa Wong
 
Ad

Recently uploaded (20)

Peer Assessment_ Unit 2 Skills Development for Live Performance - for Libby.docx
Peer Assessment_ Unit 2 Skills Development for Live Performance - for Libby.docxPeer Assessment_ Unit 2 Skills Development for Live Performance - for Libby.docx
Peer Assessment_ Unit 2 Skills Development for Live Performance - for Libby.docx
19lburrell
 
MICROBIAL GENETICS -tranformation and tranduction.pdf
MICROBIAL GENETICS -tranformation and tranduction.pdfMICROBIAL GENETICS -tranformation and tranduction.pdf
MICROBIAL GENETICS -tranformation and tranduction.pdf
DHARMENDRA SAHU
 
GENERAL QUIZ PRELIMS | QUIZ CLUB OF PSGCAS | 4 MARCH 2025 .pdf
GENERAL QUIZ PRELIMS | QUIZ CLUB OF PSGCAS | 4 MARCH 2025 .pdfGENERAL QUIZ PRELIMS | QUIZ CLUB OF PSGCAS | 4 MARCH 2025 .pdf
GENERAL QUIZ PRELIMS | QUIZ CLUB OF PSGCAS | 4 MARCH 2025 .pdf
Quiz Club of PSG College of Arts & Science
 
The role of wall art in interior designing
The role of wall art in interior designingThe role of wall art in interior designing
The role of wall art in interior designing
meghaark2110
 
Classification of mental disorder in 5th semester bsc. nursing and also used ...
Classification of mental disorder in 5th semester bsc. nursing and also used ...Classification of mental disorder in 5th semester bsc. nursing and also used ...
Classification of mental disorder in 5th semester bsc. nursing and also used ...
parmarjuli1412
 
Dastur_ul_Amal under Jahangir Key Features.pptx
Dastur_ul_Amal under Jahangir Key Features.pptxDastur_ul_Amal under Jahangir Key Features.pptx
Dastur_ul_Amal under Jahangir Key Features.pptx
omorfaruqkazi
 
YSPH VMOC Special Report - Measles Outbreak Southwest US 5-14-2025 .pptx
YSPH VMOC Special Report - Measles Outbreak  Southwest US 5-14-2025  .pptxYSPH VMOC Special Report - Measles Outbreak  Southwest US 5-14-2025  .pptx
YSPH VMOC Special Report - Measles Outbreak Southwest US 5-14-2025 .pptx
Yale School of Public Health - The Virtual Medical Operations Center (VMOC)
 
How to Manage Cross Selling in Odoo 18 Sales
How to Manage Cross Selling in Odoo 18 SalesHow to Manage Cross Selling in Odoo 18 Sales
How to Manage Cross Selling in Odoo 18 Sales
Celine George
 
MCQ PHYSIOLOGY II (DR. NASIR MUSTAFA) MCQS)
MCQ PHYSIOLOGY II (DR. NASIR MUSTAFA) MCQS)MCQ PHYSIOLOGY II (DR. NASIR MUSTAFA) MCQS)
MCQ PHYSIOLOGY II (DR. NASIR MUSTAFA) MCQS)
Dr. Nasir Mustafa
 
"Bridging Cultures Through Holiday Cards: 39 Students Celebrate Global Tradit...
"Bridging Cultures Through Holiday Cards: 39 Students Celebrate Global Tradit..."Bridging Cultures Through Holiday Cards: 39 Students Celebrate Global Tradit...
"Bridging Cultures Through Holiday Cards: 39 Students Celebrate Global Tradit...
AlionaBujoreanu
 
materi 3D Augmented Reality dengan assemblr
materi 3D Augmented Reality dengan assemblrmateri 3D Augmented Reality dengan assemblr
materi 3D Augmented Reality dengan assemblr
fatikhatunnajikhah1
 
Peer Assesment- Libby.docx..............
Peer Assesment- Libby.docx..............Peer Assesment- Libby.docx..............
Peer Assesment- Libby.docx..............
19lburrell
 
How to Configure Extra Steps During Checkout in Odoo 18 Website
How to Configure Extra Steps During Checkout in Odoo 18 WebsiteHow to Configure Extra Steps During Checkout in Odoo 18 Website
How to Configure Extra Steps During Checkout in Odoo 18 Website
Celine George
 
Final Evaluation.docx...........................
Final Evaluation.docx...........................Final Evaluation.docx...........................
Final Evaluation.docx...........................
l1bbyburrell
 
Letter to Secretary Linda McMahon from U.S. Senators
Letter to Secretary Linda McMahon from U.S. SenatorsLetter to Secretary Linda McMahon from U.S. Senators
Letter to Secretary Linda McMahon from U.S. Senators
Mebane Rash
 
Launch of The State of Global Teenage Career Preparation - Andreas Schleicher...
Launch of The State of Global Teenage Career Preparation - Andreas Schleicher...Launch of The State of Global Teenage Career Preparation - Andreas Schleicher...
Launch of The State of Global Teenage Career Preparation - Andreas Schleicher...
EduSkills OECD
 
LDMMIA 2024 Crystal Gold Lecture 1 Bonus
LDMMIA 2024 Crystal Gold Lecture 1 BonusLDMMIA 2024 Crystal Gold Lecture 1 Bonus
LDMMIA 2024 Crystal Gold Lecture 1 Bonus
LDM & Mia eStudios
 
How to Change Sequence Number in Odoo 18 Sale Order
How to Change Sequence Number in Odoo 18 Sale OrderHow to Change Sequence Number in Odoo 18 Sale Order
How to Change Sequence Number in Odoo 18 Sale Order
Celine George
 
A report on the county distress rankings in NC
A report on the county distress rankings in NCA report on the county distress rankings in NC
A report on the county distress rankings in NC
Mebane Rash
 
Aerospace Engineering Homework Help Guide – Expert Support for Academic Success
Aerospace Engineering Homework Help Guide – Expert Support for Academic SuccessAerospace Engineering Homework Help Guide – Expert Support for Academic Success
Aerospace Engineering Homework Help Guide – Expert Support for Academic Success
online college homework help
 
Peer Assessment_ Unit 2 Skills Development for Live Performance - for Libby.docx
Peer Assessment_ Unit 2 Skills Development for Live Performance - for Libby.docxPeer Assessment_ Unit 2 Skills Development for Live Performance - for Libby.docx
Peer Assessment_ Unit 2 Skills Development for Live Performance - for Libby.docx
19lburrell
 
MICROBIAL GENETICS -tranformation and tranduction.pdf
MICROBIAL GENETICS -tranformation and tranduction.pdfMICROBIAL GENETICS -tranformation and tranduction.pdf
MICROBIAL GENETICS -tranformation and tranduction.pdf
DHARMENDRA SAHU
 
The role of wall art in interior designing
The role of wall art in interior designingThe role of wall art in interior designing
The role of wall art in interior designing
meghaark2110
 
Classification of mental disorder in 5th semester bsc. nursing and also used ...
Classification of mental disorder in 5th semester bsc. nursing and also used ...Classification of mental disorder in 5th semester bsc. nursing and also used ...
Classification of mental disorder in 5th semester bsc. nursing and also used ...
parmarjuli1412
 
Dastur_ul_Amal under Jahangir Key Features.pptx
Dastur_ul_Amal under Jahangir Key Features.pptxDastur_ul_Amal under Jahangir Key Features.pptx
Dastur_ul_Amal under Jahangir Key Features.pptx
omorfaruqkazi
 
How to Manage Cross Selling in Odoo 18 Sales
How to Manage Cross Selling in Odoo 18 SalesHow to Manage Cross Selling in Odoo 18 Sales
How to Manage Cross Selling in Odoo 18 Sales
Celine George
 
MCQ PHYSIOLOGY II (DR. NASIR MUSTAFA) MCQS)
MCQ PHYSIOLOGY II (DR. NASIR MUSTAFA) MCQS)MCQ PHYSIOLOGY II (DR. NASIR MUSTAFA) MCQS)
MCQ PHYSIOLOGY II (DR. NASIR MUSTAFA) MCQS)
Dr. Nasir Mustafa
 
"Bridging Cultures Through Holiday Cards: 39 Students Celebrate Global Tradit...
"Bridging Cultures Through Holiday Cards: 39 Students Celebrate Global Tradit..."Bridging Cultures Through Holiday Cards: 39 Students Celebrate Global Tradit...
"Bridging Cultures Through Holiday Cards: 39 Students Celebrate Global Tradit...
AlionaBujoreanu
 
materi 3D Augmented Reality dengan assemblr
materi 3D Augmented Reality dengan assemblrmateri 3D Augmented Reality dengan assemblr
materi 3D Augmented Reality dengan assemblr
fatikhatunnajikhah1
 
Peer Assesment- Libby.docx..............
Peer Assesment- Libby.docx..............Peer Assesment- Libby.docx..............
Peer Assesment- Libby.docx..............
19lburrell
 
How to Configure Extra Steps During Checkout in Odoo 18 Website
How to Configure Extra Steps During Checkout in Odoo 18 WebsiteHow to Configure Extra Steps During Checkout in Odoo 18 Website
How to Configure Extra Steps During Checkout in Odoo 18 Website
Celine George
 
Final Evaluation.docx...........................
Final Evaluation.docx...........................Final Evaluation.docx...........................
Final Evaluation.docx...........................
l1bbyburrell
 
Letter to Secretary Linda McMahon from U.S. Senators
Letter to Secretary Linda McMahon from U.S. SenatorsLetter to Secretary Linda McMahon from U.S. Senators
Letter to Secretary Linda McMahon from U.S. Senators
Mebane Rash
 
Launch of The State of Global Teenage Career Preparation - Andreas Schleicher...
Launch of The State of Global Teenage Career Preparation - Andreas Schleicher...Launch of The State of Global Teenage Career Preparation - Andreas Schleicher...
Launch of The State of Global Teenage Career Preparation - Andreas Schleicher...
EduSkills OECD
 
LDMMIA 2024 Crystal Gold Lecture 1 Bonus
LDMMIA 2024 Crystal Gold Lecture 1 BonusLDMMIA 2024 Crystal Gold Lecture 1 Bonus
LDMMIA 2024 Crystal Gold Lecture 1 Bonus
LDM & Mia eStudios
 
How to Change Sequence Number in Odoo 18 Sale Order
How to Change Sequence Number in Odoo 18 Sale OrderHow to Change Sequence Number in Odoo 18 Sale Order
How to Change Sequence Number in Odoo 18 Sale Order
Celine George
 
A report on the county distress rankings in NC
A report on the county distress rankings in NCA report on the county distress rankings in NC
A report on the county distress rankings in NC
Mebane Rash
 
Aerospace Engineering Homework Help Guide – Expert Support for Academic Success
Aerospace Engineering Homework Help Guide – Expert Support for Academic SuccessAerospace Engineering Homework Help Guide – Expert Support for Academic Success
Aerospace Engineering Homework Help Guide – Expert Support for Academic Success
online college homework help
 
Ad

Answers to Problems for Advanced Engineering Mathematics, 7th Edition by Dennis Zill

  • 1. Chapter 1 Introduction to Differential Equations 1.1 Definitions and Terminology 1. Second order; linear 2. Third order; nonlinear because of (dy/dx)4 3. Fourth order; linear 4. Second order; nonlinear because of cos(r + u) 5. Second order; nonlinear because of (dy/dx)2 or 1 + (dy/dx)2 6. Second order; nonlinear because of R2 7. Third order; linear 8. Second order; nonlinear because of ẋ2 9. First order; nonlinear because of sin dy dx 10. First order; linear in x. 11. Writing the differential equation in the form x(dy/dx) + y2 = 1, we see that it is nonlinear in y because of y2. However, writing it in the form (y2 − 1)(dx/dy) + x = 0, we see that it is linear in x. 12. Writing the differential equation in the form u(dv/du) + (1 + u)v = ueu, we see that it is linear in v. However, writing it in the form (v + uv − ueu)(du/dv) + u = 0, we see that it is nonlinear in u. 13. From y = e−x/2 we obtain y = −1 2e−x/2. Then 2y + y = −e−x/2 + e−x/2 = 0. 14. From y = 6 5 − 6 5 e−20t we obtain dy/dt = 24e−20t, so that dy dt + 20y = 24e−20t + 20 6 5 − 6 5 e−20t = 24. 1 Contact me in order to access the whole complete document. WhatsApp: https://wa.me/message/2H3BV2L5TTSUF1 Email: smtb98@gmail.com Telegram: https://t.me/solutionmanual s m t b 9 8 @ g m a i l . c o m s m t b 9 8 @ g m a i l . c o m complete document is available on https://unihelp.xyz/ *** contact me if site not loaded
  • 2. 2 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS 15. From y = e3x cos 2x we obtain y = 3e3x cos 2x−2e3x sin 2x and y = 5e3x cos 2x−12e3x sin 2x, so that y − 6y + 13y = 0. 16. From y = − cos x ln(sec x + tan x) we obtain y = −1 + sin x ln(sec x + tan x) and y = tan x + cos x ln(sec x + tan x). Then y + y = tan x. 17. The domain of the function, found by solving x+2 ≥ 0, is [−2, ∞). From y = 1+2(x+2)−1/2 we have (y − x)y = (y − x)[1 + (2(x + 2)−1/2 ] = y − x + 2(y − x)(x + 2)−1/2 = y − x + 2[x + 4(x + 2)1/2 − x](x + 2)−1/2 = y − x + 8(x + 2)1/2 (x + 2)−1/2 = y − x + 8. An interval of definition for the solution of the differential equation is (−2, ∞) because y is not defined at x = −2. 18. Since tan x is not defined for x = π/2 + nπ, n an integer, the domain of y = 5 tan 5x is {x 5x = π/2 + nπ} or {x x = π/10 + nπ/5}. From y = 25 sec2 5x we have y = 25(1 + tan2 5x) = 25 + 25 tan2 5x = 25 + y2 . An interval of definition for the solution of the differential equation is (−π/10, π/10). Another interval is (π/10, 3π/10), and so on. 19. The domain of the function is {x 4 − x2 = 0} or {x x = −2 or x = 2}. From y = 2x/(4 − x2)2 we have y = 2x 1 4 − x2 2 = 2xy2 . An interval of definition for the solution of the differential equation is (−2, 2). Other intervals are (−∞, −2) and (2, ∞). 20. The function is y = 1/ √ 1 − sin x , whose domain is obtained from 1 − sin x = 0 or sin x = 1. Thus, the domain is {x x = π/2 + 2nπ}. From y = −1 2(1 − sin x)−3/2(− cos x) we have 2y = (1 − sin x)−3/2 cos x = [(1 − sin x)−1/2 ]3 cos x = y3 cos x. An interval of definition for the solution of the differential equation is (π/2, 5π/2). Another one is (5π/2, 9π/2), and so on.
  • 3. 1.1 Definitions and Terminology 3 21. Writing ln(2X − 1) − ln(X − 1) = t and differentiating implicitly we obtain 2 2X − 1 dX dt − 1 X − 1 dX dt = 1 2 2X − 1 − 1 X − 1 dX dt = 1 2X − 2 − 2X + 1 (2X − 1) (X − 1) dX dt = 1 dX dt = −(2X − 1)(X − 1) = (X − 1)(1 − 2X). – 4 –2 2 4 t – 4 –2 2 4 x Exponentiating both sides of the implicit solution we obtain 2X − 1 X − 1 = et 2X − 1 = Xet − et (et − 1) = (et − 2)X X = et − 1 et − 2 . Solving et − 2 = 0 we get t = ln 2. Thus the solution is defined on (−∞, ln 2) or on (ln 2, ∞). The graph of the solution defined on (−∞, ln 2) is dashed, and the graph of the solution defined on (ln 2, ∞) is solid. 22. Implicitly differentiating the solution, we obtain −2x2 dy dx − 4xy + 2y dy dx = 0 −x2 dy − 2xy dx + y dy = 0 2xy dx + (x2 − y)dy = 0. Using the quadratic formula to solve y2 − 2x2y − 1 = 0 for y, we get y = 2x2 ± √ 4x4 + 4 /2 = x2 ± √ x4 + 1 . Thus, two explicit solutions are y1 = x2 + √ x4 + 1 and y2 = x2 − √ x4 + 1 . Both solutions are defined on (−∞, ∞). The graph of y1(x) is solid and the graph of y2 is dashed. – 4 –2 2 4 x – 4 –2 2 4 y
  • 4. 4 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS 23. Differentiating P = c1et/ 1 + c1et we obtain dP dt = 1 + c1et c1et − c1et · c1et (1 + c1et)2 = c1et 1 + c1et 1 + c1et − c1et 1 + c1et = c1et 1 + c1et 1 − c1et 1 + c1et = P(1 − P). 24. Differentiating y = 2x2 − 1 + c1e−2x2 we obtain dy dx = 4x − 4xc1e−2x2 , so that dy dx + 4xy = 4x − 4xc1e−2x2 + 8x3 − 4x + 4c1xe−x2 = 8x3 25. From y = c1e2x + c2xe2x we obtain dy dx = (2c1 + c2)e2x + 2c2xe2x and d2y dx2 = (4c1 + 4c2)e2x + 4c2xe2x , so that d2y dx2 − 4 dy dx + 4y = (4c1 + 4c2 − 8c1 − 4c2 + 4c1)e2x + (4c2 − 8c2 + 4c2)xe2x = 0. 26. From y = c1x−1 + c2x + c3x ln x + 4x2 we obtain dy dx = −c1x−2 + c2 + c3 + c3 ln x + 8x, d2y dx2 = 2c1x−3 + c3x−1 + 8, and d3y dx3 = −6c1x−4 − c3x−2 , so that x3 d3y dx3 + 2x2 d2y dx2 − x dy dx + y = (−6c1 + 4c1 + c1 + c1)x−1 + (−c3 + 2c3 − c2 − c3 + c2)x + (−c3 + c3)x ln x + (16 − 8 + 4)x2 = 12x2 In Problems 27–30, we use the Product Rule and the derivative of an integral ((12) of this section): d dx x a g(t) dt = g(x). 27. Differentiating y = e3x x 1 e−3t t dt we obtain dy dx = e3x x 1 e−3t t dt + e−3t x · e3x or dy dx = e3x x 1 e−3t t dt + 1 x , so that x dy dx − 3xy = x e3x x 1 e−3t t dt + 1 x − 3x e3x x 1 e−3t t dt = xe3x x 1 e−3t t dt + 1 − 3xe3x x 1 e−3t t dt = 1
  • 5. 1.1 Definitions and Terminology 5 28. Differentiating y = √ x x 4 cos t √ t dt we obtain dy dx = 1 2 √ x x 4 cos t √ t dt + cos x √ x · √ x or dy dx = 1 2 √ x x 4 cos t √ t dt + cos x, so that 2x dy dx − y = 2x 1 2 √ x x 4 cos t √ t dt + cos x − √ x x 4 cos t √ t dt = √ x x 4 cos t √ t dt + 2x cos x − √ x x 4 cos t √ t dt = 2x cos x 29. Differentiating y = 5 x + 10 x x 1 sin t t dt we obtain dy dx = − 5 x2 − 10 x2 x 1 sin t t dt + sin x x · 10 x or dy dx = − 5 x2 − 10 x2 x 1 sin t t dt + 10 sin x x2 , so that x2 dy dx + xy = x2 − 5 x2 − 10 x2 x 1 sin t t dt + 10 sin x x2 + x 5 x + 10 x x 1 sin t t dt = −5 − 10 x 1 sin t t dt + 10 sin x + 5 + 10 x 1 sin t t dt = 10 sin x 30. Differentiating y = e−x2 +e−x2 x 0 et2 dt we obtain dy dx = −2xe−x2 −2xe−x2 x 0 et2 dt+ex2 ·e−x2 or dy dx = −2xe−x2 − 2xe−x2 x 0 et2 dt + 1, so that dy dx + 2xy = −2xe−x2 − 2xe−x2 x 0 et2 dt + 1 + 2x e−x2 + e−x2 x 0 et2 dt = −2xe−x2 − 2xe−x2 x 0 et2 dt + 1 + 2xe−x2 + 2xe−x2 x 0 et2 dt = 1 31. From y = −x2, x 0 x2, x ≥ 0 we obtain y = −2x, x 0 2x, x ≥ 0 so that xy − 2y = 0. 32. The function y(x) is not continuous at x = 0 since lim x→0− y(x) = 5 and lim x→0+ y(x) = −5. Thus, y(x) does not exist at x = 0. 33. Force the function y = emx into the equation y + 2y = 0 to get (emx ) + 2(emx ) = 0 memx + 2emx = 0 emx (m + 2) = 0 Now since emx 0 for all values of x, we must have m = −2 and so y = e−2x is a solution.
  • 6. 6 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS 34. Force the function y = emx into the equation 3y − 4y = 0 to get 3(emx ) − 4(emx ) = 0 3memx − 4emx = 0 emx (3m − 4) = 0 Now since emx 0 for all values of x, we must have m = 4/3 and so y = e4x/3 is a solution. 35. Force the function y = emx into the equation y − 5y + 6y = 0 to get (emx ) − 5(emx ) + 6(emx ) = 0 m2 emx − 5memx + 6emx = 0 emx (m2 − 5m + 6) = 0 emx (m − 2)(m − 3) = 0 Now since emx 0 for all values of x, we must have m = 2 and m = 3 therefore y = e2x and y = e3x are solutions. 36. Force the function y = emx into the equation 2y + 9y − 5y = 0 to get 2(emx ) + 9(emx ) − 5(emx ) = 0 2m2 emx + 9memx − 5emx = 0 emx (2m2 + 9m − 5) = 0 emx (m + 5)(2m − 1) = 0 Now since emx 0 for all values of x , we must have m = −5 and m = 1/2 therefore y = e−5x and y = ex/2 are solutions. 37. Force the function y = xm into the equation xy + 2y = 0 to get x · (xm ) + 2(xm ) = 0 x · m(m − 1)xm−2 + 2mxm−1 = 0 (m2 − m)xm−1 + 2mxm−1 = 0 xm−1 [m2 + m] = 0 xm−1 [m(m + 1)] = 0 The last line implies that m = 0 and m = −1 therefore y = x0 = 1 and y = x−1 are solutions.
  • 7. 1.1 Definitions and Terminology 7 38. Force the function y = xm into the equation 4x2y + y = 0 to get 4x2 · (xm ) + (xm ) = 0 4x2 · m(m − 1)xm−2 + xm = 0 4(m2 − m)xm + xm = 0 xm [4m2 − 4m + 1] = 0 xm [(2m − 1)(2m − 1)] = 0 The last line implies m = 1/2 therefore y = x1/2 is a solution. 39. Force the function y = xm into the equation x2y − 7xy + 15y = 0 to get x2 · (xm ) − 7x · (xm ) + 15(xm ) = 0 x2 · m(m − 1)xm−2 − 7x · mxm−1 + 15xm = 0 (m2 − m)xm − 7mxm + 15xm = 0 xm [m2 − 8m + 15] = 0 xm [(m − 3)(m − 5)] = 0 The last line implies that m = 3 and m = 5 therefore y = x3 and y = x5 are solutions. 40. Force the function y = xm into the equation x2y − 3xy + 3y = 0 to get x2 · (xm ) − 3x · (xm ) + 3(xm ) = 0 x2 · m(m − 1)(m − 2)xm−3 − 3x · m(m − 1)xm−2 + 3 · mxm−1 = 0 (m3 − 3m2 + 2m)xm−1 − 3(m2 − m)xm−1 + 3mxm−1 = 0 xm−1 [m3 − 6m2 + 8m] = 0 xm [m(m − 2)(m − 4)] = 0 The last line implies that m = 0, m = 2 and m = 4 therefore y = x0 = 1, y = x2, and y = x4 are solutions. In Problems 41–44, we substitute y = c into the differential equations and use y = 0 and y = 0. 41. Solving 5c = 10 we see that y = 2 is a constant solution. 42. Solving c2 + 2c − 3 = (c + 3)(c − 1) = 0 we see that y = −3 and y = 1 are constant solutions. 43. Since 1/(c − 1) = 0 has no solutions, the differential equation has no constant solutions.
  • 8. 8 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS 44. Solving 6c = 10 we see that y = 5/3 is a constant solution. 45. From y = (x + c1)2 we obtain dy dx 2 = (2 (x + c1))2 = 4 (x + c1)2 . Then 4y = 4(x + c1)2 . Inspection of the differential equation reveals that y = 0 is a solution of the differential equation but is not a member of the one-parameter family y = (x + c1)2. 46. From y = 3 sin (x + c1) we obtain dy dx 2 = (3 cos (x + c1))2 = 9 cos2 (x + c1). Then 9 − y2 = 9 − (3 sin (x + c1))2 = 9 − 9 sin2 (x + c1) = 9 1 − sin2 (x + c1) = 9 cos2 (x + c1). Inspection of the differential equation reveals that y = 3 and y = −3 are solutions of the differential equation but not members of the one-parameter family y = 3 sin (x + c1). 47. From x − 16 − y2 = c1 we obtain d dx x − 16 − y2 = d dx (c1) 1 − 1 2 16 − y2 −1/2 −2y dy dx = 0 1 + y 16 − y2 dy dx = 0 y 16 − y2 dy dx = −1 y dy dx = − 16 − y2. Then y dy dx + 16 − y2 = − 16 − y2 + 16 − y2 = 0. Inspection of the differential equation reveals that y = −4 and y = 4 are solutions of the given differential equation but are not members of the one-parameter family x − 16 − y2 = c1.
  • 9. 1.1 Definitions and Terminology 9 48. From y = x − (x − c1)2 we obtain dy dx 2 = (1 − 2 (x − c1))2 = 1 − 4 (x − c1) + 4 (x − c1)2 = 1 − 4x + 4c1 + 4x2 − 8xc1 + 4c2 1 −2 dy dx = −2 (1 − 2 (x − c1)) = −2 + 4x − 4c1 4y = 4 x − (x − c1)2 = 4 x − x2 + 2xc1 − c2 1 = 4x − 4x2 + 8xc1 − 4c2 1. Then dy dx 2 − 2 dy dx + 4y = 4x − 1. Inspection of the differential equation reveals that y = x is a solution of the differential equation but not a member of the one-parameter family y = x − (x − c1)2. 49. From x = e−2t + 3e6t and y = −e−2t + 5e6t we obtain dx dt = −2e−2t + 18e6t and dy dt = 2e−2t + 30e6t . Then x + 3y = (e−2t + 3e6t ) + 3(−e−2t + 5e6t ) = −2e−2t + 18e6t = dx dt and 5x + 3y = 5(e−2t + 3e6t ) + 3(−e−2t + 5e6t ) = 2e−2t + 30e6t = dy dt . 50. From x = cos 2t + sin 2t + 1 5et and y = − cos 2t − sin 2t − 1 5 et we obtain dx dt = −2 sin 2t + 2 cos 2t + 1 5 et and dy dt = 2 sin 2t − 2 cos 2t − 1 5 et and d2x dt2 = −4 cos 2t − 4 sin 2t + 1 5 et and d2y dt2 = 4 cos 2t + 4 sin 2t − 1 5 et . Then 4y + et = 4(− cos 2t − sin 2t − 1 5 et ) + et = −4 cos 2t − 4 sin 2t + 1 5 et = d2x dt2 and 4x − et = 4(cos 2t + sin 2t + 1 5 et ) − et = 4 cos 2t + 4 sin 2t − 1 5 et = d2y dt2 . 51. (y)2 + 1 = 0 has no real solutions because (y)2 + 1 is positive for all differentiable functions y = φ(x). 52. The only solution of (y)2 + y2 = 0 is y = 0, since if y = 0, y2 0 and (y)2 + y2 ≥ y2 0.
  • 10. 10 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS 53. The first derivative of f(x) = ex is ex. The first derivative of f(x) = ekx is kekx. The differential equations are y = y and y = ky, respectively. 54. Any function of the form y = cex or y = ce−x is its own second derivative. The corresponding differential equation is y −y = 0. Functions of the form y = c sin x or y = c cos x have second derivatives that are the negatives of themselves. The differential equation is y + y = 0. 55. We first note that 1 − y2 = 1 − sin2 x = √ cos2 x = | cos x|. This prompts us to consider values of x for which cos x 0, such as x = π. In this case dy dx x=π = d dx (sin x) x=π = cos x x=π = cos π = −1, but 1 − y2|x=π = 1 − sin2 π = √ 1 = 1. Thus, y = sin x will only be a solution of y = 1 − y2 when cos x 0. An interval of definition is then (−π/2, π/2). Other intervals are (3π/2, 5π/2), (7π/2, 9π/2), and so on. 56. Since the first and second derivatives of sin t and cos t involve sin t and cos t, it is plausible that a linear combination of these functions, A sin t+B cos t, could be a solution of the differential equation. Using y = A cos t − B sin t and y = −A sin t − B cos t and substituting into the differential equation we get y + 2y + 4y = −A sin t − B cos t + 2A cos t − 2B sin t + 4A sin t + 4B cos t = (3A − 2B) sin t + (2A + 3B) cos t = 5 sin t Thus 3A − 2B = 5 and 2A + 3B = 0. Solving these simultaneous equations we find A = 15 13 and B = −10 13 . A particular solution is y = 15 13 sin t − 10 13 cos t. 57. One solution is given by the upper portion of the graph with domain approximately (0, 2.6). The other solution is given by the lower portion of the graph, also with domain approximately (0, 2.6). 58. One solution, with domain approximately (−∞, 1.6) is the portion of the graph in the second quadrant together with the lower part of the graph in the first quadrant. A second solution, with domain approximately (0, 1.6) is the upper part of the graph in the first quadrant. The third solution, with domain (0, ∞), is the part of the graph in the fourth quadrant.
  • 11. 1.1 Definitions and Terminology 11 59. Differentiating (x3 + y3)/xy = 3c we obtain xy(3x2 + 3y2y) − (x3 + y3)(xy + y) x2y2 = 0 3x3 y + 3xy3 y − x4 y − x3 y − xy3 y − y4 = 0 (3xy3 − x4 − xy3 )y = −3x3 y + x3 y + y4 y = y4 − 2x3y 2xy3 − x4 = y(y3 − 2x3) x(2y3 − x3) . 60. A tangent line will be vertical where y is undefined, or in this case, where x(2y3 − x3) = 0. This gives x = 0 and 2y3 = x3. Substituting y3 = x3/2 into x3 + y3 = 3xy we get x3 + 1 2 x3 = 3x 1 21/3 x 3 2 x3 = 3 21/3 x2 x3 = 22/3 x2 x2 (x − 22/3 ) = 0. Thus, there are vertical tangent lines at x = 0 and x = 22/3, or at (0, 0) and (22/3, 21/3). Since 22/3 ≈ 1.59, the estimates of the domains in Problem 58 were close. 61. The derivatives of the functions are φ 1(x) = −x/ √ 25 − x2 and φ 2(x) = x/ √ 25 − x2, neither of which is defined at x = ±5. 62. To determine if a solution curve passes through (0, 3) we let t = 0 and P = 3 in the equation P = c1et/(1 + c1et). This gives 3 = c1/(1 + c1) or c1 = −3 2 . Thus, the solution curve P = (−3/2)et 1 − (3/2)et = −3et 2 − 3et passes through the point (0, 3). Similarly, letting t = 0 and P = 1 in the equation for the one-parameter family of solutions gives 1 = c1/(1 + c1) or c1 = 1 + c1. Since this equation has no solution, no solution curve passes through (0, 1). 63. For the first-order differential equation integrate f(x). For the second-order differential equa- tion integrate twice. In the latter case we get y = ( f(x)dx) dx + c1x + c2. 64. Solving for y using the quadratic formula we obtain the two differential equations y = 1 x 2 + 2 1 + 3x6 and y = 1 x 2 − 2 1 + 3x6 , so the differential equation cannot be put in the form dy/dx = f(x, y).
  • 12. 12 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS 65. The differential equation yy − xy = 0 has normal form dy/dx = x. These are not equivalent because y = 0 is a solution of the first differential equation but not a solution of the second. 66. Differentiating we get y = c1 + 3c2x2 and y = 6c2x. Then c2 = y/6x and c1 = y − xy/2, so y = y − xy 2 x + y 6x x3 = xy − 1 3 x2 y and the differential equation is x2y − 3xy + 3y = 0. 67. (a) Since e−x2 is positive for all values of x, dy/dx 0 for all x, and a solution, y(x), of the differential equation must be increasing on any interval. (b) lim x→−∞ dy dx = lim x→−∞ e−x2 = 0 and lim x→∞ dy dx = lim x→∞ e−x2 = 0. Since dy/dx approaches 0 as x approaches −∞ and ∞, the solution curve has horizontal asymptotes to the left and to the right. (c) To test concavity we consider the second derivative d2y dx2 = d dx dy dx = d dx e−x2 = −2xe−x2 . Since the second derivative is positive for x 0 and negative for x 0, the solution curve is concave up on (−∞, 0) and concave down on (0, ∞). (d) x y 68. (a) The derivative of a constant solution y = c is 0, so solving 5 − c = 0 we see that c = 5 and so y = 5 is a constant solution. (b) A solution is increasing where dy/dx = 5 − y 0 or y 5. A solution is decreasing where dy/dx = 5 − y 0 or y 5. 69. (a) The derivative of a constant solution is 0, so solving y(a − by) = 0 we see that y = 0 and y = a/b are constant solutions. (b) A solution is increasing where dy/dx = y(a − by) = by(a/b − y) 0 or 0 y a/b. A solution is decreasing where dy/dx = by(a/b − y) 0 or y 0 or y a/b.
  • 13. 1.1 Definitions and Terminology 13 (c) Using implicit differentiation we compute d2y dx2 = y(−by ) + y (a − by) = y (a − 2by). Solving d2y/dx2 = 0 we obtain y = a/2b. Since d2y/dx2 0 for 0 y a/2b and d2y/dx2 0 for a/2b y a/b, the graph of y = φ(x) has a point of inflection at y = a/2b. (d) y = a/b y = 0 x y 70. (a) If y = c is a constant solution then y = 0, but c2 + 4 is never 0 for any real value of c. (b) Since y = y2 + 4 0 for all x where a solution y = φ(x) is defined, any solution must be increasing on any interval on which it is defined. Thus it cannot have any relative extrema. (c) Using implicit differentiation we compute d2y/dx2 = 2yy = 2y(y2 + 4). Setting d2y/dx2 = 0 we see that y = 0 corresponds to the only possible point of inflection. Since d2y/dx2 0 for y 0 and d2y/dx2 0 for y 0, there is a point of inflection where y = 0. (d) x y
  • 14. 14 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS 71. In Mathematica use Clear[y] y[x ]:= x Exp[5x] Cos[2x] y[x] y''''[x] − 20y'''[x] + 158y''[x] − 580y'[x] +841y[x]//Simplify The output will show y(x) = e5xx cos 2x, which verifies that the correct function was entered, and 0, which verifies that this function is a solution of the differential equation. 72. In Mathematica use Clear[y] y[x ]:= 20Cos[5Log[x]]/x − 3Sin[5Log[x]]/x y[x] xˆ3 y'''[x] + 2xˆ2 y''[x] + 20x y'[x] − 78y[x]//Simplify The output will show y(x) = 20 cos(5 ln x)/x − 3 sin(5 ln x)/x, which verifies that the correct function was entered, and 0, which verifies that this function is a solution of the differential equation. 1.2 Initial-Value Problems 1. Solving −1/3 = 1/(1 + c1) we get c1 = −4. The solution is y = 1/(1 − 4e−x ). 2. Solving 2 = 1/(1 + c1e) we get c1 = −(1/2)e−1. The solution is y = 2/(2 − e−(x+1)) . 3. Letting x = 2 and solving 1/3 = 1/(4 + c) we get c = −1. The solution is y = 1/(x2 − 1). This solution is defined on the interval (1, ∞). 4. Letting x = −2 and solving 1/2 = 1/(4 + c) we get c = −2. The solution is y = 1/(x2 − 2). This solution is defined on the interval (−∞, − √ 2 ). 5. Letting x = 0 and solving 1 = 1/c we get c = 1. The solution is y = 1/(x2 + 1). This solution is defined on the interval (−∞, ∞). 6. Letting x = 1/2 and solving −4 = 1/(1/4 + c) we get c = −1/2. The solution is y = 1/(x2 − 1/2) = 2/(2x2 − 1). This solution is defined on the interval (−1/ √ 2 , 1/ √ 2 ). In Problems 7–10, we use x = c1 cos t + c2 sin t and x = −c1 sin t + c2 cos t to obtain a system of two equations in the two unknowns c1 and c2.
  • 15. 1.2 Initial-Value Problems 15 7. From the initial conditions we obtain the system c1 = −1c2 = 8 The solution of the initial-value problem is x = − cos t + 8 sin t. 8. From the initial conditions we obtain the system c2 = 0 − c1 = 1 The solution of the initial-value problem is x = − cos t. 9. From the initial conditions we obtain √ 3 2 c1 + 1 2 c2 = 1 2 − 1 2 c2 + √ 3 2 = 0 Solving, we find c1 = √ 3/4 and c2 = 1/4. The solution of the initial-value problem is x = ( √ 3/4) cos t + (1/4) sin t. 10. From the initial conditions we obtain √ 2 2 c1 + √ 2 2 c2 = √ 2 − √ 2 2 c1 + √ 2 2 c2 = 2 √ 2. Solving, we find c1 = −1 and c2 = 3. The solution of the initial-value problem is x = − cos t + 3 sin t. In Problems 11–14, we use y = c1ex + c2e−x and y = c1ex − c2e−x to obtain a system of two equations in the two unknowns c1 and c2. 11. From the initial conditions we obtain c1 + c2 = 1 c1 − c2 = 2. Solving, we find c1 = 3 2 and c2 = −1 2 . The solution of the initial-value problem is y = 3 2ex − 1 2e−x. 12. From the initial conditions we obtain ec1 + e−1 c2 = 0 ec1 − e−1 c2 = e. Solving, we find c1 = 1 2 and c2 = −1 2e2. The solution of the initial-value problem is y = 1 2 ex − 1 2 e2 e−x = 1 2 ex − 1 2 e2−x .
  • 16. 16 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS 13. From the initial conditions we obtain e−1 c1 + ec2 = 5 e−1 c1 − ec2 = −5. Solving, we find c1 = 0 and c2 = 5e−1. The solution of the initial-value problem is y = 5e−1e−x = 5e−1−x. 14. From the initial conditions we obtain c1 + c2 = 0 c1 − c2 = 0. Solving, we find c1 = c2 = 0. The solution of the initial-value problem is y = 0. 15. Two solutions are y = 0 and y = x3. 16. Two solutions are y = 0 and y = x2. (Also, any constant multiple of x2 is a solution.) 17. For f(x, y) = y2/3 we have ∂f ∂y = 2 3 y−1/3 . Thus, the differential equation will have a unique solution in any rectangular region of the plane where y = 0. 18. For f(x, y) = √ xy we have ∂f/∂y = 1 2 x/y . Thus the differential equation will have a unique solution in any region where x 0 and y 0 or where x 0 and y 0. 19. For f(x, y) = y x we have ∂f ∂y = 1 x . Thus, the differential equation will have a unique solution in any region where x = 0. 20. For f(x, y) = x + y we have ∂f ∂y = 1. Thus, the differential equation will have a unique solution in the entire plane. 21. For f(x, y) = x2/(4 − y2) we have ∂f/∂y = 2x2y/(4 − y2)2. Thus the differential equation will have a unique solution in any region where y −2, −2 y 2, or y 2. 22. For f(x, y) = x2 1 + y3 we have ∂f ∂y = −3x2y2 (1 + y3)2 . Thus, the differential equation will have a unique solution in any region where y = −1. 23. For f(x, y) = y2 x2 + y2 we have ∂f ∂y = 2x2y (x2 + y2)2 . Thus, the differential equation will have a unique solution in any region not containing (0, 0).
  • 17. 1.2 Initial-Value Problems 17 24. For f(x, y) = (y + x)/(y − x) we have ∂f/∂y = −2x/(y − x)2. Thus the differential equation will have a unique solution in any region where y x or where y x. In Problems 25–28, we identify f(x, y) = y2 − 9 and ∂f/∂y = y/ y2 − 9. We see that f and ∂f/∂y are both continuous in the regions of the plane determined by y −3 and y 3 with no restrictions on x. 25. Since 4 3, (1, 4) is in the region defined by y 3 and the differential equation has a unique solution through (1, 4). 26. Since (5, 3) is not in either of the regions defined by y −3 or y 3, there is no guarantee of a unique solution through (5, 3). 27. Since (2, −3) is not in either of the regions defined by y −3 or y 3, there is no guarantee of a unique solution through (2, −3). 28. Since (−1, 1) is not in either of the regions defined by y −3 or y 3, there is no guarantee of a unique solution through (−1, 1). 29. (a) A one-parameter family of solutions is y = cx. Since y = c, xy = xc = y and y(0) = c · 0 = 0. (b) Writing the equation in the form y = y/x, we see that R cannot contain any point on the y-axis. Thus, any rectangular region disjoint from the y-axis and containing (x0, y0) will determine an interval around x0 and a unique solution through (x0, y0). Since x0 = 0 in part (a), we are not guaranteed a unique solution through (0, 0). (c) The piecewise-defined function which satisfies y(0) = 0 is not a solution since it is not differentiable at x = 0. 30. (a) Since d dx tan (x + c) = sec2 (x + c) = 1+tan2 (x + c), we see that y = tan (x + c) satisfies the differential equation. (b) Solving y(0) = tan c = 0 we obtain c = 0 and y = tan x. Since tan x is discontinuous at x = ±π/2, the solution is not defined on (−2, 2) because it contains ±π/2. (c) The largest interval on which the solution can exist is (−π/2, π/2). 31. (a) Since d dx − 1 x + c = 1 (x + c)2 = y2 , we see that y = − 1 x + c is a solution of the differ- ential equation. (b) Solving y(0) = −1/c = 1 we obtain c = −1 and y = 1/(1−x). Solving y(0) = −1/c = −1 we obtain c = 1 and y = −1/(1+x). Being sure to include x = 0, we see that the interval of existence of y = 1/(1−x) is (−∞, 1), while the interval of existence of y = −1/(1+x) is (−1, ∞).
  • 18. 18 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS 32. (a) Solving y(0) = −1/c = y0 we obtain c = −1/y0 and y = − 1 −1/y0 + x = y01 − y0x , y0 = 0 Since we must have −1/y0 + x = 0, the largest interval of existence (which must contain 0) is either (−∞, 1/y0) when y0 0 or (1/y0, ∞) when y0 0. (b) By inspection we see that y = 0 is a solution on (−∞, ∞). 33. (a) Differentiating 3x2 − y2 = c we get 6x − 2yy = 0 or yy = 3x. (b) Solving 3x2 − y2 = 3 for y we get y = φ1(x) = 3(x2 − 1) , 1 x ∞, y = φ2(x) = − 3(x2 − 1) , 1 x ∞, y = φ3(x) = 3(x2 − 1) , −∞ x −1, y = φ4(x) = − 3(x2 − 1) , −∞ x −1. –4 –2 2 4 x –4 –2 2 4 y (c) Only y = φ3(x) satisfies y(−2) = 3. 34. (a) Setting x = 2 and y = −4 in 3x2 − y2 = c we get 12 − 16 = −4 = c, so the explicit solution is y = − 3x2 + 4 , −∞ x ∞. –4 –2 2 4 x –4 –2 2 4 y (b) Setting c = 0 we have y = √ 3x and y = − √ 3x, both defined on (−∞, ∞). In Problems 35–38, we consider the points on the graphs with x-coordinates x0 = −1, x0 = 0, and x0 = 1. The slopes of the tangent lines at these points are compared with the slopes given by y(x0) in (a) through (f). 35. The graph satisfies the conditions in (b) and (f). 36. The graph satisfies the conditions in (e). 37. The graph satisfies the conditions in (c) and (d).
  • 19. 1.2 Initial-Value Problems 19 38. The graph satisfies the conditions in (a). 39. Using the function y = c1 cos 3x + c2 sin 3x and the first boundary condition we get y(0) = c1 cos 0 + c2 sin 0 = 0 Therefore c1 = 0. Similarly, for the second boundary condition we get y(π/6) = c2 sin 3 (π/6) = −1 Therefore c2 = −1. The solution to the boundary value problem is y(x) = − sin 3x. 40. Using the function y = c1 cos 3x + c2 sin 3x and the first boundary condition we get y(0) = c1 cos 0 + c2 sin 0 = 0 Therefore c1 = 0. Thus y(x) = c2 sin 3x. Similarly for the second boundary condition we get y(π) = c2 sin 3(π) = 0 But the last line is satisfied for any choice of c2 since sin 3π = 0 therefore there are infinitely many solutions in this case. 41. The derivative of the function y = c1 cos 3x + c2 sin 3x is y = −3c1 sin 3x + 3c2 cos 3x and using the first boundary condition we get y (0) = 0 + 3c2 = 0 Therefore c2 = 0 So far then we have y = −3c1 sin 3x. Similarly, using the second boundary condition we get y (π/4) = −3c1 √ 2/2 = 0 Therefore c1 = 0. The solution is thus the trivial solution y = 0. 42. The derivative of the function y = c1 cos 3x + c2 sin 3x is y = −3c1 sin 3x + 3c2 cos 3x and using the two boundary conditions we get y(0) = c1 + 0 = 1 Therefore c1 = 1. In addition y (π) = 0 − 3c2 = 5 Therefore c2 = −5/3. The solution to this boundary value problem is y(x) = cos 3x− 5 3 sin 3x.
  • 20. 20 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS 43. Using the function y = c1 cos 3x + c2 sin 3x and the first boundary condition we get y(0) = c1 cos 0 + c2 sin 0 = 0 Therefore c1 = 0. Thus y(x) = c2 sin 3x. Similarly for the second boundary condition we get c2 sin 3(π) = 4c2 · 0 = 40 = 4 The last line is obviously a contradiction and so therefore the boundary value problem has no solution in this case. 44. The derivative of the function y = c1 cos 3x + c2 sin 3x is y = −3c1 sin 3x + 3c2 cos 3x and using the first boundary condition we get y (π/3) = 0 − 3c2 = 1 Therefore c2 = −1/3. Similarly, using the second boundary condition we get y (π) = 0 − 3c2 = 0 Therefore c2 = 0, which is a contradiction so the problem has no solution. 45. Integrating y = 8e2x + 6x we obtain y = (8e2x + 6x) dx = 4e2x + 3x2 + c. Setting x = 0 and y = 9 we have 9 = 4 + c so c = 5 and y = 4e2x + 3x2 + 5. 46. Integrating y = 12x − 2 we obtain y = (12x − 2) dx = 6x2 − 2x + c1. Then, integrating y we obtain y = (6x2 − 2x + c1) dx = 2x3 − x2 + c1x + c2. At x = 1 the y-coordinate of the point of tangency is y = −1 + 5 = 4. This gives the initial condition y(1) = 4. The slope of the tangent line at x = 1 is y(1) = −1. From the initial conditions we obtain 2 − 1 + c1 + c2 = 4 or c1 + c2 = 3 and 6 − 2 + c1 = −1 or c1 = −5. Thus, c1 = −5 and c2 = 8, so y = 2x3 − x2 − 5x + 8.
  • 21. 1.2 Initial-Value Problems 21 47. When x = 0 and y = 1 2 , y = −1, so the only plausible solution curve is the one with negative slope at (0, 1 2 ), or the blue curve. 48. If y(x) is a solution of the initial-value problem y = x2+y2. y(0) = 1, then y(0) = 02+12 = 1. Differentiating the differential equation, we have y = 2x+2yy. So y(0) = 2·0+2·1·1 = 2. 49. If the solution is tangent to the x-axis at (x0, 0), then y = 0 when x = x0 and y = 0. Substituting these values into y + 2y = 3x − 6 we get 0 + 0 = 3x0 − 6 or x0 = 2. 50. When y = 1 16 x4, y = 1 4x3 = x(1 4 x2) = xy1/2, and y(2) = 1 16(16) = 1. When y = ⎧ ⎪ ⎨ ⎪ ⎩ 0, x 0 1 16 x4, x ≥ 0 we have y = ⎧ ⎨ ⎩ 0, x 0 1 4x3, x ≥ 0 = x ⎧ ⎨ ⎩ 0, x 0 1 4x2, x ≥ 0 = xy1/2 , and y(2) = 1 16 (16) = 1. The two different solutions are the same on the interval (0, ∞), which is all that is required by Theorem 1.2.1. 51. We note that the initial condition y(0) = 0, 0 = y 0 1 √ t3 + 1 dt is satisfied only when y = 0. For any y 0, necessarily y 0 1 √ t3 + 1 dt 0 because the integrand is positive on the interval of integration. Then from (12) of Section 1.1 and the Chain Rule we have: d dx x = d dx y 0 1 √ t3 + 1 dt 1 = 1 y3 + 1 dy dx and dy dx = y3 + 1 y (0) = dy dx x=0 = (y(0))3 + 1 = √ 0 + 1 = 1. Computing the second derivative, we see that: d2y dx2 = d dx y3 + 1 = 3y2 2 y3 + 1 dy dx = 3y2 2 y3 + 1 · y3 + 1 = 3 2 y2 d2y dx2 = 3 2 y2 . This is equivalent to 2 d2y dx2 − 3y2 = 0.
  • 22. 22 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS 1.3 Differential Equations as Mathematical Models 1. dP dt = kP + r; dP dt = kP − r 2. Let b be the rate of births and d the rate of deaths. Then b = k1P and d = k2P. Since dP/dt = b − d, the differential equation is dP/dt = k1P − k2P. 3. Let b be the rate of births and d the rate of deaths. Then b = k1P and d = k2P2. Since dP/dt = b − d, the differential equation is dP/dt = k1P − k2P2. 4. dP dt = k1P − k2P2 − h, h 0 5. From the graph in the text we estimate T0 = 180◦ and Tm = 75◦. We observe that when T = 85, dT/dt ≈ −1. From the differential equation we then have k = dT/dt T − Tm = −1 85 − 75 = −0.1. 6. By inspecting the graph in the text we take Tm to be Tm(t) = 80 − 30 cos (πt/12). Then the temperature of the body at time t is determined by the differential equation dT dt = k T − 80 − 30 cos π 12 t , t 0. 7. The number of students with the flu is x and the number not infected is 1000 − x, so dx/dt = kx(1000 − x). 8. By analogy, with the differential equation modeling the spread of a disease, we assume that the rate at which the technological innovation is adopted is proportional to the number of people who have adopted the innovation and also to the number of people, y(t), who have not yet adopted it. If one person who has adopted the innovation is introduced into the population, then x + y = n + 1 and dx dt = kx(n + 1 − x), x(0) = 1. 9. The rate at which salt is leaving the tank is Rout (3 gal/min) · A 300 lb/gal = A 100 lb/min. Thus dA/dt = A/100. The initial amount is A(0) = 50.
  • 23. 1.3 Differential Equations as Mathematical Models 23 10. The rate at which salt is entering the tank is Rin = (3 gal/min) · (2 lb/gal) = 6 lb/min. Since the solution is pumped out at a slower rate, it is accumulating at the rate of (3 − 2) gal/min = 1 gal/min. After t minutes there are 300 + t gallons of brine in the tank. The rate at which salt is leaving is Rout = (2 gal/min) · A 300 + t lb/gal = 2A 300 + t lb/min. The differential equation is dA dt = 6 − 2A 300 + t . 11. The rate at which salt is entering the tank is Rin = (3 gal/min) · (2 lb/gal) = 6 lb/min. Since the tank loses liquid at the net rate of 3 gal/min − 3.5 gal/min = −0.5 gal/min, after t minutes the number of gallons of brine in the tank is 300 − 1 2t gallons. Thus the rate at which salt is leaving is Rout = A 300 − t/2 lb/gal · (3.5 gal/min) = 3.5A 300 − t/2 lb/min = 7A 600 − t lb/min. The differential equation is dA dt = 6 − 7A 600 − t or dA dt + 7 600 − t A = 6. 12. The rate at which salt is entering the tank is Rin = (cin lb/gal) · (rin gal/min) = cinrin lb/min. Now let A(t) denote the number of pounds of salt and N(t) the number of gallons of brine in the tank at time t. The concentration of salt in the tank as well as in the outflow is c(t) = x(t)/N(t). But the number of gallons of brine in the tank remains steady, is increased, or is decreased depending on whether rin = rout, rin rout, or rin rout. In any case, the number of gallons of brine in the tank at time t is N(t) = N0 + (rin − rout)t. The output rate of salt is then Rout = A N0 + (rin − rout)t lb/gal · (rout gal/min) = rout A N0 + (rin − rout)t lb/min. The differential equation for the amount of salt, dA/dt = Rin − Rout, is dA dt = cinrin − rout A N0 + (rin − rout)t or dA dt + rout N0 + (rin − rout)t A = cinrin.
  • 24. 24 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS 13. The volume of water in the tank at time t is V = Awh. The differential equation is then dh dt = 1 Aw dV dt = 1 Aw −cAh 2gh = − cAh Aw 2gh . Using Ah = π 2 12 2 = π 36 , Aw = 102 = 100, and g = 32, this becomes dh dt = − cπ/36 100 √ 64h = − cπ 450 √ h . 14. The volume of water in the tank at time t is V = 1 3 πr2h where r is the radius of the tank at height h. From the figure in the text we see that r/h = 8/20 so that r = 2 5h and V = 1 3π 2 5h 2 h = 4 75 πh3. Differentiating with respect to t we have dV/dt = 4 25πh2 dh/dt or dh dt = 25 4πh2 dV dt . From Problem 13 we have dV/dt = −cAh √ 2gh where c = 0.6, Ah = π 2 12 2 , and g = 32. Thus dV/dt = −2π √ h/15 and dh dt = 25 4πh2 − 2π √ h 15 = − 5 6h3/2 . 15. Since i = dq/dt and L d2q/dt2 + R dq/dt = E(t), we obtain L di/dt + Ri = E(t). 16. By Kirchhoff’s second law we obtain R dq dt + 1 C q = E(t). 17. From Newton’s second law we obtain m dv dt = −kv2 + mg. 18. Since the barrel in Figure 1.3.17(b) in the text is submerged an additional y feet below its equilibrium position, the number of cubic feet in the additional submerged portion is the volume of the circular cylinder: π × (radius)2 × height or π(s/2)2y. Then we have from Archimedes’ principle Upward force of water on barrel = Weight of water displaced = (62.4) × (Volume of water displaced) = (62.4)π(s/2)2 y = 15.6πs2 y. It then follows from Newton’s second law that w g d2y dt2 = −15.6πs2 y or d2y dt2 + 15.6πs2g w y = 0, where g = 32 and w is the weight of the barrel in pounds.
  • 25. 1.3 Differential Equations as Mathematical Models 25 19. The net force acting on the mass is F = ma = m d2x dt2 = −k(s + x) + mg = −kx + mg − ks. Since the condition of equilibrium is mg = ks, the differential equation is m d2x dt2 = −kx. 20. From Problem 19, without a damping force, the differential equation is m d2x/dt2 = −kx. With a damping force proportional to velocity, the differential equation becomes m d2x dt2 = −kx − β dx dt or m d2x dt2 + β dx dt + kx = 0. 21. As the rocket climbs (in the positive direction), it spends its amount of fuel and therefore the mass of the fuel changes with time. The air resistance acts in the opposite direction of the motion and the upward thrust R works in the same direction. Using Newton’s second law we get d dt (mv) = −mg − kv + R Now because the mass is variable, we must use the product rule to expand the left side of the equation. Doing so gives us the following: d dt (mv) = −mg − kv + R v × dm dt + m × dv dt = −mg − kv + R The last line is the differential equation we wanted to find. 22. (a) Since the mass of the rocket is m(t) = mp + mv + mf (t), take the time rate-of-change and get, by straight-forward calculation, d dt m(t) = d dt (mp + mv + mf (t)) = 0 + 0 + m f (t) = d dt mf (t) Therefore the rate of change of the mass of the rocket is the same as the rate of change of the mass of the fuel which is what we wanted to show. (b) The fuel is decreasing at the constant rate of λ and so from part (a) we have d dt m(t) = d dt mf (t) = −λ m(t) = −λt + c
  • 26. 26 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS Using the given condition to solve for c, m(0) = 0 + c = m0 and so m(t) = −λt + m0. The differential equation in Problem 21 now becomes v dm dt + m dv dt + kv = −mg + R −λv + (−λt + m0) dv dt + kv = −mg + R (−λt + m0) dv dt + (k − λ)v = −mg + R dv dt + k − λ −λt + m0 v = −mg −λt + m0 + R −λt + m0 dv dt + k − λ −λt + m0 v = −g + R −λt + m0 (c) From part (b) we have that d dtmf (t) = −λ and so by integrating this result we get mf (t) = −λt+c. Now at time t = 0, mf (0) = 0+c = c therefore mf (t) = −λt+mf (0) . At some later time tb we then have mf (tb) = −λtb +mf (0) = 0 and solving this equation for that time we get tb = mf (0)/λ which is what we wanted to show. 23. From g = k/R2 we find k = gR2. Using a = d2r/dt2 and the fact that the positive direction is upward we get d2r dt2 = −a = − k r2 = − gR2 r2 or d2r dt2 + gR2 r2 = 0. 24. The gravitational force on m is F = −kMrm/r2. Since Mr = 4πδr3/3 and M = 4πδR3/3, we have Mr = r3M/R3 and F = −k Mrm r2 = −k r3Mm/R3 r2 = −k mM R3 r. Now from F = ma = d2r/dt2 we have m d2r dt2 = −k mM R3 r or d2r dt2 = − kM R3 r. 25. The differential equation is dA dt = k(M − A). 26. The differential equation is dA dt = k1(M − A) − k2A. 27. The differential equation is x(t) = r − kx(t) where k 0.
  • 27. 1.3 Differential Equations as Mathematical Models 27 28. Notice from the diagram that the segment from the waterskier to the boat is tangent to the curve at the point P(x, y). For the sake of simplicity, let’s label the coordinates of the boat on the x-axis as (a, 0). The equation of this tangent line is y = y(x)(x − a) from which we get y = y (x)(x − a) y y = x − a − y y = a − x By the Pythagorean Theorem we have y2 + (a − x)2 = s2 from which we make a substitution and then solve for the derivative to get y2 + (a − x)2 = s2 y2 + (y y )2 = s2 (y y )2 = s2 − y2 y y = − s2 − y2 y = − y s2 − y2 Notice that the sign of the derivative is negative because as the boat proceeds along the positive x-axis, the y-coordinate decreases. 29. We see from the figure that 2θ + α = π. Thus y −x = tan α = tan(π − 2θ) = − tan 2θ = − 2 tan θ 1 − tan2 θ . Since the slope of the tangent line is y = tan θ we have y/x = 2y[1−(y)2] or y−y(y)2 = 2xy, which is the quadratic equation y(y)2 + 2xy − y = 0 in y. Using the quadratic formula, we get y = −2x ± 4x2 + 4y2 2y = −x ± x2 + y2 y . (x, y) x y α α θ θ θ φ x y Since dy/dx 0, the differential equation is dy dx = −x + x2 + y2 y or y dy dx − x2 + y2 + x = 0.
  • 28. 28 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS 30. The differential equation is dP/dt = kP, so from Problem 49 in Exercises 1.1, a one-parameter family of solutions is P = cekt. 31. The differential equation in (3) is dT/dt = k(T − Tm). When the body is cooling, T Tm, so T − Tm 0. Since T is decreasing, dT/dt 0 and k 0. When the body is warming, T Tm, so T − Tm 0. Since T is increasing, dT/dt 0 and k 0. 32. The differential equation in (8) is dA/dt = 6 − A/100. If A(t) attains a maximum, then dA/dt = 0 at this time and A = 600. If A(t) continues to increase without reaching a maximum, then A(t) 0 for t 0 and A cannot exceed 600. In this case, if A(t) approaches 0 as t increases to infinity, we see that A(t) approaches 600 as t increases to infinity. 33. This differential equation could describe a population that undergoes periodic fluctuations. 34. (a) As shown in Figure 1.3.23(a) in the text, the resultant of the reaction force of magnitude F and the weight of magnitude mg of the particle is the centripetal force of magnitude mω2x. The centripetal force points to the center of the circle of radius x on which the particle rotates about the y-axis. Comparing parts of similar triangles gives F cos θ = mg and F sin θ = mω2 x. (b) Using the equations in part (a) we find tan θ = F sin θ F cos θ = mω2x mg = ω2x g or dy dx = ω2x g . 35. From Problem 23, d2r/dt2 = −gR2/r2. Since R is a constant, if r = R + s, then d2r/dt2 = d2s/dt2 and, using a Taylor series, we get d2s dt2 = −g R2 (R + s)2 = −gR2 (R + s)−2 ≈ −gR2 [R−2 − 2sR−3 + · · · ] = −g + 2gs R + · · · . Thus, for R much larger than s, the differential equation is approximated by d2s/dt2 = −g. 36. (a) If ρ is the mass density of the raindrop, then m = ρV and dm dt = ρ dV dt = ρ d dt 4 3 πr3 = ρ 4πr2 dr dt = ρS dr dt . If dr/dt is a constant, then dm/dt = kS where ρ dr/dt = k or dr/dt = k/ρ. Since the radius is decreasing, k 0. Solving dr/dt = k/ρ we get r = (k/ρ)t+c0. Since r(0) = r0, c0 = r0 and r = kt/ρ + r0. (b) From Newton’s second law, d dt [mv] = mg, where v is the velocity of the raindrop. Then m dv dt + v dm dt = mg or ρ 4 3 πr3 dv dt + v(k4πr2 ) = ρ 4 3 πr3 g. Dividing by 4ρπr3/3 we get dv dt + 3k ρr v = g or dv dt + 3k/ρ kt/ρ + r0 v = g, k 0.
  • 29. 1.3 Differential Equations as Mathematical Models 29 37. We assume that the plow clears snow at a constant rate of k cubic miles per hour. Let t be the time in hours after noon, x(t) the depth in miles of the snow at time t, and y(t) the distance the plow has moved in t hours. Then dy/dt is the velocity of the plow and the assumption gives wx dy dt = k, where w is the width of the plow. Each side of this equation simply represents the volume of snow plowed in one hour. Now let t0 be the number of hours before noon when it started snowing and let s be the constant rate in miles per hour at which x increases. Then for t −t0, x = s(t + t0). The differential equation then becomes dy dt = k ws 1 t + t0 . Integrating, we obtain y = k ws [ ln(t + t0) + c ] where c is a constant. Now when t = 0, y = 0 so c = − ln t0 and y = k ws ln 1 + t t0 . Finally, from the fact that when t = 1, y = 2 and when t = 2, y = 3, we obtain 1 + 2 t0 2 = 1 + 1 t0 3 . Expanding and simplifying gives t2 0 + t0 − 1 = 0. Since t0 0, we find t0 ≈ 0.618 hours ≈ 37 minutes. Thus it started snowing at about 11:23 in the morning. 38. . (1) : dP dt = kP is linear (2) : dA dt = kA is linear (3) : dT dt = k (T − Tm) is linear (5) : dx dt = kx (n + 1 − x) is nonlinear (6) : dX dt = k (α − X) (β − X) is nonlinear (8) : dA dt = 6 − A 100 is linear (10) : dh dt = − Ah Aw 2gh is nonlinear (11) : L d2q dt2 + R dq dt + 1 C q = E(t) is linear (12) : d2s dt2 = −g is linear (14) : m dv dt = mg − kv is linear (15) : m d2s dt2 + k ds dt = mg is linear (16) : d2x dt2 − 64 L x = 0 is linear (17) : linearity or nonlinearity is determined by the manner in which W and T1 involve x. 39. At time t, when the population is 2 million cells, the differential equation P(t) = 0.15P(t) gives the rate of increase at time t. Thus, when P(t) = 2 (million cells), the rate of increase is P(t) = 0.15(2) = 0.3 million cells per hour or 300,000 cells per hour.
  • 30. 30 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS 40. Setting A(t) = −0.002 and solving A(t) = −0.0004332A(t) for A(t), we obtain A(t) = A(t) −0.0004332 = −0.002 −0.0004332 ≈ 4.6 grams. Chapter 1 in Review 1. d dx c1ekx = k y c1ekx ; dy dx = ky 2. d dx (5 + c1e−2x ) = −2c1e−2x = −2( y 5 + c1e−2x −5); dy dx = −2(y − 5) or dy dx = −2y + 10 3. d dx (c1 cos kx + c2 sin kx) = −kc1 sin kx + kc2 cos kx; d2 dx2 (c1 cos kx + c2 sin kx) = −k2 c1 cos kx − k2 c2 sin kx = −k2 ( y c1 cos kx + c2 sin kx); d2y dx2 = −k2 y or d2y dx2 + k2 y = 0 4. d dx (c1 cosh kx + c2 sinh kx) = kc1 sinh kx + kc2 cosh kx; d2 dx2 (c1 cosh kx + c2 sinh kx) = k2 c1 cosh kx + k2 c2 sinh kx = k2 ( y c1 cosh kx + c2 sinh kx); d2y dx2 = k2 y or d2y dx2 − k2 y = 0 5. y = c1ex + c2xex; y = c1ex + c2xex + c2ex; y = c1ex + c2xex + 2c2ex; y + y = 2(c1ex + c2xex) + 2c2ex = 2(c1ex + c2xex + c2ex) = 2y; y − 2y + y = 0 6. y = −c1ex sin x + c1ex cos x + c2ex cos x + c2ex sin x; y = −c1ex cos x − c1ex sin x − c1ex sin x + c1ex cos x − c2ex sin x + c2ex cos x + c2ex cos x + c2ex sin x = −2c1ex sin x + 2c2ex cos x; y − 2y = −2c1ex cos x − 2c2ex sin x = −2y; y − 2y + 2y = 0 7. a, d 8. c 9. b 10. a, c 11. b 12. a, b, d 13. A few solutions are y = 0, y = c, and y = ex. 14. Easy solutions to see are y = 0 and y = 3. 15. The slope of the tangent line at (x, y) is y, so the differential equation is y = x2 + y2.
  • 31. Chapter 1 in Review 31 16. The rate at which the slope changes is dy/dx = y, so the differential equation is y = −y or y + y = 0. 17. (a) The domain is all real numbers. (b) Since y = 2/3x1/3, the solution y = x2/3 is undefined at x = 0. This function is a solution of the differential equation on (−∞, 0) and also on (0, ∞). 18. (a) Differentiating y2 − 2y = x2 − x + c we obtain 2yy − 2y = 2x − 1 or (2y − 2)y = 2x − 1. (b) Setting x = 0 and y = 1 in the solution we have 1 − 2 = 0 − 0 + c or c = −1. Thus, a solution of the initial-value problem is y2 − 2y = x2 − x − 1. (c) Solving the equation y2 − 2y − (x2 − x − 1) = 0 by the quadratic formula we get y = (2 ± 4 + 4(x2 − x − 1) )/2 = 1 ± √ x2 − x = 1 ± x(x − 1) . Since x(x − 1) ≥ 0 for x ≤ 0 or x ≥ 1, we see that neither y = 1 + x(x − 1) nor y = 1 − x(x − 1) is differentiable at x = 0. Thus, both functions are solutions of the differential equation, but neither is a solution of the initial-value problem. 19. Setting x = x0 and y = 1 in y = −2/x + x, we get 1 = − 2 x0 + x0 or x2 0 − x0 − 2 = (x0 − 2)(x0 + 1) = 0. Thus, x0 = 2 or x0 = −1. Since x = 0 in y = −2/x + x, we see that y = −2/x + x is a solution of the initial-value problem xy + y = 2x, y(−1) = 1, on the interval (−∞, 0) and y = −2/x + x is a solution of the initial-value problem xy + y = 2x, y(2) = 1, on the interval (0, ∞). 20. From the differential equation, y(1) = 12+[y(1)]2 = 1+(−1)2 = 2 0, so y(x) is increasing in some neighborhood of x = 1. From y = 2x+2yy we have y(1) = 2(1)+2(−1)(2) = −2 0, so y(x) is concave down in some neighborhood of x = 1. 21. (a) –3 –2 1 2 3 x –3 –2 1 2 3 y –1 –1 –3 –2 1 2 3 x –3 –2 1 y = x2 + c1 y = –x2 + c2 2 3 y –1 –1 (b) When y = x2 + c1, y = 2x and (y)2 = 4x2. When y = −x2 + c2, y = −2x and (y)2 = 4x2.
  • 32. 32 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS (c) Pasting together x2, x ≥ 0, and −x2, x ≤ 0, we get y = ⎧ ⎨ ⎩ −x2, x ≤ 0 x2, x 0 . 22. The slope of the tangent line is y (−1,4) = 6 √ 4 + 5(−1)3 = 7. 23. Differentiating y = x sin x + x cos x we get y = x cos x + sin x − x sin x + cos x and y = −x sin x + cos x + cos x − x cos x − sin x − sin x = −x sin x − x cos x + 2 cos x − 2 sin x. Thus y + y = −x sin x − x cos x + 2 cos x − 2 sin x + x sin x + x cos x = 2 cos x − 2 sin x. An interval of definition for the solution is (−∞, ∞). 24. Differentiating y = x sin x + (cos x) ln(cos x) we get y = x cos x + sin x + cos x − sin x cos x − (sin x) ln (cos x) = x cos x + sin x − sin x − (sin x) ln (cos x) = x cos x − (sin x) ln (cos x) and y = −x sin x + cos x − sin x − sin x cos x − (cos x) ln (cos x) = −x sin x + cos x + sin2 x cos x − (cos x) ln (cos x) = −x sin x + cos x + 1 − cos2 x cos x − (cos x) ln (cos x) = −x sin x + cos x + sec x − cos x − (cos x) ln (cos x) = −x sin x + sec x − (cos x) ln (cos x). Thus y + y = −x sin x + sec x − (cos x) ln(cos x) + x sin x + (cos x) ln (cos x) = sec x. To obtain an interval of definition we note that the domain of ln x is (0, ∞), so we must have cos x 0. Thus, an interval of definition is (−π/2, π/2).
  • 33. Chapter 1 in Review 33 25. Differentiating y = sin (ln x) we obtain y = cos (ln x)/x and y = −[sin (ln x) + cos (ln x)]/x2. Then x2 y + xy + y = x2 − sin (ln x) + cos (ln x) x2 + x cos (ln x) x + sin (ln x) = 0. An interval of definition for the solution is (0, ∞). 26. Differentiating y = cos (ln x) ln (cos (ln x)) + (ln x) sin (ln x) we obtain y = cos (ln x) 1 cos (ln x) − sin (ln x) x + ln (cos (ln x)) − sin (ln x) x + ln x cos (ln x) x + sin (ln x) x = − ln (cos (ln x)) sin (ln x) x + (ln x) cos (ln x) x and y = −x ln (cos (ln x)) cos (ln x) x + sin (ln x) 1 cos (ln x) − sin (ln x) x 1 x2 + ln (cos (ln x)) sin (ln x) 1 x2 + x (ln x) − sin (ln x) x + cos (ln x) x 1 x2 − (ln x) cos (ln x) 1 x2 = 1 x2 − ln (cos (ln x)) cos (ln x) + sin2 (ln x) cos (ln x) + ln (cos (ln x)) sin (ln x) − (ln x) sin (ln x) + cos (ln x) − (ln x) cos (ln x) . Then x2 y + xy + y = − ln (cos (ln x)) cos (ln x) + sin2 (ln x) cos (ln x) + ln (cos (ln x)) sin (ln x) − (ln x) sin (ln x) + cos (ln x) − (ln x) cos (ln x) − ln (cos (ln x)) sin (ln x) + (ln x) cos (ln x) + cos (ln x) ln (cos (ln x)) + (ln x) sin (ln x) = sin2 (ln x) cos (ln x) + cos (ln x) = sin2 (ln x) + cos2 (ln x) cos (ln x) = 1 cos (ln x) = sec (ln x). To obtain an interval of definition, we note that the domain of ln x is (0, ∞), so we must have cos (ln x) 0. Since cos x 0 when −π/2 x π/2, we require −π/2 ln x π/2. Since ex is an increasing function, this is equivalent to e−π/2 x eπ/2. Thus, an interval of definition is (e−π/2, eπ/2). (Much of this problem is more easily done using a computer algebra system such as Mathematica or Maple.)
  • 34. 34 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS In Problems 27–30 we use (12) of Section 1.1 and the Product Rule. 27. . y = ecos x x 0 te− cos t dt dy dx = ecos x xe− cos x − sin xecos x x 0 te− cos t dt dy dx + (sin x) y = ecos x xe− cos x − sin xecos x x 0 te− cos t dt + sin x ecos x x 0 te− cos t dt = x − sin xecos x x 0 te− cos t dt + sin xecos x x 0 te− cos t dt = x 28. . y = ex2 x 0 et−t2 dt dy dx = ex2 ex−x2 + 2xex2 x 0 et−t2 dt dy dx − 2xy = ex2 ex−x2 + 2xex2 x 0 et−t2 dt − 2x ex2 x 0 et−t2 dt = ex 29. . y = x x 1 e−t t dt y = x e−x x + x 1 e−t t dt = e−x + x 1 e−t t dt y = −e−x + e−x x x2 y + x2 − x y + (1 − x) y = −x2 e−x + xe−x + x2 e−x + x2 x 1 e−t t dt − xe−x − x x 1 e−t t dt + x x 1 e−t t dt − x2 x 1 e−t t dt = 0
  • 35. Chapter 1 in Review 35 30. . y = sin x x 0 et2 cos t dt − cos x x 0 et2 sin t dt y = sin x ex2 cos x + cos x x 0 et2 cos t dt − cos x ex2 sin x + sin x x 0 et2 sin t dt = cos x x 0 et2 cos t dt + sin x x 0 et2 sin t dt y = cos x ex2 cos x − sin x x 0 et2 cos t dt + sin x ex2 sin x + cos x x 0 et2 sin t dt = ex2 cos2 x + sin2 x − ⎛ ⎜ ⎜ ⎝ y sin x x 0 et2 cos t dt − cos x x 0 et2 sin t dt ⎞ ⎟ ⎟ ⎠ = ex2 − y y + y = ex2 − y + y = ex2 31. Using implicit differentiation we get x3 y3 = x3 + 5 3x2 · y3 + x3 · 3y2 dy dx = 3x2 3x2y3 3x2y2 + x33y2 3x2y2 dy dx = 3x2 3x2y2 y + x dy dx = 1 y2 32. Using implicit differentiation we get (x − 7)2 + y2 = 1 2(x − 7) + 2y dy dx = 0 2y dy dx = −2(x − 7) y dy dx = −(x − 7) y dy dx 2 = (x − 7)2 dy dx 2 = (x − 7)2 y2
  • 36. 36 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS Now from the original equation, isolating the first term leads to (x−7)2 = 1−y2. Continuing from the last line of our proof we now have dy dx 2 = (x − 7)2 y2 = 1 − y2 y2 = 1 y2 − 1 Adding 1 to both sides leads to the desired result. 33. Using implicit differentiation we get y3 + 3y = 2 − 3x 3y2 y + 3y = −3 y2 y + y = −1 (y2 + 1)y = −1 y = −1 y2 + 1 Differentiating the last line and remembering to use the quotient rule on the right side leads to y = 2yy (y2 + 1)2 Now since y = −1 y2 + 1 we can write the last equation as y = 2y (y2 + 1)2 y = 2y (y2 + 1)2 −1 (y2 + 1) = 2y −1 y2 + 1 3 = 2y(y )3 which is what we wanted to show. 34. Using implicit differentiation we get y = e−xy y = e−xy (−y − xy ) ye−xy + xe−xy y + y = 0 ye−xy + (xe−xy + 1)y = 0 Now since y = e−xy, substitute this into the last line to get yy + (xy + 1)y = 0 or (1 + xy)y + y2 = 0 which is what we wanted to show.
  • 37. Chapter 1 in Review 37 35. Substituting y = c1 + cos3x, y = −3 sin 3x, and y = −9 cos 3x into the left-hand side of the differential equation gives y + 9y = −9 cos 3x + 9 (c1 + cos 3x) = −9 cos 3x + 9c1 + 9 cos 3x = 9c1. Setting the result equal to the right-hand side of the differential equation yields 9c1 = 5 or c1 = 5 9 . Thus y = 5 9 + cos 3x. 36. Substituting y = c1 + c2x and y = c2 into the left-hand side of the differential equation gives y + 2y = c2 + 2 (c1 + c2x) = c2 + 2c1 + 2c2x. Setting the result equal to the right-hand side of the differential equation yields c2 + 2c1 + 2c2x = 3x (c2 + 2c1) + 2c2x = 0 + 3x. Therefore, 2c2 = 3 or c2 = 3 2 , and 3 2 + 2c1 = 0 or c1 = − 3 2 · 1 2 = −3 4. Thus y = − 3 4 + 3 2 x. 37. Differentiating and isolating c gives y = −ce−x + 4 ce−x = 4 − y c = ex 4 − y . Substituting the last equation into y = ce−x + 4x − 6 and simplifying yields y = ex 4 − y e−x + 4x − 6 y = 4 − y + 4x − 6 y + y = 4x − 2. Alternatively, isolating c in y = ce−x + 4x − 6 gives ce−x = y − 4x + 6 c = ex (y − 4x + 6) . Differentiating the last equation yields 0 = ex (y − 4x + 6) + ex y − 4 0 = ex y + y − 4x + 2 0 = y + y − 4x + 2 y + y = 4x − 2.
  • 38. 38 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS 38. The slope dy dx = 6 − 2x. Integrating with respect to x we have y = 6x − x2 + c1. Since the graph must pass through the origin (0, 0), we have y(0) = 6 (0) − (0)2 + c1 0 = c1. Thus y = 6x − x2. In Problems 39–42, y = c1e−3x + c2ex + 4x is given as a two-parameter family of solutions of the second-order differential equation y + 2y − 3y = −12x + 8. 39. If y(0) = 0 and y(0) = 0, then c1 + c2 = 0 −3c1 + c2 = −4 Subtracting the second equation from the first gives us 4c1 = 4 or c1 = 1, and thus c2 = −1. Therefore y = e−3x − ex + 4x. 40. If y(0) = 5 and y(0) = −11, then c1 + c2 = 5 −3c1 + c2 = −15 Subtracting the second equation from the first gives us 4c1 = 20 or c1 = 5, and thus c2 = 0. Therefore y = 5e−3x + 4x. 41. If y(1) = −2 and y(1) = 4, then c1e−3 + c2e = −6 −3c1e−3 + c2e = 0 Subtracting the second equation from the first gives us 4c1 = −6 or c1 = − 3 2 e3, and thus c2 = − 9 2 e−1. Therefore y = − 3 2 e−3x+3 − 9 2 ex−1 + 4x. 42. If y(−1) = 1 and y(−1) = 1, then c1e3 + c2e−1 = 5 −3c1e3 + c2e−1 = −3 Subtracting the second equation from the first gives us 4c1 = 8 or c1 = 2e−3, and thus c2 = 3e. Therefore y = 2e−3x−3 + 3ex+1 + 4x.
  • 39. Chapter 1 in Review 39 43. We are to use the Leibniz’s rule d dx v(x) u(x) F(x, t) dt = F(x, v(x)) dv dx − F(x, u(x)) du dx + v(x) u(x) ∂ ∂x F(x, t) dt Since y(x) = 1 3 x 0 f(t) sin (3x − 3t) dt, take F(x, t) = f(t) sin (3x − 3t), u(x) = 0, and v(x) = x to get d dx y(x) = d dx 1 3 x 0 f(t) sin (3x − 3t) dt = 1 3 F(x, x) · 1 − F(x, 0) · 0 + x 0 ∂ ∂x f(t) sin (3x − 3t) dt = 1 3 f(x) sin (3x − 3x) + x 0 3f(t) cos (3x − 3t) dt = x 0 f(t) cos (3x − 3t) dt Apply the Leibniz’s rule a second time to y = x 0 f(t) cos (3x − 3t) dt by taking F(x, t) = f(t) cos (3x − 3t), u(x) = 0, and v(x) = x to get d dx y (x) = d dx x 0 f(t) cos (3x − 3t) dt = F(x, x) · 1 − F(x, 0) · 0 + x 0 ∂ ∂x f(t) cos (3x − 3t) dt = f(x) cos (3x − 3x) + x 0 −3f(t) sin (3x − 3t) dt = f(x) − 3 x 0 f(t) sin (3x − 3t) dt Therefore y = f(x)−3 x 0 f(t) sin (3x − 3t) dt. Now substituting y and y into the differential equation we get y + 9y = f(x) − 3 x 0 f(t) sin (3x − 3t) dt + 9 · 1 3 x 0 f(t) sin (3x − 3t) dt = f(x)
  • 40. 40 CHAPTER 1 INTRODUCTION TO DIFFERENTIAL EQUATIONS 44. We are to use the Leibniz’s rule d dx v(x) u(x) F(x, t) dt = F(x, v(x)) dv dx − F(x, u(x)) du dx + v(x) u(x) ∂ ∂x F(x, t) dt Since y(x) = π 0 ex cos t dt, take F(x, t) = ex cos t, u(x) = 0, and v(x) = π to get d dx y(x) = d dx π 0 ex cos t dt = F(x, π) · 0 − F(x, 0) · 0 + π 0 ∂ ∂x ex cos t dt = π 0 cos tex cos t dt Applying the Leibniz’s rule a second time to y (x) = π 0 cos tex cos t dt by taking F(x, t) = cos tex cos t, u(x) = 0, and v(x) = π to get d dx y (x) = π 0 cos tex cos t dt = F(x, π) · 0 − F(x, 0) · 0 + π 0 ∂ ∂x cos tex cos t dt = π 0 cos2 tex cos t dt An alternative form of y can be obtained by integrating by parts with respect to t: y = sin t · ex cos t π 0 + π 0 xecos t sin2 t dt = π 0 xex cos t sin2 t dt = x π 0 ex cos t sin2 t dt. We use this last form of y instead of the first in the differential equation: xy + y = x π 0 ex cos t cos2 t dt + x π 0 ex cos t sin2 t dt = x π 0 ⎛ ⎜ ⎝ 1 cos2 t + sin2 t ⎞ ⎟ ⎠ ex cos t dt = y x π 0 ex cos t dt = xy xy + y = xy Therefore, xy + y − xy = 0.
  • 41. Chapter 1 in Review 41 45. From the graph we see that estimates for y0 and y1 are y0 = −3 and y1 = 0. 46. The differential equation is dh dt = − cA0 Aw 2gh . Using A0 = π(1/24)2 = π/576, Aw = π(2)2 = 4π, and g = 32, this becomes dh dt = − cπ/576 4π √ 64h = − c 288 √ h . 47. Let x(t) represent the height of the top of the rope at any time t with the positive direction upward as indicated. The weight of the portion of the rope off the ground is given by W = (x ft) · (1 lb/ft) = x. The mass of the rope is m = W/g = x/32. The net force is F = 5 − W = 5 − x. Now by Newton’s second law we get F = d dt (mv) = d dt x 32 · v = 5 − x Now expand the derivative and remember that v = dx/dt to get d dt x 32 · v = 5 − x 1 32 x dv dt + v dx dt = 5 − x x d dt dx dt + dx dt dx dt = 160 − 32x x d2x dt2 + dx dt 2 + 32x = 160
  翻译: