0% found this document useful (0 votes)
241 views52 pages

Differencial Equations Boyce Chapter 2 Solution

This document contains solutions to initial value problems for first order differential equations. Key points: - The integrating factor method is used to find general solutions to the differential equations. Initial conditions are then applied to determine specific solutions. - Solutions are analyzed to determine their long-term behavior as t approaches infinity. Critical values of initial conditions may cause a change between increasing/decreasing or diverging/converging behavior. - Specific solutions are found that satisfy given initial value problems and pass through requested points. Stationary points are located and classified for one solution.
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
241 views52 pages

Differencial Equations Boyce Chapter 2 Solution

This document contains solutions to initial value problems for first order differential equations. Key points: - The integrating factor method is used to find general solutions to the differential equations. Initial conditions are then applied to determine specific solutions. - Solutions are analyzed to determine their long-term behavior as t approaches infinity. Critical values of initial conditions may cause a change between increasing/decreasing or diverging/converging behavior. - Specific solutions are found that satisfy given initial value problems and pass through requested points. Stationary points are located and classified for one solution.
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 52

13

C HA P T E R
2
First Order Dierential Equations
2.1
5.(a)
(b) If y(0) > 2, solutions eventually have positive slopes, and hence increase with-
out bound. If y(0) 2, solutions have negative slopes and decrease without
bound.
(c) The integrating factor is (t) = e

3dt
= e
3t
. The dierential equation can be
written as e
3t
y

3e
3t
y = 4e
2t
, that is, (e
3t
y)

= 4e
2t
. Integration of both
sides of the equation results in the general solution y(t) = 2e
t
+c e
3t
. It follows
that all solutions will increase exponentially if c > 0 and will decrease exponentially
14 Chapter 2. First Order Dierential Equations
if c 0. Letting c = 0 and then t = 0, we see that the boundary of these behaviors
is at y(0) = 2.
9.(a)
(b) All solutions eventually have positive slopes, and hence increase without bound.
(c) The integrating factor is (t) = e

(1/2) dt
= e
t/2
. The dierential equation can
be written as e
t/2
y

+e
t/2
y/2 = 4t e
t/2
/2 , that is, (e
t/2
y/2)

= 2t e
t/2
. Integration
of both sides of the equation results in the general solution y(t) = 4t 8 +c e
t/2
.
All solutions approach the specic solution y
0
(t) = 4t 8 .
10.(a)
(b) For y > 0 , the slopes are all positive, and hence the corresponding solutions
increase without bound. For y < 0 , almost all solutions have negative slopes, and
hence solutions tend to decrease without bound.
(c) First divide both sides of the equation by t (t > 0). From the resulting standard
form, the integrating factor is (t) = e

(1/t) dt
= 1/t . The dierential equation
can be written as y

/t y/t
2
= t e
2t
, that is, ( y/t)

= t e
2t
. Integration leads to
the general solution y(t) = te
2t
/2 +c t . For c = 0 , solutions diverge, as implied
by the direction eld. For the case c = 0 , the specic solution is y(t) = te
2t
/2,
which evidently approaches zero as t .
2.1 15
12.(a)
(b) All solutions eventually have positive slopes, and hence increase without bound.
(c) The integrating factor is (t) = e
t/2
. The dierential equation can be written
as e
t/2
y

+e
t/2
y/2 = 4e
t/2
t
2
/2 , that is, (e
t/2
y/2)

= 2e
t/2
t
2
. Integration of both
sides of the equation results in the general solution y(t) = 4t
2
16t + 32 +c e
t/2
.
It follows that all solutions converge to the specic solution 4t
2
16t + 32 .
14. The integrating factor is (t) = e
3t
. After multiplying both sides by (t),
the equation can be written as (e
3t
y)

= t . Integrating both sides of the equation


results in the general solution y(t) = t
2
e
3t
/2 +c e
3t
. Invoking the specied con-
dition, we require that e
3
/2 +c e
3
= 0 . Hence c = 1/2 , and the solution to
the initial value problem is y(t) = (t
2
1)e
3t
/2 .
16. The integrating factor is (t) = e

(3/t) dt
= t
3
. Multiplying both sides by (t),
the equation can be written as (t
3
y)

= cos t . Integrating both sides of the equation


results in the general solution y(t) = sin t/t
3
+c t
3
. Substituting t = and setting
the value equal to zero gives c = 0 . Hence the specic solution is y(t) = sin t/t
3
.
17. The integrating factor is (t) = e
4t
, and the dierential equation can be
written as (e
4t
y)

= 1 . Integrating, we obtain e
4t
y(t) = t +c . Invoking the
specied initial condition results in the solution y(t) = (t + 2)e
4t
.
19. After writing the equation in standard form, we nd that the integrating
factor is (t) = e

(5/t) dt
= t
5
. Multiplying both sides by (t), the equation can be
written as (t
5
y)

= t e
t
. Integrating both sides results in t
5
y(t) = (t + 1)e
t
+
c . Letting t = 1 and setting the value equal to zero gives c = 0 . Hence the
specic solution of the initial value problem is y(t) = (t
4
+t
5
)e
t
.
16 Chapter 2. First Order Dierential Equations
22.(a)
The solutions eventually increase or decrease, depending on the initial value a .
The critical value seems to be a
0
= 2 .
(b) The integrating factor is (t) = e
t/2
, and the general solution of the dierential
equation is y(t) = 2e
t/4
+c e
t/4
. Invoking the initial condition y(0) = a , the
solution may also be expressed as y(t) = 2e
t/4
+ (a + 2) e
t/2
. The critical value
is a
0
= 2 .
(c) For a
0
= 2 , the solution is y(t) = 2e
t/4
, which diverges to as t .
23.(a)
Solutions appear to grow innitely large in absolute value, with signs depending on
the initial value y(0) = a
0
. The direction eld appears horizontal for a
0
1/8 .
(b) Dividing both sides of the given equation by 3, the integrating factor is (t) =
e
2t/3
. Multiplying both sides of the original dierential equation by (t) and inte-
grating results in y(t) = (2 e
2t/3
2 e
t/2
+a(4 + 3) e
2t/3
)/(4 + 3). The quali-
tative behavior of the solution is determined by the terms containing e
2t/3
: 2 e
2t/3
+
a(4 + 3) e
2t/3
. The nature of the solutions will change when 2 +a(4 + 3) = 0 .
Thus the critical initial value is a
0
= 2/(4 + 3) .
(c) In addition to the behavior described in part (a), when y(0) = 2/(4 + 3),
the solution is y(t) = (2 e
t/2
)/(4 + 3), and that specic solution will converge
to y = 0 .
2.1 17
24.(a)
As t 0 , solutions increase without bound if y(1) = a > 0.8 , and solutions de-
crease without bound if y(1) = a < 0.8 .
(b) The integrating factor is (t) = e

(t+1)/t dt
= t e
t
. The general solution of the
dierential equation is y(t) = 2t e
t
+c e
t
/t . Since y(1) = a, we have that 2 +
c = ae. That is, c = ae 2. Hence the solution can also be expressed as y(t) =
2t e
t
+ (ae 2) e
t
/t . For small values of t , the second term is dominant. Setting
ae 2 = 0 , the critical value of the parameter is a
0
= 2/e .
(c) When a = 2/e , the solution is y(t) = 2t e
t
, which approaches 0 as t 0 .
27. The integrating factor is (t) = e

(1/2) dt
= e
t/2
. Therefore the general solution
is y(t) = (4 cos t + 8 sin t)/5 +c e
t/2
. Invoking the initial condition, the specic so-
lution is y(t) = (4 cos t + 8 sin t 9 e
t/2
)/5 . Dierentiating, it follows that y

(t) =
(4 sin t + 8 cos t + 4.5 e
t/2
)/5 and y

(t) = (4 cos t 8 sin t 2.25 e


t/2
)/5. Set-
ting y

(t) = 0 , the rst solution is t


1
= 1.3643 , which gives the location of the rst
stationary point. Since y

(t
1
) < 0 , the rst stationary point in a local maximum.
The coordinates of the point are (1.3643 , 0.82008).
28. The integrating factor is (t) = e

(2/3) dt
= e
2t/3
, and the dierential equa-
tion can be written as (e
2t/3
y)

= e
2t/3
t e
2t/3
/2 . The general solution is y(t) =
(21 6t)/8 +c e
2t/3
. Imposing the initial condition, we have y(t) = (21 6t)/8 +
(y
0
21/8)e
2t/3
. Since the solution is smooth, the desired intersection will be
a point of tangency. Taking the derivative, y

(t) = 3/4 (2y


0
21/4)e
2t/3
/3.
Setting y

(t) = 0 , the solution is t


1
= (3/2) ln [(21 8y
0
)/9]. Substituting into the
solution, the respective value at the stationary point is y(t
1
) = 3/2 + (9/4) ln 3
(9/8) ln(21 8y
0
). Setting this result equal to zero, we obtain the required initial
value y
0
= (21 9 e
4/3
)/8 1.643 .
29.(a) The integrating factor is (t) = e
t/4
, and the dierential equation can be
written as (e
t/4
y)

= 3 e
t/4
+ 2 e
t/4
cos 2t. After integration, we get that the general
solution is y(t) = 12 + (8 cos 2t + 64 sin 2t)/65 +ce
t/4
. Invoking the initial condi-
tion, y(0) = 0 , the specic solution is y(t) = 12 + (8 cos 2t + 64 sin 2t 788 e
t/4
)/65.
As t , the exponential term will decay, and the solution will oscillate about
18 Chapter 2. First Order Dierential Equations
an average value of 12 , with an amplitude of 8/

65 .
(b) Solving y(t) = 12, we obtain the desired value t 10.0658.
31. The integrating factor is (t) = e
2t
, and the dierential equation can be
written as (e
2t
y)

= 3t e
2t
+ 2 e
t
. The general solution is y(t) = 3t/2 3/4
2 e
t
+c e
2t
. Imposing the initial condition, y(t) = 3t/2 3/4 2 e
t
+ (y
0
+ 11/4) e
2t
.
Now as t , the term containing e
2t
will dominate the solution. Its sign will
determine the divergence properties. Hence the critical value of the initial con-
dition is y
0
= 11/4. The corresponding solution, y(t) = 3t/2 3/4 2 e
t
, will
also decrease without bound.
Note on Problems 34-37 :
Let g(t) be given, and consider the function y(t) = y
1
(t) +g(t), in which y
1
(t) 0
as t . Dierentiating, y

(t) = y

1
(t) +g

(t) . Letting a be a constant, it follows


that y

(t) +ay(t) = y

1
(t) +ay
1
(t) +g

(t) +ag(t). Note that the hypothesis on the


function y
1
(t) will be satised, if y

1
(t) +ay
1
(t) = 0 . That is, y
1
(t) = c e
at
. Hence
y(t) = c e
at
+g(t), which is a solution of the equation y

+ay = g

(t) +ag(t). For


convenience, choose a = 1 .
34. Here g(t) = 5 , and we consider the linear equation y

+y = 5 . The integrating
factor is (t) = e
t
, and the dierential equation can be written as (e
t
y)

= 5e
t
. The
general solution is y(t) = 5 +c e
t
.
36. Here g(t) = 3t + 1. Consider the linear equation y

+y = 3 + 3t + 1. The
integrating factor is (t) = e
t
, and the dierential equation can be written as
(e
t
y)

= (3t + 4)e
t
. The general solution is y(t) = 3t + 1 +c e
t
.
37. g(t) = 1 t
2
. Consider the linear equation y

+y = 1 2t t
2
.The integrating
factor is (t) = e
t
, and the equation can be written as (e
t
y)

= (1 2t t
2
)e
t
. The
general solution is y(t) = 1 t
2
+c e
t
.
38.(a) Dierentiating y and using the fundamental theorem of calculus we obtain
that y

= Ae

p(t)dt
(p(t)), and then y

+p(t)y = 0.
(b) Dierentiating y we obtain that
y

= A

(t)e

p(t)dt
+A(t)e

