Math 115 HW #9 Solutions

1. Solve the differential equation

dy
sin x + (cos x)y = x sin(x2 ).
dx

Answer: Dividing everything by sin x yields the linear first order equation

cos x x sin(x2 )
y0 + y= .
sin x sin x
cos x x sin(x2 )
Here P (x) = sin x and Q(x) = sin x , so
Z Z
cos x
P (x)dx = dx = ln | sin x|
sin x
(doing a u-substitution with u = sin x). Thus,
R
µ(x) = e P (x)dx
= eln |sinx| = | sin x|.

We can drop the absolute value signs (they would cancel in the final expression for y, anyway)
and let µ(x) = sin x. Then

x sin(x2 )
Z Z Z
1 1
µ(x)Q(x)dx = sin x dx = x sin(x2 )dx = (− cos(x2 ) + C) = − cos(x2 ) + C 0
sin x 2 2

(by a u-substitution with u = x2 ). Therefore,

1 cos(x2 ) C0
Z  
1 1 1 2 0
y= µ(x)Q(x)dx = − cos(x ) + C = − + .
µ(x) sin x 2 2 sin x sin x

2. Solve the initial-value problem

y 0 = x + y, y(0) = 2.

Answer: Re-write in the standard form for a linear equation:

y 0 − y = x.

Here P (x) = −1 and Q(x) = x, so

Z Z
P (x)dx = −dx = −x.

Hence R
µ(x) = e P (x)dx
= e−x .
In turn, Z Z
µ(x)Q(x)dx = e−x xdx = −xe−x − e−x + C

1
 (integrating by parts with u = x, dv = e−x dx). Therefore,
Z
1 1
µ(x)Q(x)dx = −x −xe−x − e−x + C = ex −xe−x − e−x + C = −x−1+Cex .



y=
µ(x) e
Therefore, plugging in x = 0, we see that

2 = −(0) − 1 + Ce0 = −1 + C,

so C = 3. Thus, the solution to the initial-value problem is

y = −x − 1 + 3ex .

3. Solve the initial-value problem

xy 0 = y + x2 sin x, y(π) = 0.

Answer: Re-write as
xy 0 − y = x2 sin x.
Then, after dividing everything by x, we get the standard form of a linear first-order equation:
1
y0 − y = x sin x.
x
Here P (x) = − x1 and Q(x) = x sin x, so
Z Z  
1 1
P (x)dx = − dx = − ln |x| = ln .
x |x|
Thus, “ ”
R
P (x)dx ln 1
|x|
1
µ(x) = e =e = .
|x|
In the statement of the problem I should have specified that x > 0, so we can ignore the
absolute value signs and let µ(x) = x1 . Therefore,
Z Z Z
1
µ(x)Q(x)dx = x sin xdx = sin xdx = − cos x + C.
x
Therefore, Z
1 1
y= µ(x)Q(x)dx = 1 (− cos x + C) = −x cos x + Cx.
µ(x) x
Plugging in x = π yields
0 = −π cos π + Cπ,
so
Cπ = π cos π = −π,
meaning that C = −1. Thus, the solutions to the initial-value problem is

y = −x cos x − x.

2
4. The figure shows a circuit containing an electromotive force, a capacitor with a capacitance
of C farads (F ), and a resistor with a resistance of R ohms (Ω). Ohm’s Law says that the
voltage drop across the resistor is RI. The voltage drop across the capacitor is Q/C, where
Q is the charge (in coulombs), so in this case Kirchhoff’s Law gives
Q
RI + = E(t).
C
By definition I = dQ/dt, so we have
dQ 1
R + Q = E(t).
dt C
Suppose the resistance is 5 Ω, the capacitance is 0.05 F , a battery gives a constant voltage of
60 V , and the initial charge is Q(0) = 0 coulombs. Find the charge and the current at time t.
C

E R

Answer: Plugging in the values for R, C, and E, we see that

dQ 1
5 + Q = 60,
dt 0.05
or
dQ
5
+ 20Q = 60.
dt
Dividing everything by 5 yields the linear equation in standard form
dQ
+ 4Q = 12.
dt
Here P (t) = 4 and Q(t) = 12, so
Z Z
P (t)dt = 4dt = 4t.

Thus, R
P (t)dt
µ(t) = e = e4t .
In turn, Z Z Z
4t
e4t 4dt = 3 e4t + C .


µ(t)Q(t)dt = e 12dt = 3

Therefore,
C0
Z
1
Q(t) = µ(t)Q(t)dt = 3 + .
µ(t) e4t

3
 Using what we know about the situation when t = 0, we can solve for C 0 :
C0
0=3+ = 3 + C 0,
e4·0
so C 0 = −3.
Hence,
3
Q(t) = 3 − .
e4t
5. An object with mass m is dropped from rest and we assume that the air resistance is pro-
portional to the speed of the object. If s(t) is the distance dropped after t seconds, then the
speed is v = s0 (t) and the acceleration is a = v 0 (t). If g is the acceleration due to gravity, then
the downward force on the object is mg − cv, where c is a positive constant, and Newton’s
Second Law gives
dv
m = mg − cv.
dt
(a) Solve this for v (note that this is a linear differential equation in v).
Answer: Re-write as
dv
m + cv = mg.
dt
Then dividing everything by m yields the linear equation in standard form
dv c
+ v = g.
dt m
c
Here P (t) = m and Q(t) = g, so
Z Z
c c
P (t)dt = dt = t.
m m
Thus,
c
R
P (t)dt
µ(t) = e = e m t.
In turn, this means that
Z Z
c mg c t
µ(t)Q(t)dt = e m t gdt = e m + C.
c
Therefore,
Z
1 1  mg c t  mg c
v(t) = µ(t)Q(t)dt = c
t
e m + C = + Ce− m t .
µ(t) em c c
Since the object starts at rest, v(0) = 0, so
mg c mg
0= + Ce m ·0 = + C.
c c
Hence, C = − mg
c , so
mg mg − c t mg  c