p(t)dt
(p(t)).
If this satises the dierential equation then
y

+p(t)y = A

(t)e

p(t)dt
= g(t)
and the required condition follows.
(c) Let us denote (t) = e

p(t)dt
. Then clearly A(t) =
_
(t)g(t)dt, and after sub-
stitution y =
_
(t)g(t)dt (1/(t)), which is just Eq. (33).
2.2 19
40. We assume a solution of the form y = A(t)e

(1/t) dt
= A(t)e
ln t
= A(t)t
1
,
where A(t) satises A

(t) = 2t sin 2t. This implies that


A(t) = t cos 2t +
sin 2t
2
+c
and the solution is
y = cos 2t +
sin 2t
2t
+
c
t
.
41. First rewrite the dierential equation as
y

+
2
t
y =
2 sin t
t
.
Assume a solution of the form y = A(t)e

(2/t) dt
= A(t)t
2
, where A(t) satises
the ODE A

(t) = 2t sin t. It follows that A(t) = 2 sin t 2t cos t +c and thus


y = (2 sin t 2t cos t +c)/t
2
.
2.2
Problems 1 through 20 follow the pattern of the examples worked in this section.
The rst eight problems, however, do not have an initial condition, so the integra-
tion constant c cannot be found.
2. For x = 1 , the dierential equation may be written as y dy =
_
3x
2
/(1 +x
3
)

dx.
Integrating both sides, with respect to the appropriate variables, we obtain the re-
lation y
2
/2 = ln

1 +x
3

+c . That is, y(x) =


_
2 ln |1 +x
3
| +c .
3. The dierential equation may be written as y
2
dy = cos xdx. Integrating
both sides of the equation, with respect to the appropriate variables, we obtain the
relation y
1
= sin x +c . That is, (c + sin x)y = 1, in which c is an arbitrary
constant. Solving for the dependent variable, explicitly, y(x) = 1/(c + sin x) .
5. Write the dierential equation as cos
2
4y dy = cos
2
xdx, which also can be writ-
ten as sec
2
4y dy = cos
2
xdx. Integrating both sides of the equation, with respect to
the appropriate variables, we obtain the relation tan 4y = 2 sin x cos x + 2x +c .
7. The dierential equation may be written as (2y +e
y
)dy = (3x e
x
)dx. In-
tegrating both sides of the equation, with respect to the appropriate variables, we
obtain the relation y
2
+e
y
= 3x
2
/2 + e
x
+c .
8. Write the dierential equation as (2 +y
2
)dy = x
3
dx. Integrating both sides of
the equation, we obtain the relation 2y +y
3
/3 = x
4
/4 +c.
9.(a) The dierential equation is separable, with y
2
dy = (1 2x)dx. Integration
yields y
1
= x x
2
+c. Substituting x = 0 and y = 1/12, we nd that c = 12.
Hence the specic solution is y = 1/(x
2
x 12).
20 Chapter 2. First Order Dierential Equations
(b)
(c) Note that x
2
x 12 = (x + 3)(x 4) . Hence the solution becomes singular
at x = 3 and x = 4 , so the interval of existence is (3, 4).
11.(a) Rewrite the dierential equation as xe
x
dx = 2y dy . Integrating both sides
of the equation results in xe
x
e
x
= y
2
+c . Invoking the initial condition, we
obtain c = 3/4 . Hence y
2
= e
x
xe
x
3/4. The explicit form of the solution is
y(x) =
_
e
x
xe
x
3/4 . The positive sign is chosen, since y(0) = 1/2.
(b)
(c) The function under the radical becomes negative near x 0.96 and x 0.58.
12.(a) Write the dierential equation as r
2
dr =
1
d . Integrating both sides of
the equation results in the relation r
1
= ln +c . Imposing the condition r(1) =
3 , we obtain c = 1/3 . The explicit form of the solution is r = 3/(1 3 ln ).
2.2 21
(b)
(c) Clearly, the solution makes sense only if > 0 . Furthermore, the solution
becomes singular when ln = 1/3 , that is, =
3

e .
14.(a) Write the dierential equation as y
3
dy = x(1 +x
2
)
1/2
dx. Integrating
both sides of the equation, with respect to the appropriate variables, we obtain
the relation y
2
/2 =

1 +x
2
+c . Imposing the initial condition, we obtain
c = 3/2 . Hence the specic solution can be expressed as y
2
= 3 2

1 +x
2
.
The explicit form of the solution is y(x) = 1/
_
3 2

1 +x
2
. The positive sign is
chosen to satisfy the initial condition.
(b)
(c) The solution becomes singular when 2

1 +x
2
= 3 . That is, at x =

5 /2 .
16.(a) Rewrite the dierential equation as 4y
3
dy = x(x
2
+ 1)dx. Integrating both
sides of the equation results in y
4
= (x
2
+ 1)
2
/4 +c. Imposing the initial condition,
we obtain c = 0. Hence the solution may be expressed as (x
2
+ 1)
2
4y
4
= 0. The
explicit form of the solution is y(x) =
_
(x
2
+ 1)/2. The sign is chosen based on
y(0) = 1/

2.
22 Chapter 2. First Order Dierential Equations
(b)
(c) The solution is valid for all x R.
18.(a) Write the dierential equation as (3 + 4y)dy = (e
x
e
x
)dx. Integrating
both sides of the equation, with respect to the appropriate variables, we obtain the
relation 3y + 2y
2
= (e
x
+e
x
) +c . Imposing the initial condition, y(0) = 2 , we
obtain c = 16. Thus, the solution can be expressed as 3y + 2y
2
= (e
x
+e
x
) + 16.
Now by completing the square on the left hand side, 2(y + 3/4)
2
= (e
x
+e
x
) +
137/8. Hence the explicit form of the solution is y(x) = 3/4 +
_
137/16 cosh x.
(b)
(c) Note the 137 16 cosh x 0 as long as |x| > 2.8371 (approximately). Hence
the solution is valid on the interval 2.8371 < x < 2.8371.
20.(a) Rewrite the dierential equation as y
2
dy = arcsin x/

1 x
2
dx. Integrat-
ing both sides of the equation results in y
3
/3 = (arcsin x)
2
/2 +c. Imposing the
condition y(0) = 1, we obtain c = 1/3. The explicit form of the solution is y(x) =
(3(arcsin x)
2
/2 + 1)
1/3
.
2.2 23
(b)
(c) Since arcsin x is dened for 1 x 1, this is the interval of existence.
22. The dierential equation can be written as (3y
2
6)dy = 3x
2
dx. Integrating
both sides, we obtain y
3
6y = x
3
+c. Imposing the initial condition, the specic
solution is y
3
6y = x
3
1. Referring back to the dierential equation, we nd
that y

as y

2. The respective values of the abscissas are x 1.67,


1.88 . Hence the solution is valid for 1.67 < x < 1.88 .
24. Write the dierential equation as (3 + 2y)dy = (2 e
x
)dx. Integrating both
sides, we obtain 3y +y
2
= 2x e
x
+c. Based on the specied initial condition, the
solution can be written as 3y +y
2
= 2x e
x
+ 1. Completing the square, it follows
that y(x) = 3/2 +
_
2x e
x
+ 13/4. The solution is dened if 2x e
x
+ 13/4
0, that is, 1.5 x 2 (approximately). In that interval, y

= 0 for x = ln 2. It
can be veried that y

(ln 2) < 0. In fact, y

(x) < 0 on the interval of denition.


Hence the solution attains a global maximum at x = ln 2.
26. The dierential equation can be written as (1 +y
2
)
1
dy = 4(1 +x)dx. Inte-
grating both sides of the equation, we obtain arctan y = 4x + 2x
2
+c. Imposing
the given initial condition, the specic solution is arctan y = 4x + 2x
2
. Therefore,
y = tan(4x + 2x
2
). The solution is valid on the interval 0.537 < x < 0.336. Re-
ferring back to the dierential equation, the solution is stationary at x = 1. This
is not on the interval of existence, and there is no global minimum for the solution.
28.(a) Write the dierential equation as y
1
(4 y)
1
dy = t(1 +t)
1
dt . Integrat-
ing both sides of the equation, we obtain ln |y| ln |y 4| = 4t 4 ln |1 +t| +c .
Taking the exponential of both sides |y/(y 4)| = c e
4t
/(1 +t)
4
. It follows that as
t , |y/(y 4)| = |1 + 4/(y 4)| . That is, y(t) 4 .
(b) Setting y(0) = 2 , we obtain that c = 1. Based on the initial condition, the
solution may be expressed as y/(y 4) = e
4t
/(1 +t)
4
. Note that y/(y 4) < 0 ,
for all t 0. Hence y < 4 for all t 0. Referring back to the dierential equation,
it follows that y

is always positive. This means that the solution is monotone


increasing. We nd that the root of the equation e
4t
/(1 +t)
4
= 399 is near t =
2.844 .
24 Chapter 2. First Order Dierential Equations
(c) Note the y(t) = 4 is an equilibrium solution. Examining the local direction
eld we see that if y(0) > 0 , then the corresponding solutions converge to y =
4 . Referring back to part (a), we have y/(y 4) = [y
0
/(y
0
4)] e
4t
/(1 +t)
4
, for
y
0
= 4 . Setting t = 2 , we obtain y
0
/(y
0
4) = (3/e
2
)
4
y(2)/(y(2) 4). Now since
the function f(y) = y/(y 4) is monotone for y < 4 and y > 4 , we need only
solve the equations y
0
/(y
0
4) = 399(3/e
2
)
4
and y
0
/(y
0
4) = 401(3/e
2
)
4
. The
respective solutions are y
0
= 3.6622 and y
0
= 4.4042 .
32.(a) Observe that (2x
2
+ 3y
2
)/2xy = (y/x)
1
+ (3/2)(y/x). Hence the dieren-
tial equation is homogeneous.
(b) The substitution y = xv results in v +xv

= (2x
2
+ 3x
2
v
2
)/2x
2
v . The trans-
formed equation is v

= (2 +v
2
)/2xv . This equation is separable, with general
solution v
2
+ 2 = c x. In terms of the original dependent variable, the solution is
2x
2
+y
2
= c x
3
.
(c) The integral curves are symmetric with respect to the origin.
34.(a) Observe that (4x + 3y)/(2x +y) = 2 (y/x) [2 + (y/x)]
1
. Hence the
dierential equation is homogeneous.
(b) The substitution y = xv results in v +xv

= 2 v/(2 +v). The transformed


equation is v

= (v
2
+ 5v + 4)/(2 +v)x. This equation is separable, with general
solution (v + 4)
2
|v + 1| = c/x
3
. In terms of the original dependent variable, the
solution is (4x +y)
2
|x +y| = c.
2.2 25
(c) The integral curves are symmetric with respect to the origin.
36.(a) Divide by x
2
to see that the equation is homogeneous. Substituting y = xv ,
we obtain xv

= 1 + 2v + 2v
2
. The resulting dierential equation is separable.
(b) Write the equation as dv/(1 + 2v + 2v
2
) = x
1
dx. Integrating both sides of
the equation, we obtain the general solutionarctan(2v + 1) = ln |x| +c . In terms
of the original dependent variable, the solution is arctan((2y +x)/x) ln |x| = c.
(c) The integral curves are symmetric with respect to the origin.
37.(a) The dierential equation can be expressed as y

= (1/2)(y/x)
1
2(y/x).
Hence the equation is homogeneous. The substitution y = xv results in xv

=
(1 6v
2
)/2v. Separating variables, we have 2vdv/(1 6v
2
) = dx/x.
(b) Integrating both sides of the transformed equation yields (ln |1 6v
2
|)/6 =
ln |x| +c, that is, 1 6v
2
= c/ |x|
6
. In terms of the original dependent variable,
the general solution is 6y
2
= x
2
c/ |x|
4
.
26 Chapter 2. First Order Dierential Equations
(c) The integral curves are symmetric with respect to the origin.
38.(a) The dierential equation can be expressed as y

= (3/2)(y/x) (y/x)
1
.
Hence the equation is homogeneous. The substitution y = xv results in xv