v(t) = − e m = 1 − e− m t .
c c c

4
 (b) What is the terminal velocity of the falling object?
Answer: We can solve for the terminal velocity by computing
mg  c
 mg  c
 mg mg
lim v(t) = lim 1 − e− m t = lim 1 − e− m t = (1 − 0) = .
t→∞ t→∞ c c t→∞ c c
mg
Hence, the terminal velocity of the falling object is c .
(c) Find the distance the object has fallen after t seconds.
Rt
Answer: Remember that v(t) = s0 (t), so s(t) = 0 v(τ )dτ . Thus, we compute that the
distance traveled after t seconds is
Z t Z t
mg  c
 mg h m c it mg m2 g − c t m2 g
s(t) = v(τ )dτ = 1 − e− m τ dτ = τ + e− m τ = t+ 2 e m − 2 .
0 0 c c c 0 c c c
6. Each system of differential equations is a model for two species that either compete for the
same resources or cooperate for mutual benefit (flowering plants and insect pollinators, for
instance). Decide whether each system describes competition or cooperation and explain why
it is a reasonable model. (Hint: think about the effect an increase in one species has on the
growth rate of the other.)
dx 2
(a) dt = 0.12x − 0.0006x + 0.00001xy
dy
dt = 0.08y + 0.00004xy
Answer: Increasing the y population causes an increase in dx dt , the rate of growth of x.
Likewise, an increase in the x population causes an increase in dy dt , the rate of growth
of y. Therefore, this system of equations describes the populations of two species which
are cooperating with each other.
dx
(b) dt = 0.15x − 0.0002x2 − 0.0006xy
dy
dt = 0.2y − 0.00008y 2 − 0.0002xy
Answer: Increasing the y population causes a decrease in the rate of growth of x, and
likewise an increase in the x population causes a decrease in the rate of growth of the y
population. Thus, these equations describe the populations of two species which are in
competition with each other.
7. Populations of aphids and ladybugs are modeled by the equations
dA
= 2A − 0.01AL
dt
dL
= −0.5L + 0.0001AL
dt
(a) Find the equilibrium solutions and explain their significance.
Answer: An equilibrium solution is a solutions for which both A and L are constant.
If A is constant, then dA
dt = 0, so the first equation becomes

0 = 2A − 0.01AL = A(2 − 0.01L),

dL
meaning either A = 0 or L = 200. If A = 0 and L is constant, then dt = 0, so the
second equation reduces to
0 = −0.5L + 0.0001(0)L = −0.5L,

5
 so L = 0. Thus, (A, L) = (0, 0) is one equilibrium solution. On the other hand, if
L = 200 and A is constant, then

0 = −0.5(200) + 0.0001A(200) = −100 + 0.02A,

so A = 5000. Therefore, (A, L) = (5000, 200) is the other equilibrium solution.

The trivial solution A = L = 0 just says that, if there are no aphids and no ladybugs to
start out with, there will continue to be none forever.
The solution (A, L) = (5000, 200) says that when there are 5000 aphids and 200 ladybugs,
the populations are just the right size that there are no changes in the size of either
population.
(b) Find an expression for dL/dA.
Answer: If L were a function of A, then we would have that
dL dL dA
= .
dt dA dt
Therefore,
dL
dL dt −0.5L + 0.0001AL
= dA
= .
dA dt
2A − 0.01AL

(c) The direction field for the differential equation in part (b) is shown [omitted from this
answer key]. Use it to sketch a phase portrait.
Answer: Below is a phase portrait of the system:

6
(d) Suppose that at time t = 0 there are 1000 aphids and 200 ladybugs. Draw the corre-
sponding phase trajectory and use it to describe how both populations change.
Answer: Below is the phase trajectory determined by this initial condition:

This says that in the beginning there aren’t enough aphids to support 200 ladybugs,
so the ladybug population decreases rapidly, while the aphid population increases very
slowly. As more ladybugs die, the rate of decrease of the population of ladybugs slows,
while the rate of growth of the aphid population increases.
After the ladybug population bottoms out at 100, the ladybug population starts increas-
ing; the aphid population is still growing rapidly, eventually getting to a maximum of
almost 15,000 aphids. Of course, with all these aphids to eat, the ladybug population
increases rapidly, eventually causing the aphids to start dying off.
The ladybugs reach a maximum population of more than 350, but by this point the
falling aphid population means that this is unsustainable. As the aphid population
bottoms out at 1000, the ladybug population falls through 200, starting the cycle anew.
(e) Use part (d) to make rough sketches of the aphid and ladybug populations as functions
of t. How are the graphs related to each other?
Answer: Both graphs have the same period and the graph of L peaks about a quarter
of a cycle after the graph of A.