=
(v
2
2)/2v, that is, 2vdv/(v
2
2) = dx/x.
(b) Integrating both sides of the transformed equation yields ln

v
2
2

= ln |x| +c,
that is, v
2
2 = c |x|. In terms of the original dependent variable, the general
solution is y
2
= c x
2
|x| + 2x
2
.
(c) The integral curves are symmetric with respect to the origin.
2.3
1. Let Q(t) be the amount of dye in the tank at time t. Clearly, Q(0) = 150 g.
The dierential equation governing the amount of dye is Q

(t) = 2Q(t)/150. The


solution of this separable equation is Q(t) = Q(0)e
t/75
= 150e
t/75
. We need the
time T such that Q(T) = 1.5 g. This means we have to solve 1.5 = 150e
T/75
and
we obtain that T = 75 ln(1/100) = 75 ln 100 345.4 min.
5.(a) Let Q be the amount of salt in the tank. Salt enters the tank of water at a
rate of 2 (1/4)(1 + (1/2) sin t) = 1/2 + (1/4) sin t dkg/min. It leaves the tank at a
2.3 27
rate of 2 Q/100 dkg/min. Hence the dierential equation governing the amount of
salt at any time is
dQ
dt
=
1
2
+
1
4
sin t
Q
50
.
The initial amount of salt is Q
0
= 50 dkg. The governing dierential equation
is linear, with integrating factor (t) = e
t/50
. Write the equation as (e
t/50
Q)

=
e
t/50
(1/2 + (1/4) sin t). The specic solution is Q(t) = 25 + (12.5 sin t 625 cos t +
63150 e
t/50
)/2501 dkg.
(b)
(c) The amount of salt approaches a steady state, which is an oscillation of approx-
imate amplitude 1/4 about a level of 25 dkg.
6.(a) Using the Principle of Conservation of Energy, the speed v of a particle falling
from a height h is given by
1
2
mv
2
= mgh.
(b) The outow rate is (outow cross-section area)(outow velocity): a

2gh.
At any instant, the volume of water in the tank is V (h) =
_
h
0
A(u)du. The time rate
of change of the volume is given by dV/dt = (dV/dh)(dh/dt) = A(h)dh/dt. Since
the volume is decreasing, dV/dt = a

2gh.
(c) With A(h) = 4, a = 0.01 , = 0.6 , the dierential equation for the wa-
ter level h is 4(dh/dt) = 0.006

2gh, with solution h(t) = 0.00000125gt


2

0.0015
_
2gh(0) t +h(0). Setting h(0) = 4 and g = 9.8 , h(t) = 0.000011025 t
2

0.0132816 t + 4, resulting in h(t) = 0 for t 602.3 s.


7.(a) The equation governing the value of the investment is dS/dt = r S . The value
of the investment, at any time, is given by S(t) = S
0
e
rt
. Setting S(T) = 2S
0
, the
required time is T = ln(2)/r .
(b) For the case r = .08 , T 8.66 yr.
28 Chapter 2. First Order Dierential Equations
(c) Referring to part (a), r = ln(2)/T . Setting T = 8 , the required interest rate is
to be approximately r = 8.66%.
12.(a) Using Eq.(15) we have dS/dt 0.005S = (800 + 10t), S(0) = 150, 000. Us-
ing an integrating factor and integration by parts we obtain that S(t) = 560, 000
410, 000e
0.005t
+ 2000t. Setting S(t) = 0 and solving numerically for t yields t =
146.54 months.
(b) The solution we obtained in part (a) with a general initial condition S(0) =
S
0
is S(t) = 560, 000 560, 000e
0.005t
+S
0
e
0.005t
+ 2000t. Solving the equation
S(240) = 0 yields S
0
= 246, 758.
13.(a) Let Q

= r Q. The general solution is Q(t) = Q


0
e
rt
. Based on the
denition of half-life, consider the equation Q
0
/2 = Q
0
e
5730 r
. It follows that
5730 r = ln(1/2), that is, r = 1.2097 10
4
per year.
(b) The amount of carbon-14 is given by Q(t) = Q
0
e
1.209710
4
t
.
(c) Given that Q(T) = Q
0
/4 , we have the equation 1/4 = e
1.209710
4
T
. Solving
for the decay time, the apparent age of the remains is approximately T = 11, 460
years.
15.(a) The dierential equation dy/dt = r(t) y k is linear, with integrating factor
(t) = e

r(t)dt
. Write the equation as (y)

= k (t) . Integration of both sides


yields the general solution y =
_
k
_
()d +y
0
(0)

/(t) . In this problem, the


integrating factor is (t) = e
(cos tt)/5
.
(b) The population becomes extinct, if y(t

) = 0 , for some t = t

. Referring to
part (a), we nd that y(t

) = 0 when
_
t

0
e
(cos )/5
d = 5 e
1/5
y
c
.
It can be shown that the integral on the left hand side increases monotonically, from
zero to a limiting value of approximately 5.0893 . Hence extinction can happen only
if 5 e
1/5
y
0
< 5.0893 . Solving 5e
1/5
y
c
= 5.0893 yields y
c
= 0.8333.
2.3 29
(c) Repeating the argument in part (b), it follows that y(t

) = 0 when
_
t

0
e
(cos )/5
d =
1
k
e
1/5
y
c
.
Hence extinction can happen only if e
1/5
y
0
/k < 5.0893 , so y
c
= 4.1667 k .
(d) Evidently, y
c
is a linear function of the parameter k .
17.(a) The solution of the governing equation satises u
3
= u
3
0
/( 3 u
3
0
t + 1 ). With
the given data, it follows that u(t) = 2000/
3
_
6 t/125 + 1 .
(b)
(c) Numerical evaluation results in u(t) = 600 for t 750.77 s.
22.(a) The dierential equation for the upward motion is mdv/dt = v
2
mg,
in which = 1/1325. This equation is separable, with m/(v
2
+mg) dv = dt .
Integrating both sides and invoking the initial condition, v(t) = 44.133 tan(0.425
0.222 t). Setting v(t
1
) = 0 , the ball reaches the maximum height at t
1
= 1.916 s .
Integrating v(t) , the position is given by x(t) = 198.75 ln [cos(0.222 t 0.425)] +
48.57 . Therefore the maximum height is x(t
1
) = 48.56 m.
(b) The dierential equation for the downward motion is mdv/dt = +v
2
mg .
This equation is also separable, with m/(mg v
2
) dv = dt . For convenience,
set t = 0 at the top of the trajectory. The new initial condition becomes v(0) = 0 .
Integrating both sides and invoking the initial condition, we obtain ln((44.13
v)/(44.13 +v)) = t/2.25. Solving for the velocity, v(t) = 44.13(1 e
t/2.25
)/(1 +
e
t/2.25
). Integrating v(t), we obtain x(t) = 99.29 ln(e
t/2.25
/(1 +e
t/2.25
)
2
) + 186.2.
To estimate the duration of the downward motion, set x(t
2
) = 0, resulting in
t
2
= 3.276 s. Hence the total time that the ball spends in the air is t
1
+t
2
= 5.192 s.
30 Chapter 2. First Order Dierential Equations
(c)
24.(a) Setting v
2
= v(dv/dx), we obtain dv/dx = v.
(b) The speed v of the sled satises ln(v/v
0
) = x. Noting that the unit conversion
factors cancel, solution of ln(24/240) = 600 results in = ln(10)/600 m
1

0.0038376 m
1
3.8376 km
1
.
(c) Solution of dv/dt = v
2
can be expressed as 1/v 1/v
0
= t. The elapsed
time is
t = (1/24 1/240)/(ln(10)/600)) 3600 35.18 s.
25.(a) Measure the positive direction of motion upward . The equation of motion
is given by mdv/dt = k v mg . The initial value problem is dv/dt = kv/m
g , with v(0) = v
0
. The solution is v(t) = mg/k + (v
0
+mg/k)e
kt/m
. Setting
v(t
m
) = 0, the maximum height is reached at time t
m
= (m/k) ln [(mg +k v
0
)/mg].
Integrating the velocity, the position of the body is
x(t) = mg t/k +
_
(
m
k
)
2
g +
mv
0
k
_
(1 e
kt/m
).
Hence the maximum height reached is
x
m
= x(t
m
) =
mv
0
k
g(
m
k
)
2
ln
_
mg +k v
0
mg
_
.
(b) Recall that for 1 , ln(1 +) =
2
/2 +
3
/3
4
/4 +. . ..
(c) The dimensions of the quantities involved are [k] = MT
1
, [v
0
] = LT
1
, [m] =
M and [g] = LT
2
. This implies that kv
0
/mg is dimensionless.
31.(a) Both equations are linear and separable. Initial conditions: v(0) = ucos A
and w(0) = usin A. We obtain the solutions v(t) = (ucos A)e
rt
and w(t) = g/r +
(usin A+g/r)e
rt
.
(b) Integrating the solutions in part (a), and invoking the initial conditions, the
coordinates are x(t) = ucos A(1 e
rt
)/r and
y(t) =
gt
r
+
g +ur sin A+hr
2
r
2
(
u
r
sin A+
g
r
2
)e
rt
.
2.3 31
(c)
(d) Let T be the time that it takes the ball to go 350 ft horizontally. Then from
above, e
T/5
= (u cos A70)/u cos A. At the same time, the height of the ball
is given by
y(T) = 160T + 803 + 5usin A
(800 + 5usin A)(ucos A70)
ucos A
.
Hence A and u must satisfy the equality
800 ln
_
u cos A70
u cos A
_
+ 803 + 5usin A
(800 + 5usin A)(ucos A70)
ucos A
= 10
for the ball to touch the top of the wall. To nd the optimal values for u and A,
consider u as a function of A and use implicit dierentiation in the above equation
to nd that
du
dA
=
u(u
2
cos A70u 11200 sin A)
11200 cos A
.
Solving this equation simultaneously with the above equation yields optimal values
for u and A: u 145.3 ft/s, A 0.644 rad.
32.(a) Solving equation (i), y

(x) =
_
(k
2
y)/y

1/2
. The positive answer is chosen,
since y is an increasing function of x.
(b) Let y = k
2
sin
2
t. Then dy = 2k
2
sin t cos tdt. Substituting into the equation in
part (a), we nd that
2k
2
sin t cos tdt
dx
=
cos t
sin t
.
Hence 2k
2
sin
2
tdt = dx.
(c) Setting = 2t, we further obtain k
2
sin
2
(/2) d = dx. Integrating both sides
of the equation and noting that t = = 0 corresponds to the origin, we obtain the
solutions x() = k
2
( sin )/2 and (from part (b)) y() = k
2
(1 cos )/2.
32 Chapter 2. First Order Dierential Equations
(d) Note that y/x = (1 cos )/( sin ). Setting x = 1 , y = 2 , the solution of
the equation (1 cos )/( sin ) = 2 is 1.401 . Substitution into either of
the expressions yields k 2.193 .
2.4
2. Rewrite the dierential equation as y

+ 1/(t(t 5)) y = 0. It is evident that the


coecient 1/t(t 5) is continuous everywhere except at t = 0 , 5 . Since the initial
condition is specied at t = 3 , Theorem 2.4.1 assures the existence of a unique
solution on the interval 0 < t < 5 .
3. The function tan t is discontinuous at odd multiples of /2. Since 3/2 < 2 <
5/2, the initial value problem has a unique solution on the interval (3/2, 5/2).
5. p(t) = 2t/(16 t
2
) and g(t) = 3t
2
/(16 t
2
). These functions are discontinuous
at x = 4 . The initial value problem has a unique solution on the interval (4 , 4).
6. The function ln t is dened and continuous on the interval (0 , ) . At t = 1,
ln t = 0, so the normal form of the dierential equation has a singularity there.
Also, cot t is not dened at integer multiples of , so the initial value problem will
have a solution on the interval (1, ).
7. The function f(t , y) is continuous everywhere on the plane, except along the
straight line y = 2t/5 . The partial derivative f/y = 16t/(2t + 5y)
2
has the
same region of continuity.
9. The function f(t , y) is discontinuous along the coordinate axes, and on the
hyperbola t
2
y
2
= 1 . Furthermore,
f
y
=
1
y(1 t
2
+y
2
)
2
y ln |ty|
(1 t
2
+y
2
)
2
has the same points of discontinuity.
10. f(t , y) is continuous everywhere on the plane. The partial derivative f/y is
also continuous everywhere.
12. The function f(t , y) is discontinuous along the lines t = k for k = 0, 1, 2, . . .
and y = 1 . The partial derivative f/y = cot t(2y +y
2
)/(1 +y)
2
has the same
region of continuity.
14. The equation is separable, with dy/y
2
= 4t dt . Integrating both sides, the
solution is given by y(t) = y
0
/(1 y
0
2t
2
). For y
0
> 0 , solutions exist as long as
t
2
< 1/2y
0
. For y
0
0 , solutions are dened for all t .
15. The equation is separable, with dy/y
5
= dt . Integrating both sides and
invoking the initial condition, y(t) = y
0
/
4
_
4y
4
0
t + 1 . Solutions exist as long as
2.4 33
4y
4
0
t + 1 > 0 , that is, 4y
4
0
t > 1 . If y
0
= 0 , solutions exist for t > 1/4y
4
0
. If
y
0
= 0 , then the solution y(t) = 0 exists for all t .
16. The function f(t , y) is discontinuous along the straight lines t = 1 and y = 0 .
The partial derivative f/y is discontinuous along the same lines. The equa-
tion is separable, with y dy = t
2
dt/(1 +t
3
). Integrating and invoking the initial
condition, the solution is y(t) =
_
(2/3) ln

1 +t
3

+y
2
0

1/2
. Solutions exist as long
as (2/3) ln

1 +t
3

+y
2
0
0, that is, y
2
0
(2/3) ln

1 +t
3

. For all y
0
(it can be
veried that y
0
= 0 yields a valid solution, even though Theorem 2.4.2 does not
guarantee one) , solutions exist as long as

1 +t
3

e
3y
2
0
/2
. From above, we must
have t > 1 . Hence the inequality may be written as t
3
e
3y
2
0
/2
1 . It follows
that the solutions are valid for (e
3y
2
0
/2
1)
1/3
< t < .
18.
Based on the direction eld, and the dierential equation, for y
0
< 0 , the slopes
eventually become negative, and hence solutions tend to . For y
0
> 0, solutions
increase without bound if t
0
< 0 . Otherwise, the slopes eventually become negative,
and solutions tend to zero. Furthermore, y
0
= 0 is an equilibrium solution. Note
that slopes are zero along the curves y = 0 and ty = 4 .
19.
For initial conditions (t
0
, y
0
) satisfying ty < 4 , the respective solutions all tend to
zero . For y
0
16, the solutions tend to 0; for y
0
> 16, the solutions tend to .
Also, y
0
= 0 is an equilibrium solution.
34 Chapter 2. First Order Dierential Equations
20.
Solutions with t
0
< 0 all tend to . Solutions with initial conditions (t
0
, y
0
)
to the right of the parabola t = 1 +y
2
asymptotically approach the parabola as
t . Integral curves with initial conditions above the parabola (and y
0
> 0)
also approach the curve. The slopes for solutions with initial conditions below the
parabola (and y
0
< 0) are all negative. These solutions tend to .
21.(a) No. There is no value of t
0
0 for which (2/3)(t t
0
)
2/3
satises the con-
dition y(1) = 1.
(b) Yes. Let t
0
= 1/2 in Eq.(19).
(c) For t
0
> 0, |y(2)| (4/3)
3/2
1.54.
24. The assumption is

(t) +p(t)(t) = 0. But then c

(t) +p(t)c(t) = 0 as well.


26.(a) Recalling Eq.(33) in Section 2.1,
y =
1
(t)
_
t
t0
(s)g(s) ds +
c
(t)
.
It is evident that y
1
(t) = 1/(t) and y
2
(t) = (1/(t))
_
t
t0
(s)g(s) ds.
(b) By denition, 1/(t) = e

p(t)dt
. Hence y

1
= p(t)/(t) = p(t)y
1
. That is,
y

1
+p(t)y
1
= 0.
(c) y

2
= (p(t)/(t))
_
t
0
(s)g(s) ds +(t)g(t)/(t) = p(t)y
2
+g(t). This implies
that y

2
+p(t)y
2
= g(t).
30. Since n = 3, set v = y
2
. It follows that v

= 2y
3
y

and y

= (y
3
/2)v

. Sub-
stitution into the dierential equation yields (y
3
/2)v

y = y
3
, which further
results in v

+ 2v = 2. The latter dierential equation is linear, and can be writ-


ten as (ve
2t
)

= 2e
2t
. The solution is given by v(t) = / +ce
2t
. Converting
back to the original dependent variable, y = v
1/2
= (/ +ce
2t
)
1/2
.
31. Since n = 3, set v = y
2
. It follows that v

= 2y
3
y

and y

= (y
3
/2)v

. The
dierential equation is written as (y
3
/2)v

(cos t +T)y = y
3
, which upon
2.5 35
further substitution is v

+ 2(cos t +T)v = 2. This ODE is linear, with integrat-


ing factor (t) = e
2

(cos t+T)dt
= e
2sin t+2Tt
. The solution is
v(t) = 2e
(2sin t+2Tt)
_
t
0
e
2sin +2T
d +ce
(2sin t+2Tt)
.
Converting back to the original dependent variable, y = v
1/2
.
33. The solution of the initial value problem y

1
+ 4y
1
= 0, y
1
(0) = 1 is y
1
(t) = e
4t
.
Therefore y(1

) = y
1
(1) = e
4
. On the interval (1, ), the dierential equation
is y

2
+y
2
= 0, with y
2
(t) = ce
t
. Therefore y(1
+
) = y
2
(1) = ce
1
. Equating the
limits y(1

) = y(1
+
), we require that c = e
3
. Hence the global solution of the
initial value problem is
y(t) =
_
e
4t
, 0 t 1
e
3t
, t > 1
.
Note the discontinuity of the derivative
y

(t) =
_
4e
4t
, 0 < t < 1
e
3t
, t > 1
.
2.5
1.
36 Chapter 2. First Order Dierential Equations
0
For y
0
0 , the only equilibrium point is y

= 0, and y

= ay +by
2
> 0 when y > 0,
hence the equilibrium solution y = 0 is unstable.
2.
0
-a/b
The equilibrium points are y

= a/b and y

= 0, and y

> 0 when y > 0 or y <


a/b, and y

< 0 when a/b < y < 0, therefore the equilibrium solution y = a/b
is asymptotically stable and the equilibrium solution y = 0 is unstable.
2.5 37
4.
0
The only equilibrium point is y

= 0, and y

> 0 when y > 0, y

< 0 when y < 0,


hence the equilibrium solution y = 0 is unstable.
6.
38 Chapter 2. First Order Dierential Equations
0
The only equilibrium point is y

= 0, and y

> 0 when y < 0, y

< 0 when y > 0,


hence the equilibrium solution y = 0 is asymptotically stable.
8.
2
The only equilibrium point is y

= 2, and y

< 0 for y = 2. As long as y


0
= 2,
the corresponding solution is monotone decreasing. Hence the equilibrium solution
y = 2 is semistable.
2.5 39
10.
0
2
-2
The equilibrium points are y

= 0, 2 , and y

> 0 for y < 2 or 0 < y < 2 and


y

< 0 for 2 < y < 0 or y > 2. The equilibrium solution y = 0 is unstable, and the
remaining two are asymptotically stable.
12.
40 Chapter 2. First Order Dierential Equations
0
3
-3
The equilibrium points are y

= 0 , 3, and y

< 0 when y < 3 or y > 3, and


y

> 0 for 3 < y < 0 or 0 < y < 3. The equilibrium solutions y = 3 and y = 3
are unstable and asymptotically stable, respectively. The equilibrium solution y = 0
is semistable.
13.
0
3
The equilibrium points are y

= 0, 3. y

> 0 for all y except y = 0 and y = 3. Both


equilibrium solutions are semistable.
15.(a) Inverting Eq.(11), Eq.(13) shows t as a function of the population y and the
carrying capacity K. With y
0
= K/3,
t =
1
r
ln

(1/3) [1 (y/K)]
(y/K) [1 (1/3)]

.
2.5 41
Setting y = 2y
0
,
=
1
r
ln

(1/3) [1 (2/3)]
(2/3) [1 (1/3)]

.
That is, = (ln 4)/r. If r = 0.03 per year, 46.21 years.
(b) In Eq.(13), set y
0
/K = and y/K = . As a result, we obtain
T =
1
r
ln

[1 ]
[1 ]

.
Given = 0.1, = 0.9 and r = 0.03 per year, 146.48 years.
19.(a) The rate of increase of the volume is given by rate of ow inrate of ow out.
That is, dV/dt = k a

2gh. Since the cross section is constant, dV/dt = Adh/dt.


Hence the governing equation is dh/dt = (k a

2gh)/A.
(b) Setting dh/dt = 0, the equilibrium height is h
e
= (1/2g)(k/a)
2
. Furthermore,
since dh/dt < 0 for h > h
e
and dh/dt > 0 for h < h
e
, it follows that the equilibrium
height is asymptotically stable.
22.(a) The equilibrium points are at y

= 0 and y

= 1. Since f

(y) = 2y ,
the equilibrium solution y = 0 is unstable and the equilibrium solution y = 1 is
asymptotically stable.
(b) The dierential equation is separable, with [y(1 y)]
1
dy = dt . Integrating
both sides and invoking the initial condition, the solution is
y(t) =
y
0
e
t
1 y
0
+y
0
e
t
=
y
0
y
0
+ (1 y
0
)e
t
.
It is evident that (independent of y
0
) lim
t
y(t) = 0 and lim
t
y(t) = 1 .
23.(a) y(t) = y
0
e
t
.
(b) From part (a), dx/dt = xy
0
e
t
. Separating variables, dx/x = y
0
e
t
dt.
Integrating both sides, the solution is x(t) = x
0
e
y0(1e
t
)/
.
(c) As t , y(t) 0 and x(t) x
0
e
y0/
. Over a long period of time,
the proportion of carriers vanishes. Therefore the proportion of the population that
escapes the epidemic is the proportion of susceptibles left at that time, x
0
e
y0/
.
26.(a) For a < 0 , the only critical point is at y = 0 , which is asymptotically stable.
For a = 0 , the only critical point is at y = 0 , which is asymptotically stable. For
a > 0 , the three critical points are at y = 0 ,

a . The critical point at y = 0 is


unstable, whereas the other two are asymptotically stable.
(b) Below, we graph solutions in the case a = 1, a = 0 and a = 1 respectively.
42 Chapter 2. First Order Dierential Equations
2.5 43
(c)
27.(a) f(y) = y(a y); f

(y) = a 2y. For a < 0, the critical points are at y = a
and y = 0. Observe that f

(a) > 0 and f

(0) < 0 . Hence y = a is unstable and
y = 0 asymptotically stable. For a = 0 , the only critical point is at y = 0 , which is
semistable since f(y) = y
2
is concave down. For a > 0 , the critical points are at
y = 0 and y = a . Observe that f

(0) > 0 and f

(a) < 0 . Hence y = 0 is unstable
and y = a asymptotically stable.
(b) Below, we graph solutions in the case a = 1, a = 0 and a = 1 respectively.
44 Chapter 2. First Order Dierential Equations
(c)
2.6
1. M(x, y) = 4x + 3 and N(x, y) = 6y 1 . Since M
y
= N
x
= 0 , the equation is
exact. Integrating M with respect to x, while holding y constant, yields (x, y) =
2x
2
+ 3x +h(y) . Now
y
= h

(y) , and equating with N results in the possible


function h(y) = 3y
2
y . Hence (x, y) = 2x
2
+ 3x + 3y
2
y , and the solution is
dened implicitly as 2x
2
+ 3x + 3y
2
y = c .
2. M(x, y) = 3x y and N(x, y) = x 3y . Note that M
y
= N
x
, and hence the
dierential equation is not exact.
4. First divide both sides by (4xy + 4). We now have M(x, y) = y and N(x, y) = x.
Since M
y
= N
x
= 0 , the resulting equation is exact. Integrating M with respect
to x, while holding y constant, results in (x, y) = xy +h(y) . Dierentiating with
respect to y ,
y
= x +h

(y) . Setting
y
= N, we nd that h

(y) = 0 , and hence


h(y) = 0 is acceptable. Therefore the solution is dened implicitly as xy = c . Note
that if xy + 1 = 0 , the equation is trivially satised.
6. Write the equation as (ax by)dx + (bx cy)dy = 0. Now M(x, y) = ax by
and N(x, y) = bx cy. Since M
y
= N
x
, the dierential equation is not exact.
2.6 45
8. M(x, y) = e
x
sin y + 2y and N(x, y) = 2x +e
x
sin y . Note that M
y
= N
x
, and
hence the dierential equation is not exact.
10. M(x, y) = y/x + 4x and N(x, y) = ln x 3. Since M
y
= N
x
= 1/x, the given
equation is exact. Integrating N with respect to y , while holding x constant,
results in (x, y) = y ln x 3y +h(x) . Dierentiating with respect to x,
x
=
y/x +h

(x). Setting
x
= M, we nd that h

(x) = 4x, and hence h(x) = 2x


2
.
Therefore the solution is dened implicitly as 2x
2
+y ln x 3y = c .
11. M(x, y) = x ln y +xy and N(x, y) = y ln x +xy. Note that M
y
= N
x
, and
hence the dierential equation is not exact.
13. M(x, y) = 2x y and N(x, y) = 2y x. Since M
y
= N
x
= 1, the equation is
exact. Integrating M with respect to x, while holding y constant, yields (x, y) =
x
2
xy +h(y). Now
y
= x +h

(y). Equating
y
with N results in h

(y) =
2y, and hence h(y) = y
2
. Thus (x, y) = x
2
xy +y
2
, and the solution is given
implicitly as x
2
xy +y
2
= c . Invoking the initial condition y(1) = 4 , the specic
solution is x
2
xy +y
2
= 13. The explicit form of the solution is y(x) = (x +

52 3x
2
)/2. Hence the solution is valid as long as 3x
2
52 .
16. M(x, y) = y e
2xy
+ 3x and N(x, y) = bxe
2xy
. Note that M
y
= e
2xy
+ 2xy e
2xy
,
and N
x
= b e
2xy
+ 2bxy e
2xy
. The given equation is exact, as long as b = 1 . In-
tegrating N with respect to y , while holding x constant, results in (x, y) =
e
2xy
/2 +h(x) . Now dierentiating with respect to x,
x
= y e
2xy
+h

(x). Set-
ting
x
= M, we nd that h

(x) = 3x, and hence h(x) = 3x


2
/2 . We conclude that
(x, y) = e
2xy
/2 + 3x
2
/2 . Hence the solution is given implicitly as e
2xy
+ 3x
2
= c .
17. Note that is of the form (x, y) = f(x) +g(y), since each of the integrands
is a function of a single variable. It follows that
x
= f

(x) and
y
= g

(y). That
is,
x
= M(x, y
0
) and
y
= N(x
0
, y). Furthermore,

xy
(x
0
, y
0
) =
M
y
(x
0
, y
0
) and

2

yx
(x
0
, y
0
) =
N
x
(x
0
, y
0
) ,
based on the hypothesis and the fact that the point (x
0
, y
0
) is arbitrary,
xy
=
yx
and M
y
(x, y) = N
x
(x, y).
18. Observe that (M(x))
y
= (N(y))
x
= 0 .
20. M
y
= y
1
cos y y
2
sin y and N
x
= 3 e
x
(cos x + sin x)/y . Multiplying
both sides by the integrating factor (x, y) = y e
x
, the given equation can be written
as (e
x
sin y 3y sin x)dx + (e
x
cos y + 3 cos x)dy = 0 . Let

M = M and

N = N.
Observe that

M
y
=

N
x
, and hence the latter ODE is exact. Integrating

N with
respect to y , while holding x constant, results in (x, y) = e
x
sin y + 3y cos x +
h(x) . Now dierentiating with respect to x,
x
= e
x
sin y 3y sin x +h

(x). Set-
ting
x
=

M, we nd that h

(x) = 0 , and hence h(x) = 0 is feasible. Hence the


solution of the given equation is dened implicitly by e
x
sin y + 3y cos x = c.
46 Chapter 2. First Order Dierential Equations
21. M
y
= 1 and N
x
= 2 . Multiply both sides by the integrating factor (x, y) = y
to obtain y
2
dx + (2xy 3y
2
e
y
)dy = 0. Let

M = yM and

N = yN. It is easy to
see that

M
y
=

N
x
, and hence the latter ODE is exact. Integrating

M with respect
to x yields (x, y) = xy
2
+h(y) . Equating
y
with

N results in h

(y) = 3y
2
e
y
,
and hence h(y) = 3e
y
(y
2
2y + 2). Thus (x, y) = xy
2
3e
y
(y
2
2y + 2), and
the solution is dened implicitly by xy
2
3e
y
(y
2
2y + 2) = c .
24. The equation M +Ny

= 0 has an integrating factor if (M)


y
= (N)
x
,
that is,
y
M
x
N = N
x
M
y
. Suppose that N
x
M
y
= R(xM yN), in
which R is some function depending only on the quantity z = xy . It follows that
the modied form of the equation is exact, if
y
M
x
N = R(xM yN) =
R(xM yN). This relation is satised if
y
= (x)R and
x
= (y)R. Now
consider = (xy). Then the partial derivatives are
x
=

y and
y
=

x. Note
that

= d/dz . Thus must satisfy

(z) = R(z). The latter equation is sepa-


rable, with d = R(z)dz , and (z) =
_
R(z)dz . Therefore, given R = R(xy), it is
possible to determine = (xy) which becomes an integrating factor of the dier-
ential equation.
28. The equation is not exact, since N
x
M
y
= 2y 1 . However, (N
x
M
y
)/M =
(2y 1)/y is a function of y alone. Hence there exists = (y) , which is a solution
of the dierential equation

= (2 1/y). The latter equation is separable, with


d/ = 2 1/y . One solution is (y) = e
2yln y
= e
2y
/y . Now rewrite the given
ODE as e
2y
dx + (2xe
2y
1)dy = 0 . This equation is exact, and it is easy to see
that (x, y) = xe
2y
y. Therefore the solution of the given equation is dened
implicitly by xe
2y
y = c .
30. The given equation is not exact, since N
x
M
y
= 8x
3
/y
3
+ 6/y
2
. But note that
(N
x
M
y
)/M = 2/y is a function of y alone, and hence there is an integrating fac-
tor = (y). Solving the equation

= (2/y), an integrating factor is (y) = y


2
.
Now rewrite the dierential equation as (4x
3
+ 3y)dx + (3x + 2y
3
)dy = 0. By in-
spection, (x, y) = x
4
+ 3xy +y
4
/2, and the solution of the given equation is de-
ned implicitly by x
4
+ 3xy +y
4
/2 = c .
32. Multiplying both sides of the ODE by = [xy(2x +y)]
1
, the given equation is
equivalent to
_
(3x +y)/(2x
2
+xy)

dx +
_
(x +y)/(2xy +y
2
)

dy = 0 . Rewrite the
dierential equation as
_
2
x
+
2
2x +y
_
dx +
_
1
y
+
1
2x +y
_
dy = 0 .
It is easy to see that M
y
= N
x
. Integrating M with respect to x, while keep-
ing y constant, results in (x, y) = 2 ln |x| + ln |2x +y| +h(y) . Now taking the
partial derivative with respect to y ,
y
= (2x +y)
1
+h

(y) . Setting
y
= N,
we nd that h

(y) = 1/y , and hence h(y) = ln |y| . Therefore (x, y) = 2 ln |x| +


ln |2x +y| + ln |y|, and the solution of the given equation is dened implicitly by
2x
3
y +x
2
y
2
= c .
2.7 47
2.7
2. The Euler formula is given by y
n+1
= y
n
+h(3y
n
1) = (1 + 3h)y
n
h.
(a) 1.2, 1.46, 1.798, 2.2374
(b) 1.215, 1.49934, 1.87537, 2.37268
(c) 1.22365, 1.52232, 1.92119, 2.45386
(d) The dierential equation is linear with solution y(t) = (1 + 2e
3t
)/3. The values
are 1.23324, 1.54808, 1.97307, 2.54674.
5.
All solutions seem to converge to y = 25/4 .
7.
All solutions seem to converge to a specic function.
48 Chapter 2. First Order Dierential Equations
8.
Solutions with initial conditions |y(0)| > 2.5 seem to diverge. On the other hand,
solutions with initial conditions |y(0)| < 2.5 seem to converge to zero. Also, y = 0
is an equilibrium solution.
10.
Solutions with positive initial conditions increase without bound. Solutions with
negative initial conditions decrease without bound. Note that y = 0 is an equilib-
rium solution.
11. The Euler formula is y
n+1
= y
n
2h

y
n
+ 5h. The initial value is y
0
= 2.
(a) 2.95148, 3.65042, 4.18248, 4.59561, 4.92048, 5.17816
(b) 2.93808, 3.63167, 4.16190, 4.57499, 4.90074, 5.15979
(c) 2.93159, 3.62253, 4.15183, 4.56486, 4.89102, 5.15072
(d) 2.92776, 3.73397, 4.14586, 4.55884, 4.88523, 5.14531
12. The Euler formula is y
n+1
= (1 + 4h)y
n
ht
n
y
2
n
. The initial value is (t
0
, y
0
) =
(0, 0.5).
(a) 2.42831, 4.49378, 3.13540, 2.26560, 1.76850, 1.44951
2.7 49
(b) 2.65748, 4.46004, 3.14322, 2.27392, 1.77288, 1.45195
(c) 2.78710, 4.44811, 3.14607, 2.27804, 1.77516, 1.45322
(d) 2.86988, 4.44232, 3.14756, 2.28049, 1.77655, 1.45399
14. The Euler formula is y
n+1
= (1 ht
n
)y
n
+hy
3
n
/10, with (t
0
, y
0
) = (0, 1).
(a) 0.950517, 0.687550, 0.369188, 0.145990, 0.0421429, 0.00872877
(b) 0.938298, 0.672145, 0.362640, 0.147659, 0.0454100, 0.0104931
(c) 0.932253, 0.664778, 0.359567, 0.148416, 0.0469514, 0.0113722
(d) 0.928649, 0.660463, 0.357783, 0.148848, 0.0478492, 0.0118978
17. The Euler formula is y
n+1
= y
n
+h(y
2
n
+ 2t
n
y
n
)/(3 +t
2
n
). The initial point is
(t
0
, y
0
) = (1 , 2). Using this iteration formula with the specied h values, the value
of the solution at t = 2.5 is somewhere between 18 and 19. At t = 3 there is no
reliable estimate.
19.(a)
(b) The iteration formula is y
n+1
= y
n
+hy
2
n
ht
2
n
. The critical value
0
appears
to be between 0.67 and 0.68. For y
0
>
0
, the iterations diverge.
20.(a) The ODE is linear, with general solution y(t) = t +ce
t
. Invoking the spec-
ied initial condition, y(t
0
) = y
0
, we have y
0
= t
0
+ce
t0
. Hence c = (y
0
t
0
)e
t0
.
Thus the solution is given by (t) = (y
0
t
0
)e
tt0
+t.
(b) The Euler formula is y
n+1
= (1 +h)y
n
+h ht
n
. Now set k = n + 1 .
(c) We have y
1
= (1 +h)y
0
+h ht
0
= (1 +h)y
0
+ (t
1
t
0
) ht
0
. Rearranging
the terms, y
1
= (1 +h)(y
0
t
0
) +t
1
. Now suppose that y
k
= (1 +h)
k
(y
0
t
0
) +
t
k
, for some k 1. Then y
k+1
= (1 +h)y
k
+h ht
k
. Substituting for y
k
, we nd
50 Chapter 2. First Order Dierential Equations
that
y
k+1
= (1 +h)
k+1
(y
0
t
0
) + (1 +h)t
k
+h ht
k
= (1 +h)
k+1
(y
0
t
0
) +t
k
+h.
Noting that t
k+1
= t
k
+h, the result is veried.
(d) Substituting h = (t t
0
)/n, with t
n
= t , y
n
= (1 + (t t
0
)/n)
n
(y
0
t
0
) + t.
Taking the limit of both sides, and using the fact that lim
n
(1 +a/n)
n
= e
a
,
pointwise convergence is proved.
21. The exact solution is y(t) = e
2t
. The Euler formula is y
n+1
= (1 + 2h)y
n
. It is
easy to see that y
n
= (1 + 2h)
n
y
0
= (1 + 2h)
n
. Given t > 0 , set h = t/n. Taking
the limit, we nd that lim
n
y
n
= lim
n
(1 + 2t/n)
n
= e
2t
.
23. The exact solution is y(t) = t/2 +e
2t
. The Euler formula is y
n+1
= (1 +
2h)y
n
+h/2 ht
n
. Since y
0
= 1 , y
1
= (1 + 2h) +h/2 = (1 + 2h) +t
1
/2 . It is
easy to show by mathematical induction, that y
n
= (1 + 2h)
n
+t
n
/2 . For t >
0 , set h = t/n and thus t
n
= t . Taking the limit, we nd that lim
n
y
n
=
lim
n
[(1 + 2t/n)
n
+t/2] = e
2t
+t/2. Hence pointwise convergence is proved.
2.8
2. Let z = y 3 and = t + 1 . It follows that dz/d = (dz/dt)(dt/d) = dz/dt .
Furthermore, dz/dt = dy/dt = 4 y
3
. Hence dz/d = 4 (z + 3)
3
. The new ini-
tial condition is z(0) = 0 .
3.(a) The approximating functions are dened recursively by

n+1
(t) =
_
t
0
3 [
n
(s) + 1] ds .
Setting
0
(t) = 0 ,
1
(t) = 3t . Continuing,
2
(t) = 9t
2
/2 + 3t ,
3
(t) = 9t
3
/2 +
9t
2
/2 + 3t ,
4
(t) = 9t
4
/4 + 9t
3
/2 + 9t
2
/2 + 3t , . . . . Based upon these we conjec-
ture that
n
(t) =

n
k=1
3
k
t
k
/k! and use mathematical induction to verify this form
for
n
(t). First, let n = 1, then
n
(t) = 3t, so it is certainly true for n = 1. Then,
using Eq.(7) again we have

n+1
(t) =
_
t
0
3 [
n
(s) + 1] ds =
_
t
0
3
_
n

k=1
3
k
k !
s
k
+ 1
_
ds =
n+1

k=1
3
k
k !
t
k
,
and we have veried our conjecture.
2.8 51
(b)
(c) Recall from calculus that e
at
= 1 +

k=1
a
k
t
k
/k!. Thus
(t) =

k=1
3
k
k !
t
k
= e
3t
1 .
(d)
From the plot it appears that
4
is a good estimate for |t| < 1/2.
5.(a) The approximating functions are dened recursively by

n+1
(t) =
_
t
0
[
n
(s)/3 +s] ds .
Setting
0
(t) = 0,
1
(t) = t
2
/2. Continuing,
2
(t) = t
2
/2 t
3
/18,
3
(t) = t
2
/2
t
3
/18 +t
4
/216,
4
(t) = t
2
/2 t
3
/18 +t
4
/216 t
5
/3240, . . .. Based upon these we
conjecture that
n
(t) =

n
k=1
9(1/3)
k+1
t
k+1
/(k + 1)! and use mathematical in-
duction to verify this form for
n
(t).
52 Chapter 2. First Order Dierential Equations
(b)
(c) Recall from calculus that e
at
= 1 +

k=1
a
k
t
k
/k!. Thus
(t) =

k=1
9
(1/3)
k+1
k + 1 !
t
k+1
= 9e
t/3
+ 3t 9 .
(d)
From the plot it appears that
4
is a good estimate for |t| < 2.
6.(a) The approximating functions are dened recursively by

n+1
(t) =
_
t
0
[
n
(s) + 1 s] ds .
Setting
0
(t) = 0,
1
(t) = t t
2
/2,
2
(t) = t t
3
/6,
3
(t) = t t
4
/24,
4
(t) = t
t
5
/120, . . . . Based upon these we conjecture that
n
(t) = t t
n+1
/(n + 1)! and
use mathematical induction to verify this form for
n
(t).
2.8 53
(b)
(c) Clearly (t) = t.
(d)
From the plot it appears that
4
is a good estimate for |t| < 1.
8.(a) The approximating functions are dened recursively by

n+1
(t) =
_
t
0
_
s
2

n
(s) s

ds .
Set
0
(t) = 0. The iterates are given by
1
(t) = t
2
/2 ,
2
(t) = t
2
/2 t
5
/10 ,

3
(t) = t
2
/2 t
5
/10 t
8
/80 ,
4
(t) = t
2
/2 t
5
/10 t
8
/80 t
11
/880 ,. . . . Upon
inspection, it becomes apparent that

n
(t) = t
2
_
1
2
+
t
3
2 5
+
t
6
2 5 8
+. . . +
(t
3
)
n1
2 5 8 . . . [2 + 3(n 1)]
_
=
= t
2
n

k=1
(t
3
)
k1
2 5 8 . . . [2 + 3(k 1)]
.
54 Chapter 2. First Order Dierential Equations
(b)
(c) Using the identity
n
(t) =
1
(t) + [
2
(t)
1
(t)] + [
3
(t)
2
(t)] +. . . + [
n
(t)

n1
(t)], consider the series
1
(t) +

k=1
[
k+1
(t)
k
(t)]. Fix any t value now.
We use the Ratio Test to prove the convergence of this series:

k+1
(t)
k
(t)

k
(t)
k1
(t)

(t
2
)(t
3
)
k
25(2+3k)
(t
2
)(t
3
)
k1
25(2+3(k1))

=
|t|
3
2 + 3k
.
The limit of this quantity is 0 for any xed t as k , and we obtain that
n
(t)
is convergent for any t.
9.(a) The approximating functions are dened recursively by

n+1
(t) =
_
t
0
_
2s
2
+
2
n
(s)

ds .
Set
0
(t) = 0. The rst three iterates are given by
1
(t) = 2t
3
/3,
2
(t) = 2t
3
/3 +
4t
7
/63,
3
(t) = 2t
3
/3 + 4t
7
/63 + 16t
11
/2079 + 16t
15
/59535 .
(b)
The iterates appear to be converging.
2.8 55
12.(a) The approximating functions are dened recursively by

n+1
(t) =
_
t
0
_
3s
2
+ 4s + 2
2(
n
(s) 1)
_
ds .
Note that 1/(2y 2) = (1/2)

6
k=0
y
k
+O(y
7
). For computational purposes, use
the geometric series sum to replace the above iteration formula by

n+1
(t) =
1
2
_
t
0
_
(3s
2
+ 4s + 2)
6

k=0

k
n
(s)
_
ds .
Set
0
(t) = 0. The rst four approximations are given by
1
(t) = t t
2
t
3
/2,

2
(t) = t t
2
/2 +t
3
/6 +t
4
/4 t
5
/5 t
6
/24 +. . .,
3
(t) = t t
2
/2 +t
4
/12
3t
5
/20 + 4t
6
/45 +. . .,
4
(t) = t t
2
/2 +t
4
/8 7t
5
/60 +t
6
/15 +. . .
(b)
The approximations appear to be converging to the exact solution, which can be
found by separating the variables: (t) = 1

1 + 2t + 2t
2
+t
3
.
14.(a)
n
(0) = 0, for every n 1 . Let a (0 , 1]. Then
n
(a) = 2na e
na
2
=
2na/e
na
2
. Using lHospitals rule, lim
z
2az/e
az
2
= lim
z
1/ze
az
2
= 0. Hence
lim
n

n
(a) = 0 .
(b)
_
1
0
2nxe
nx
2
dx = e
nx
2

1
0
= 1 e
n
. Therefore,
lim
n
_
1
0

n
(x)dx =
_
1
0
lim
n

n
(x)dx.
15. Let t be xed, such that (t , y
1
), (t , y
2
) D. Without loss of generality, assume
that y
1
< y
2
. Since f is dierentiable with respect to y, the mean value theorem as-
serts that there exists (y
1
, y
2
) such that f(t , y
1
) f(t , y
2
) = f
y
(t , )(y
1
y
2
).
This means that |f(t , y
1
) f(t , y
2
)| = |f
y
(t , )| |y
1
y
2
|. Since, by assumption,
f/y is continuous in D, f
y
attains a maximum K on any closed and bounded
subset of D. Hence |f(t , y
1
) f(t , y
2
)| K |y
1
y
2
|.
16. For a suciently small interval of t,
n1
(t),
n
(t) D. Since f satises a
Lipschitz condition, |f(t,
n
(t)) f(t,
n1
(t))| K|
n
(t)
n1
(t)|. Here K =
max |f
y
|.
56 Chapter 2. First Order Dierential Equations
17.(a)
1
(t) =
_
t
0
f(s , 0)ds . Hence |
1
(t)|
_
|t|
0
|f(s , 0)| ds
_
|t|
0
Mds = M|t| , in
which M is the maximum value of |f(t , y)| on D.
(b) By denition,
2
(t)
1
(t) =
_
t
0
[f(s ,
1
(s)) f(s , 0)] ds . Taking the absolute
value of both sides, |
2
(t)
1
(t)|
_
|t|
0
|[f(s ,
1
(s)) f(s , 0)]| ds . Based on the
results in Problems 16 and 17,
|
2
(t)
1
(t)|
_
|t|
0
K|
1
(s) 0| ds KM
_
|t|
0
|s| ds .
Evaluating the last integral, we obtain that |
2
(t)
1
(t)| MK|t|
2
/2 .
(c) Suppose that
|
i
(t)
i1
(t)|
MK
i1
|t|
i
i!
for some i 1 . By denition,

i+1
(t)
i
(t) =
_
t
0
[f(s,
i
(s)) f(s,
i1
(s))] ds .
It follows that
|
i+1
(t)
i
(t)|
_
|t|
0
|f(s,
i
(s)) f(s,
i1
(s))| ds

_
|t|
0
K|
i
(s)
i1
(s)| ds
_
|t|
0
K
MK
i1
|s|
i
i!
ds =
=
MK
i
|t|
i+1
(i + 1)!

MK
i
h
i+1
(i + 1)!
.
Hence, by mathematical induction, the assertion is true.
18.(a) Use the triangle inequality, |a +b| |a| +|b| .
(b) For |t| h, |
1
(t)| Mh, and |
n
(t)
n1
(t)| MK
n1
h
n
/(n!) . Hence
|
n
(t)| M
n

i=1
K
i1
h
i
i !
=
M
K
n

i=1
(Kh)
i
i !
.
(c) The sequence of partial sums in (b) converges to M(e
Kh
1)/K. By the com-
parison test, the sums in (a) also converge. Since individual terms of a convergent
series must tend to zero, |
n
(t)
n1
(t)| 0 , and it follows that the sequence
|
n
(t)| is convergent.
19.(a) Let (t) =
_
t
0
f(s , (s))ds and (t) =
_
t
0
f(s , (s))ds . Then by linearity of
the integral, (t) (t) =
_
t
0
[f(s , (s)) f(s , (s))] ds .
(b) It follows that |(t) (t)|
_
t
0
|f(s , (s)) f(s , (s))| ds .
2.9 57
(c) We know that f satises a Lipschitz condition, |f(t , y
1
) f(t , y
2
)| K |y
1
y
2
|,
based on |f/y| K in D. Therefore,
|(t) (t)|
_
t
0
|f(s , (s)) f(s , (s))| ds
_
t
0
K|(s) (s)| ds.
2.9
1. Writing the equation for each n 0 , y
1
= 0.8 y
0
, y
2
= 0.8 y
1
= (0.8)
2
y
0
,
y
3
= 0.8 y
2
= (0.8)
3
y
0
and so on, it is apparent that y
n
= (0.8)
n
y
0
. The terms
constitute an alternating series, which converge to zero, regardless of y
0
.
3. Write the equation for each n 0, y
1
=

2 y
0
, y
2
=
_
3/2 y
1
, y
3
=
_
4/3 y
2
, . . .
Upon substitution, we nd that y
2
=
_
(2 3)/2 y
0
=

3y
0
, y
3
=
_
(4 3 2)/(3 2) y
0
=

4y
0
, . . . It can be proved by mathematical induction, that y
n
=

n + 1y
0
. This
sequence is divergent, except for y
0
= 0 .
4. Writing the equation for each n 0 , y
1
= y
0
, y
2
= y
1
, y
3
= y
2
, y
4
= y
3
,
and so on. It can be shown that
y
n
=
_
y
0
, for n = 4k or n = 4k 1
y
0
, for n = 4k 2 or n = 4k 3
The sequence is convergent only for y
0
= 0 .
6. Writing the equation for each n 0 ,
y
1
= 0.5 y
0
+ 4
y
2
= 0.5 y
1
+ 4 = 0.5(0.5 y
0
+ 4) + 4 = (0.5)
2
y
0
+ 4 + (0.5)4
y
3
= 0.5 y
2
+ 4 = 0.5(0.5 y
1
+ 4) + 4 = (0.5)
3
y
0
+ 4
_
1 + (0.5) + (0.5)
2

.
.
.
y
n
= (0.5)
n
y
0
+ (8/3) [1 (0.5)
n
]
which follows from Eq.(13) and (14). The sequence is convergent for all y
0
, and in
fact y
n
8/3.
8. Let y
n
be the balance at the end of the nth month. Then y
n+1
= (1 +r/12)y
n
+
25. We have y
n
=
n
[y
0
25/(1 )] + 25/(1 ), in which = (1 +r/12). Here
r is the annual interest rate, given as 7%. Thus y
48
= (1.00583333)
48
[1000 + 12 25/r]
12 25/r = $2, 702.28.
9. Let y
n
be the balance due at the end of the nth month. The appropriate
dierence equation is y
n+1
= (1 +r/12) y
n
P . Here r is the annual interest rate
and P is the monthly payment. The solution, in terms of the amount borrowed,
is given by y
n
=
n
[y
0
+P/(1 )] P/(1 ), in which = (1 +r/12) and y
0
=
9, 000 . To gure out the monthly payment P, we require that y
36
= 0. That is,
58 Chapter 2. First Order Dierential Equations

36
[y
0
+P/(1 )] = P/(1 ). After the specied amounts are substituted, we
nd that P = $282.03.
11. Let y
n
be the balance due at the end of the nth month. The appropriate dier-
ence equation is y
n+1
= (1 +r/12) y
n
P, in which r = .06 and P is the monthly
payment. The initial value of the mortgage is y
0
= $100, 000. Then the balance
due at the end of the n-th month is y
n
=
n
[y
0
+P/(1 )] P/(1 ), where
= (1 +r/12). In terms of the specied values, y
n
= (1.005)
n
[10
5
12P/r] +
12P/r. Setting n = 30 12 = 360 , and y
360
= 0 , we nd that P = $599.55. For the
monthly payment corresponding to a 20 year mortgage, set n = 240 and y
240
= 0
to nd that P = $716.43. The total amount paid during the term of the loan is
360 599.55 = $215, 838.00 for the 30-year loan and is 240 716.43 = $171, 943.20
for the 20-year loan.
12. Let y
n
be the balance due at the end of the nth month, with y
0
the initial value
of the mortgage. The appropriate dierence equation is y
n+1
= (1 +r/12) y
n
P,
in which r = 0.08 and P = $1000 is the maximum monthly payment. Given that the
life of the mortgage is 20 years, we require that y
240
= 0. The balance due at the end
of the n-th month is y
n
=
n
[y
0
+P/(1 )] P/(1 ). In terms of the specied
values for the parameters, the solution of (1.00666667)
240
[y
0
12 1000/0.08] =
12 1000/0.08 is y
0
= $119, 554.29.
19.(a)
2
= (
2

1
)/(
3

2
) = (3.449 3)/(3.544 3.449) = 4.7263 .
(b) di= (|
2
|/) 100 = (|4.6692 4.7363|/4.6692) 100 1.22%.
(c) Assuming (
3

2
)/(
4

3
) = ,
4
3.5643
(d) A period 16 solution appears near 3.565 .
(e) Note that (
n+1

n
) =
1
n
(
n

n1
). With the assumption that
n
= , we
have (
n+1

n
) =
1
(
n

n1
), which is of the form y
n+1
= y
n
, n 3 . It
2.9 59
follows that (
k

k1
) =
3k
(
3

2
) for k 4 . Then

n
=
1
+ (
2

1
) + (
3

2
) + (
4

3
) +. . . + (
n

n1
)
=
1
+ (
2

1
) + (
3

2
)
_
1 +
1
+
2
+. . . +
3n

=
1
+ (
2

1
) + (
3

2
)
_
1
4n
1
1
_
.
Hence lim
n

n
=
2
+ (
3

2
)
_

1
_
. Substitution of the appropriate values
yields
lim
n

n
= 3.5699
PROBLEMS
1. The equation is linear. It can be written in the form y

+ 2y/x = x
3
, and the
integrating factor is (x) = e

(2/x) dx
= e
2 ln x
= x
2
. Multiplication by (x) yields
x
2
y

+ 2yx = (yx
2
)

= x
5
. Integration with respect to x and division by x
2
gives
that y = x
4
/6 +c/x
2
.
5. The equation is exact. Algebraic manipulations give the symmetric form of
the equation, (2xy +y
2
+ 3)dx + (x
2
+ 2xy)dy = 0. We can check that M
y
= 2x +
2y = N
x
, so the equation is really exact. Integrating M with respect to x gives that
(x, y) = x
2
y +xy
2
+ 3x +g(y), then
y
= x
2
+ 2xy +g

(y) = x
2
+ 2xy, so we get
that g

(y) = 0, so we obtain that g(y) = 0 is acceptable. Therefore the solution is


dened implicitly as x
2
y +xy
2
+ 3x = c.
6. The equation is linear. It can be written in the form y

+ (1 + (1/x))y = 1/x
and the integrating factor is (x) = e

1+(1/x) dx
= e
x+ln x
= xe
x
. Multiplication by
(x) yields xe
x
y

+ (xe
x
+e
x
)y = (xe
x
y)

= e
x
. Integration with respect to x and
division by xe
x
shows that the general solution of the equation is y = 1/x +c/(xe
x
).
The initial condition implies that 0 = 1 +c/e, which means that c = e and the
solution is y = 1/x e/(xe
x
) = x
1
(1 e
1x
).
7. The equation is separable. Separation of variables gives the dierential equation
y(2 + 3y)dy = (5x
4
+ 1)dx, and then after integration we obtain that the solution
is x
5
+x y
2
y
3
= c.
8. The equation is linear. It can be written in the form y

+ 2y/x = cos x/x


2
and the
integrating factor is (x) = e

(2/x) dx
= e
2 ln x
= x
2
. Multiplication by (x) gives
x
2
y

+ 2xy = (x
2
y)

= cos x, and after integration with respect to x and division by


x
2
we obtain the general solution y = (c + sin x)/x
2
. The initial condition implies
that c = 4 sin 2 and the solution becomes y = (4 sin 2 + sin x)/x
2
.
11. The equation is exact. It is easy to check that M
y
= 1 = N
x
. Integrating
M with respect to x gives that (x, y) = x
4
/4 +xy +g(y), then
y
= x +g

(y) =
60 Chapter 2. First Order Dierential Equations
x +e
y
, which means that g

(y) = e
y
, so we obtain that g(y) = e
y
. Therefore the
solution is dened implicitly as x
4
/4 +xy +e
y
= c.
13. The equation is separable. Factoring the right hand side leads to the equa-
tion y

= (1 +y
2
)(1 + 3x
2
). We separate the variables to obtain dy/(1 +y
2
) =
(1 + 3x
2
)dx, then integration gives us arctan y = x +x
3
+c. The solution is y =
tan(x +x
3
+c).
14. The equation is exact. We can check that M
y
= 1 = N
x
. Integrating M with
respect to x gives that (x, y) = x
2
/2 +xy +g(y), then
y
= x +g

(y) = x + 2y,
which means that g

(y) = 2y, so we obtain that g(y) = y


2
. Therefore the general
solution is dened implicitly as x
2
/2 +xy +y
2
= c. The initial condition gives us
c = 10, so the solution is x
2
+ 2xy + 2y
2
= 20.
15. The equation is separable. Separation of variables leads us to the equation
dy
y
=
1 e
x
1 +e
x
dx.
Note that 1 +e
x
2e
x
= 1 e
x
. We obtain that
ln |y| =
_
1 e
x
1 +e
x
dx =
_
1
2e
x
1 +e
x
dx = x 2 ln(1 +e
x
) + c.
This means that y = ce
x
(1 +e
x
)
2
, which also can be written as y = c/ cosh
2
(x/2)
after some algebraic manipulations.
16. The equation is exact. The symmetric form is (e
x
cos y +e
2y
cos x)dx +
(e
x
sin y + 2e
2y
sin x)dy = 0. We can check that M
y
= e
x
sin y + 2e
2y
cos x =
N
x
. Integrating M with respect to x gives that (x, y) = e
x
cos y +e
2y
sin x +
g(y), then
y
= e
x
sin y + 2e
2y
sin x +g

(y) = e
x
sin y + 2e
2y
sin x, so we get
that g

(y) = 0, so we obtain that g(y) = 0 is acceptable. Therefore the solution is


dened implicitly as e
x
cos y +e
2y
sin x = c.
17. The equation is linear. The integrating factor is (x) = e

3 dx
= e
3x
, which
turns the equation into e
3x
y

3e
3x
y = (e
3x
y)

= e
x
. We integrate with respect
to x to obtain e
3x
y = e
x
+c, and the solution is y = ce
3x
+e
4x
after multiplication
by e
3x
.
18. The equation is linear. The integrating factor is (x) = e

2 dx
= e
2x
, which
gives us e
2x
y

+ 2e
2x
y = (e
2x
y)

= e
x
2
. The antiderivative of the function on the
right hand side can not be expressed in a closed form using elementary functions,
so we have to express the solution using integrals. Let us integrate both sides of
this equation from 0 to x. We obtain that the left hand side turns into
_
x
0
(e
2s
y(s))

ds = e
2x
y(x) e
0
y(0) = e
2x
y 4.
The right hand side gives us
_
x
0
e
s
2
ds. So we found that
y = e
2x
_
x
0
e
s
2
ds + 4e
2x
.
2.9 61
19. The equation is exact. Algebraic manipulations give us the symmetric form
(y
3
+ 2y 4x
3
)dx + (2x + 3xy
2
)dy = 0. We can check that M
y
= 3y
2
+ 2 = N
x
.
Integrating M with respect to x gives that (x, y) = xy
3
+ 2xy x
4
+g(y), then

y
= 3xy
2
+ 2x +g

(y) = 2x + 3xy
2
, which means that g

(y) = 0, so we obtain that


g(y) = 0 is acceptable. Therefore the solution is xy
3
+ 2xy x
4
= c.
20. The equation is separable, because y

= e
x+y
= e
x
e
y
. Separation of variables
yields the equation e
y
dy = e
x
dx, which turns into e
y
= e
x
+c after integration
and we obtain the implicitly dened solution e
x
+e
y
= c.
22. The equation is separable. Separation of variables turns the equation into
(y
2
+ 2)dy = (x
2
1)dx, which, after integration, gives y
3
/3 + 2y = x
3
/3 x +c.
The initial condition yields c = 5/3, and the solution is y
3
+ 6y x
3
+ 3x = 5.
23. The equation is linear. Division by t gives y

+ (1 + (1/t))y = e
2t
/t, so the
integrating factor is (t) = e

(1+(1/t))dt
= e
t+ln t
= te
t
. The equation turns into
te
t
y

+ (te
t
+e
t
)y = (te
t
y)

= e
3t
. Integration therefore leads to te
t
y = e
3t
/3 +c
and the solution is y = e
2t
/(3t) +ce
t
/t.
24. The equation is exact. We can check that M
y
= 2 cos y sin xcos x = N
x
. In-
tegrating M with respect to x gives that (x, y) = sin y sin
2
x +g(y), then
y
=
cos y sin
2
x +g

(y) = cos y sin


2
x, which means that g

(y) = 0, so we obtain that


g(y) = 0 is acceptable. Therefore the solution is dened implicitly as sin y sin
2
x = c.
25. The equation is exact. We can check that
M
y
=
2x
y
2

x
2
y
2
(x
2
+y
2
)
2
= N
x
.
Integrating M with respect to x gives that (x, y) = x
2
/y + arctan(y/x) +g(y),
then
y
= x
2
/y
2
+x/(x
2
+y
2
) +g

(y) = x/(x
2
+y
2
) x
2
/y
2
, which means that
g

(y) = 0, so we obtain that g(y) = 0 is acceptable. Therefore the solution is dened


implicitly as x
2
/y + arctan(y/x) = c.
28. The equation can be made exact by choosing an appropriate integrating factor.
We can check that (M
y
N
x
)/N = (2 1)/x = 1/x depends only on x, so (x) =
e

(1/x)dx
= e
ln x
= x is an integrating factor. After multiplication, the equation
becomes (2yx + 3x
2
)dx +x
2
dy = 0. This equation is exact now, because M
y
=
2x = N
x
. Integrating M with respect to x gives that (x, y) = yx
2
+x
3
+g(y),
then
y
= x
2
+g

(y) = x
2
, which means that g

(y) = 0, so we obtain that g(y) = 0


is acceptable. Therefore the solution is dened implicitly as x
3
+x
2
y = c.
29. The equation is homogeneous. (See Section 2.2, Problem 30) We can see that
y

=
x +y
x y
=
1 (y/x)
1 + (y/x)
.
We substitute u = y/x, which means also that y = ux and then y

= u

x +u =
62 Chapter 2. First Order Dierential Equations
(1 u)/(1 +u), which implies that
u

x =
1 u
1 +u
u =
1 2u u
2
1 +u
,
a separable equation. Separating the variables yields
1 +u
1 2u u
2
du =
dx
x
,
and then integration gives ln(1 2u u
2
)/2 = ln |x| +c. Substituting u = y/x
back into this expression and using that
ln(1 2(y/x) (y/x)
2
)/2 ln |x| = ln(|x|
_
1 2y/x (y/x)
2
) = ln(
_
x
2
2yx y
2
)
we obtain that the solution is x
2
2yx y
2
= c.
30. The equation is homogeneous. (See Section 2.2, Problem 30) Algebraic manip-
ulations show that it can be written in the form
y

=
3y
2
+ 2xy
2xy +x
2
=
3(y/x)
2
+ 2(y/x)
2(y/x) + 1
.
Substituting u = y/x gives that y = ux and then
y

= u

x +u =
3u
2
+ 2u
2u + 1
,
which implies that
u

x =
3u
2
+ 2u
2u + 1
u =
u
2
+u
2u + 1
,
a separable equation. We obtain that (2u + 1)du/(u
2
+u) = dx/x, which in turn
means that ln(u
2
+u) = ln |x| + c. Therefore, u
2
+u = cx and then substituting
u = y/x gives us the solution (y
2
/x
3
) + (y/x
2
) = c.
31. The equation can be made exact by choosing an appropriate integrating
factor. We can check that (M
y
N
x
)/M = (3x
2
+y)/(y(3x
2
+y)) = 1/y de-
pends only on y, so (y) = e

(1/y)dy
= e
ln y
= y is an integrating factor. After
multiplication, the equation becomes (3x
2
y
2
+y
3
)dx + (2x
3
y + 3xy
2
)dy = 0. This
equation is exact now, because M
y
= 6x
2
y + 3y
2
= N
x
. Integrating M with re-
spect to x gives that (x, y) = x
3
y
2
+y
3
x +g(y), then
y
= 2x
3
y + 3y
2
x +g

(y) =
2x
3
y + 3xy
2
, which means that g

(y) = 0, so we obtain that g(y) = 0 is acceptable.


Therefore the general solution is dened implicitly as x
3
y
2
+xy
3
= c. The initial
condition gives us 4 8 = c = 4, and the solution is x
3
y
2
+xy
3
= 4.
33. Let y
1
be a solution, i.e. y

1
= q
1
+q
2
y
1
+q
3
y
2
1
. Now let y = y
1
+ (1/v) also be
a solution. Dierentiating this expression with respect to t and using that y is also
a solution we obtain y

= y

1
(1/v
2
)v

= q
1
+q
2
y +q
3
y
2
= q
1
+q
2
(y
1
+ (1/v)) +
q
3
(y
1
+ (1/v))
2
. Now using that y
1
was also a solution we get that (1/v
2
)v

=
q
2
(1/v) + 2q
3
(y
1
/v) +q
3
(1/v
2
), which, after some simple algebraic manipulations
turns into v

= (q
2
+ 2q
3
y
1
)v q
3
.
2.9 63
35.(a) The equation is y

= (1 y)(x +by) = x + (b x)y by


2
. We set y = 1 +
(1/v) and dierentiate: y

= v
2
v

= x + (b x)(1 + (1/v)) b(1 + (1/v))


2
, which,
after simplication, turns into v

= (b +x)v +b.
(b) When x = at, the equation is v

(b +at)v = b, so the integrating factor is


(t) = e
btat
2
/2
. This turns the equation into (v(t))

= b(t), so v(t) =
_
b(t)dt,
and then v = (b
_
(t)dt)/(t).
36. Substitute v = y

, then v

= y

. The equation turns into t


2
v

+ 2tv = (t
2
v)

=
2, which yields t
2
v = 2t +c
1
, so y

= v = (2/t) + (c
1
/t
2
). Integrating this expres-
sion gives us the solution y = 2 ln t (c
1
/t) +c
2
.
37. Set v = y

, then v

= y

. The equation with this substitution is tv

+v =
(tv)

= 4, which gives tv = 4t +c
1
, so y

= v = 4 + (c
1
/t). Integrating this expres-
sion yields the solution y = 4t +c
1
ln t +c
2
.
38. Set v = y

, so v

= y

. The equation is v

+ 2tv
2
= 0, which is a separable
equation. Separating the variables we obtain dv/v
2
= 2tdt, so 1/v = t
2
+c,
and then y

= v = 1/(t
2
+c
1
). Now depending on the value of c
1
, we have the
following possibilities: when c
1
= 0, then y = 1/t +c
2
, when 0 < c
1
= k
2
, then
y = (1/k) arctan(t/k) +c
2
, and when 0 > c
1
= k
2
then
y = (1/2k) ln |(t k)/(t +k)| +c
2
.
We also divided by v = y

when we separated the variables, and v = 0 (which is


y = c) is also a solution.
39. Substitute v = y

and v

= y

. The equation is 2t
2
v

+v
3
= 2tv. This is a
Bernoulli equation (See Section 2.4, Problem 27), so the substitution z = v
2
yields
z

= 2v
3
v

, and the equation turns into 2t


2
v

v
3
+ 1 = 2t/v
2
, i.e. into 2t
2
z

/2 +
1 = 2tz, which in turn simplies to t
2
z

+ 2tz = (t
2
z)

= 1. Integration yields t
2
z =
t +c, which means that z = (1/t) + (c/t
2
). Now y

= v =
_
1/z = t/

t +c
1
and another integration gives
y =
2
3
(t 2c
1
)

t +c
1
+c
2
.
The substitution also loses the solution v = 0, i.e. y = c.
40. Set v = y

, then v

= y

. The equation reads v

+v = 2e
t
, which is a linear
equation with integrating factor (t) = e
t
. This turns the equation into e
t
v

+e
t
v =
(e
t
v)

= 2, which means that e


t
v = 2t +c and then y

= v = 2te
t
+ce
t
. Another
integration yields the solution y = 2te
t
+c
1
e
t
+c
2
.
41. Let v = y

and v

= y

. The equation is t
2
v

= v
2
, which is a separable equation.
Separating the variables we obtain dv/v
2
= dt/t
2
, which gives us 1/v = (1/t) +
c
1
, and then y

= v = t/(1 +c
1
t). Now when c
1
= 0, then y = t
2
/2 +c
2
, and when
c
1
= 0, then y = t/c
1
(ln |1 +c
1
t|)/c
2
1
+c
2
. Also, at the separation we divided by
v = 0, which also gives us the solution y = c.
64 Chapter 2. First Order Dierential Equations
43. Set y

= v(y). Then y

= v

(y)(dy/dt) = v

(y)v(y). We obtain the equation


v

v +y = 0, where the dierentiation is with respect to y. This is a separable


equation which simplies to vdv = ydy. We obtain that v
2
/2 = y
2
/2 +c, so
y

= v(y) =
_
c y
2
. We separate the variables again to get dy/
_
c y
2
= dt,
so arcsin(y/

c) = t +d, which means that y =

c sin(t +d) = c
1
sin(t +c
2
).
44. Set y

= v(y). Then y

= v

(y)(dy/dt) = v

(y)v(y). We obtain the equation


v

v +yv
3
= 0, where the dierentiation is with respect to y. Separation of variables
turns this into dv/v
2
= ydy, which gives us y

= v = 2/(c
1
+y
2
). This implies
that (c
1
+y
2
)dy = 2dt and then the solution is dened implicitly as c
1
y +y
3
/3 =
2t +c
2
. Also, y = c is a solution which we lost when divided by y

= v = 0.
46. Set y

= v(y). Then y

= v

(y)(dy/dt) = v

(y)v(y). We obtain the equa-


tion yv

v v
3
= 0, where the dierentiation is with respect to y. This separa-
ble equation gives us dv/v
2
= dy/y, which means that 1/v = ln |y| +c, and then
y

= v = 1/(c ln |y|). We separate variables again to obtain (c ln |y|)dy = dt,


and then integration yields the implicitly dened solution cy (y ln |y| y) = t +d.
Also, y = c is a solution which we lost when we divided by v = 0.
49. Set y

= v(y). Then y

= v

(y)(dy/dt) = v

(y)v(y). We obtain the equation


v

v 3y
2
= 0, where the dierentiation is with respect to y. Separation of variables
gives vdv = 3y
2
dy, and after integration this turns into v
2
/2 = y
3
+c. The initial
conditions imply that c = 0 here, so (y

)
2
= v
2
= 2y
3
. This implies that y

2y
3/2
(the sign is determined by the initial conditions again), and this separable equation
now turns into y
3/2
dy =

2dt. Integration yields 2y


1/2
=

2t +d, and the


initial conditions at this point give that d =

2. Algebraic manipulations nd
that y = 2(1 t)
2
.
50. Set v = y

, then v

= y

. The equation with this substitution turns into


the equation (1 +t
2
)v

+ 2tv = ((1 +t
2
)v)

= 3t
2
. Integrating this we get that
(1 +t
2
)v = 3t
1
+c, and c = 5 from the initial conditions. This means that
y

= v = 3/(t(1 +t
2
)) 5/(1 +t
2
). The partial fraction decomposition of the rst
expression shows that y

= 3/t 3t/(1 +t
2
) 5/(1 +t
2
) and then another inte-
gration here gives us that y = 3 ln t (3/2) ln(1 +t
2
) 5 arctan t +d. The initial
conditions identify d = 2 + (3/2) ln 2 + 5/4, and we obtained the solution.

You might also like

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy