ODE-Project First-Order Linear Equations

Section 1.5 First-Order Linear Equations

Objectives

To understand that any first-order linear differential equation

$y^{'} + p (x) y = g (x)$

can be solved by multiplying each side of the equation by an integrating factor.
To understand the existence and uniqueness of solutions to first-order initial value problems.

A first-order differential equation is an equation of the form

\begin{matrix} (1.5.1) & \frac{d x}{d t} + p (t) x = q (t) . \end{matrix}

🔗

This equation will not be separable if

p (t)

is not a constant. We shall have to find a new approach to solving such an equation. We could, of course, use a numerical algorithm to solve (1.5.1); however, we can always find an algebraic solution to a first-order linear differential equation. Moreover, the fact that we can obtain such a solution analytically will prove very useful when we investigate more complicated equations and systems of equations.

🔗

Subsection 1.5.1 Mine Tailings

🔗

In any mining operation, tailings are what is left after everything of value has been extracted. For example, in a hard rock mining operation, ore is often pulverized and then processed using chemicals to extract certain minerals of value. Soft rock mining operations such as coal mining or extracting oil from tar sands might use solvents or water to extract any commodity of value. The material that is left over after the minerals, coal, or oil is extracted can often present huge environmental challenges. There are different ways of processing mine tailings, but one way is to store them in a pond, especially if water is used in the mining operation. This method allows any particles that are suspended in the water to settle to the bottom of the pond. The water can then be treated and recycled.

🔗

Suppose that we have a gold mining operation and we are storing our tailings in a pond that has an initial volume of 20,000 cubic meters. When we begin our operation, the tailings pond is filled with clean water. The pond has a stream flowing into it, and water is also pumped out of the pond. Chemicals are used in processing gold ore. These chemicals such as sodium cyanide can be highly poisonous and dangerous to the environment, and the water must be treated before it is released into the watershed. Suppose that 1000 cubic meters per day flow into the pond from stream and 1000 cubic meters are pumped from the pond each day to be processed and recycled. Thus, the water level of the pond remains constant.

🔗

At time

t = 0,

the water from stream becomes contaminated with chemicals from the mining operation, say at a rate of 5 kilograms of chemicals per 1000 cubic meters. We will assume that water in our tailings pond is well mixed so that the concentration of chemicals through out the pond is fairly uniform. In addition, any particulate matter pumped into the pond from the stream settles to the bottom of the pond at a rate of 50 cubic meters per day. Thus, the volume of our tailings pond is reduced by 50 cubic meters each day, and our tailings pond will become full after 400 days of operation. We shall assume that the particulate matter and the chemicals are included in the 1000 cubic meters that flow into the pond from the stream each day.

🔗

We wish to find a differential equation that will model the amount of chemicals in the tailings pond at any particular time. Let

x (t)

be the amount of chemicals in the pond at time

t .

Then

d x / d t

is the difference between the rate at which the chemicals are entering the pond and the rate at which the chemicals leave the pond.

\frac{d x}{d t} = rate in - rate out .

🔗

Since water flows into the pond from the stream at a rate of 1000 cubic meters per day, the rate at which the chemicals enter the pond is 5 kilograms per day. On the other hand, the rate at which the chemicals leave the pond will depend on the amount of chemicals in the pond at time

t .

The volume of the pond is decreasing due to sediment, and at time

t

it is

V (t) = 20000 - 50 t .

Thus, the concentration of chemicals in the pond at time

t

x / (20000 - 50 t),

and the rate at which the chemicals are flowing out of the pond to be recycled is

1000 (\frac{x}{20000 - 50 t}) = \frac{20 x}{400 - t} .

🔗

Hence, the differential equation that models the amount of chemical in the tailings pond at time

t

\begin{matrix} (1.5.2) & \frac{d x}{d t} = 5 - \frac{20 x}{400 - t} . \end{matrix}

🔗

Of course, we will have to cease mining operations once the pond is full, since there will only be water in the pond if

V (t) = 20000 - 50 t \geq 0;

that is, when

0 \leq t < 400 .

🔗

Notice that equation (1.5.2) is not autonomous. In fact, it is not even separable. We will have to use a different approach to find a solution. First, we will rewrite the equation in the form

\frac{d x}{d t} + \frac{20}{400 - t} x = 5.

🔗

If we multiply both sides of this equation by

(400 - t)^{- 20},

we obtain

\begin{matrix} (1.5.3) & (400 - t)^{- 20} \frac{d x}{d t} + 20 (400 - t)^{- 21} x = 5 (400 - t)^{- 20} . \end{matrix}

🔗

We now make the crucial observation that the product rule applies to the lefthand side of our equation,

(400 - t)^{- 20} \frac{d x}{d t} + 20 (400 - t)^{- 21} x = \frac{d}{d t} ((400 - t)^{- 20} x) .

🔗

Thus, equation (1.5.3) becomes

\frac{d}{d t} ((400 - t)^{- 20} x) = 5 (400 - t)^{- 20} .

🔗

Integrating both sides of this equation, we have

(400 - t)^{- 20} x = 5 \int (400 - t)^{- 20} d t = \frac{5 (400 - t)^{- 19}}{19} + C,

🔗

where

C

is an arbitrary constant. Solving for

x,

we obtain

x (t) = \frac{5}{19} (400 - t) + C (400 - t)^{20} .

🔗

Since

x (0) = 0,

we can quickly determine that

C = - (5 / 19) 400^{- 19}

and that the solution to our initial value problem is

x (t) = \frac{5}{19} (400 - t) [1 - {(\frac{400 - t}{400})}^{19}] .

🔗

The graph of the solution to our differential equation (Figure 1.5.1) fits the situation. Initially, there are no chemicals in the pond, but

x (t)

quickly increases. However, the amount of chemicals decreases as the pond begins to fill with sediment. Eventually, there are no chemicals at

t = 400 .

described in detail following the image — 🔗
Figure 1.5.1. Chemical waste in a tailings pond

🔗

Subsection 1.5.2 First-Order Linear Equations

🔗

The differential equation

\frac{d x}{d t} + \frac{20}{400 - t} x = 5

🔗

is an example of a first-order linear differential equation. More specifically, a first-order linear differential equation is an equation that can be written in the form

\frac{d x}{d t} + p (t) x = q (t) .

🔗

Let us first show how to solve first-order linear equations when the coefficient functions are constant. If

\begin{aligned} \frac{d x}{d t} + p x & = q, \\ x (t_{0}) & = x_{0}, \end{aligned}

🔗

then we first multiply both sides of the equation by

e^{p t}

to obtain

e^{p t} \frac{d x}{d t} + e^{p t} p x = q e^{p t} .

🔗

The left-hand side of the equation is exact. That is,

e^{p t} \frac{d x}{d t} + e^{p t} p x = \frac{d}{d t} (x (t) e^{p t}) .

🔗

If we integrate both sides of

\frac{d}{d t} (x (t) e^{p t}) = q e^{p t},

🔗

then

x (t) e^{p t} = \frac{q}{p} e^{p t} + C .

🔗

If we apply the initial condition, we can determine

C,

C = (x_{0} - \frac{q}{p}) e^{p t_{0}} .

🔗

Thus, the solution that we seek is

x (t) = \frac{q}{p} + (x_{0} - \frac{q}{p}) e^{p (t_{0} - t)} .

Example 1.5.2.

Suppose we wish to solve the initial value problem

\begin{aligned} (1.5.4) & \frac{d x}{d t} - 2 x & = 3, \\ (1.5.5) & x (0) & = 1. \end{aligned}

Multiplying both sides of equation (1.5.4) by

e^{- 2 t},

we obtain

\frac{d}{d t} (e^{- 2 t} x) = e^{- 2 t} \frac{d x}{d t} - 2 e^{- 2 t} x = 3 e^{- 2 t} .

Integrating both sides of this last equation, gives us the following

e^{- 2 t} x = 3 \int e^{- 2 t} d t + C = - \frac{3}{2} e^{- 2 t} + C .

Applying the initial condition

x (0) = 1,

we can conclude that

C = 5 / 2,

and

x (t) = \frac{5}{2} e^{2 t} - \frac{3}{2} .

Example 1.5.3.

Now let us solve a first-order linear differential equation where the coefficients are not constant. Suppose that

y^{'} + 2 x y = e^{- x^{2}},

where

y (0) = 1 .

We will multiply both sides of the equation by

e^{x^{2}}

and use the product rule to obtain

\frac{d}{d x} (e^{x^{2}} y) = e^{x^{2}} y^{'} + 2 x e^{x^{2}} y = e^{x^{2}} e^{- x^{2}} = 1.

Integrating both sides, we get

e^{x^{2}} y = \int d x + C .

Thus, the general solution is

y = (x + C) e^{- x^{2}} .

Using the initial condition to solve for

C,

we find that

C = 1

and

y = (x + 1) e^{- x^{2}} .

🔗

Surprisingly, the strategy in Example 1.5.3 will always work. Suppose that

\begin{matrix} (1.5.6) & x^{'} + p (t) x = q (t) . \end{matrix}

🔗

If we choose

P (t)

such that

P^{'} (t) = p (t)

and multiply both sides of the equation by

e^{P (t)},

then

{(e^{P (t)} x)}^{'} = e^{P (t)} (x^{'} (t) + p (t) x (t)) = e^{P (t)} q (t) .

🔗

Integrating both sides,

e^{P (t)} x = \int e^{P (t)} q (t) d t + C

🔗

x = \frac{1}{e^{P (t)}} (\int e^{P (t)} q (t) d t + C) .

🔗

The Fundamental Theorem of Calculus tells us that

P (t) = \int p (t) d t .

We say that

μ (t) = e^{\int p (t) d t} = \exp (\int p (t) d t)

🔗

is an integrating factor for the differential equation (1.5.6).

Example 1.5.4.

Consider the initial value problem

\begin{matrix} \frac{d y}{d x} - \frac{2}{x} y = 2 x^{2} \\ y (- 2) = 4. \end{matrix}

Our integrating factor is

μ (x) = \exp (- \int \frac{2}{x} d x) = e^{- 2 \ln x} = e^{\ln x^{- 2}} = x^{- 2} .

Multiplying both of our differential equation by

μ (x) = x^{- 2},

we obtain

x^{- 2} (\frac{d y}{d x} - \frac{2}{x} y) = 2

\frac{d}{d x} (x^{- 2} y) = x^{- 2} \frac{d y}{d x} - 2 x^{- 3} y = 2.

We can now integrate this equation to get

x^{- 2} y = 2 x + C .

The initial condition

y (- 2) = 4

allows us to find

C = 5 .

Therefore, the solution to our initial value problem is

y = 2 x^{3} + 5 x^{2} .

🔗

Activity 1.5.1. Finding Solutions to First-Order Linear Differential Equations.

🔗

Solve each of the following initial value problems.

🔗

(a)

🔗

y^{'} = 3 y + t;

y (0) = 1

🔗

(b)

🔗

x^{'} = t x + 2 t;

x (1) = 2

🔗

(c)

🔗

y^{'} = 2 x - 3 y + e^{- x};

y (0) = 1

🔗

(d)

🔗

\frac{d u}{d t} = \frac{t - 1}{t} u + t^{2};

u (1) = 0

🔗

(e)

🔗

y^{'} = \sin x - y;

y (0) = 1

🔗

Subsection 1.5.3 Mixing Models

🔗

Many applications involve the mixing of two or more substances together. As we mentioned previously, we can model how petroleum products are mixed together in a refinery, how various ingredients are mixed together in a brewery, or how greenhouse gases move across various layers of the earth’s atmosphere.

Example 1.5.5.

Suppose that a 100-gallon tank initially contains 50 gallons of salt water containing five pounds of salt. A brine mixture containing one pound of salt per gallon flows into the top of the tank at a rate of 5 gallons per minute. A well mixed solution leaves the tank at rate of 4 gallons per minute. We wish to know how much salt is in the tank, when the tank is full.

To construct our model, we will let

t

be the time (measured in minutes) and set up a differential equation that will measure how fast the amount of salt at time

t,

x (t),

is changing. We have the initial condition

x (0) = 5,

and

\begin{aligned} \frac{d x}{d t} & = rate of salt flowing in - rate of salt flowing out \\ = \underset{in flow}{\underset{⏟}{5}} - \underset{out flow}{\underset{⏟}{4 \frac{x}{V (t)}}}, \end{aligned}

where

V (t)

is the volume at time

t .

The expression

x / V (t)

is the amount of salt in one gallon at time

t .

We have

V (t) = 50 + t,

since the tank starts with 50 gallons and five gallons are pumped into the tank per minute while four gallons leave the tank during the same time interval. Thus, our differential equation becomes

\frac{d x}{d t} = 5 - \frac{4}{50 + t} x .

Our equation is linear since we can rewrite it as

\begin{matrix} (1.5.7) & \frac{d x}{d t} + \frac{4}{50 + t} x = 5. \end{matrix}

An integrating factor for this differential equation is

μ (t) = \exp (\int \frac{4}{50 + t} d t) = e^{4 \ln (50 + t)} = (50 + t)^{4} .

Therefore, if we multiply both sides of equation (1.5.7) by

μ (t),

we get

(50 + t)^{4} \frac{d x}{d t} + 4 (t + 50)^{3} x = 5 (50 + t)^{4} .

We can now apply the product rule to obtain

\frac{d}{d t} [(50 + t)^{4} x] = 5 (50 + t)^{4} .

Integrating both sides and simplifying gives us

x = t + 50 + \frac{C}{(t + 50)^{4}} .

Our initial condition,

x (0) = 5

tells us that

C = - 281,250,000

and

x (t) = t + 50 - \frac{281250000}{(t + 50)^{4}} .

Thus, when the tank is full,

t = 50

and the amount of salt in the tank is

x (50) = 97.188

pounds.

🔗

Activity 1.5.2. A Mixing Problem.

🔗

Suppose that a tank contains 1000 gallons of a solution consisting of 200 pounds of salt dissolved in water. Pure water is pumped into the tank at a rate of 6 gallons per minute. At the same time, the tank is drained at the same rate. Assume that the brine mixture is kept well stirred.

🔗

(a)

🔗

Set up an initial value problem to model the amount of salt in the tank at time

t .

🔗

(b)

🔗

How long will it take until there is only 20 pounds of salt left in the tank?.

🔗

Subsection 1.5.4 Finance Models

🔗

There are a number of problems in finance that can be modeled using differential equations. Let

P (t)

be the balance of an account at time

t

and suppose that the account pays interest at a rate of

r

percent per year compounded continuously. Suppose that we also allow withdrawals of

W

dollars per year. The net increase in the balance between times

t

and

t + Δ t

can now be described as

P (t + Δ t) - P (t) \approx r P Δ t - W Δ t

🔗

Thus,

P^{'} (t) = lim_{Δ t \to 0} \frac{P (t + Δ t) - P (t)}{Δ t} = r P - W .

🔗

We can solve the equation

P^{'} = r P - W

🔗

by multiplying both sides of the equation by the integrating factor

μ (t) = \exp (- \int r d t) = e^{- r t} .

🔗

to obtain

\frac{d}{d t} [e^{- r t} P] = - W e^{- r t} .

🔗

Integrating both sides of this equation, we have

e^{- r t} P = \frac{W}{r} e^{- r t} + C

🔗

P = \frac{W}{r} + C e^{r t} .

🔗

If we know the initial balance in the account, say

P (0) = P_{0},

we can determine

C .

That is,

P_{0} = \frac{W}{r} + C

🔗

C = P_{0} - \frac{W}{r} .

🔗

Thus, the solution to our initial value problem is

P = \frac{W}{r} + (P_{0} - \frac{W}{r}) e^{r t} .

Example 1.5.6.

Suppose that your parents have established a money market account with a balance of $50,000 that they will use to help you pay for your college education. The account receives an average annual interest of 4%. You estimate that your tuition, room and board, and other college expenses to be $20,000 per year.

We model this financial situation with the differential equation

\begin{aligned} \frac{d P}{d t} & = 0.04 P - 20000 \\ P (0) & = 50000. \end{aligned}

Rewriting the differential equation as

d P / d t - 0.04 P = - 20000,

our integrating factor becomes

μ (t) = e^{- 0.04 t},

and

e^{- 0.04 t} \frac{d P}{d t} - 0.04 e^{- 0.04 t} P = \frac{d}{d t} [e^{- 0.04 t} P] = - 20000 e^{- 0.04 t} .

Thus,

e^{- 0.04 t} P = \frac{20000}{0.04} e^{- 0.04 t} + C = 500000 e^{- 0.04 t} + C .

The solution to this initial value problem is

P (t) = \frac{20000}{0.04} + (50000 - \frac{20000}{0.04}) e^{0.04 t} = 500000 - 450000 e^{0.04 t} .

Your parents have been quite generous but have told you that you must be responsible for the balance of the cost of your education.

🔗

Activity 1.5.3. Paying for College.

🔗

Suppose that new parents want to start a college fund for their child. They are willing to invest $2000 per year at a rate of 4%.

🔗

(a)

🔗

Find an initial value problem that models the parents’ investment.

🔗

(b)

🔗

How much will be in the college fund when their child turns 18?

🔗

(c)

🔗

What would they need to invest per year to have $80,000 in the college fund when their child turns 18?

🔗

Subsection 1.5.5 Existence and Uniqueness of Solutions

🔗

Several questions about the existence and uniqueness of solutions to first-order linear differential equations now arise.

Does an initial value problem always have a solution?
Is the solution unique?
Is the solution globally defined or does it only hold for a small interval?

🔗

We can use the following theorem to answer these questions.

🔗

Theorem 1.5.7.

🔗

y^{'} + p (t) y = g (t)

🔗

is a differential equation such that

y (t_{0}) = y_{0},

and

p (t)

and

g (t)

are continuous on the open interval

I = (α, β),

then there exists a unique function

y = ϕ (t)

satisfying the differential equation and the initial condition on

I .

Proof.

μ (t) = \exp (\int p (t) d t),

then

\frac{d}{d t} (μ (t) y) = μ (t) (\frac{d y}{d t} + p (t) y) = μ (t) g (t) .

Integrating both sides of this equation and solving for

y,

we have

y = \frac{1}{μ (t)} (\int μ (t) g (t) d t + C) .

Since

y (t_{0}) = y_{0},

the constant

C

is uniquely determined. Notice that we have used continuity to guarantee that the integrals exist.

🔗

Subsection 1.5.6 Important Lessons

🔗

Any first-order linear differential equation

$y^{'} + p (x) y = g (x)$

can be solved by multiplying each side of the equation by an integrating factor

$μ (s) = e^{\int p (x) d x} .$

In this case, we get

${(e^{P (x)} y)}^{'} = e^{P (x)} (y^{'} (x) + p (x) y (x)) = e^{P (x)} g (x) .$

Integrating both sides,

$e^{P (x)} y = \int e^{P (x)} g (x) d x$

or

$y = \frac{1}{e^{P (x)}} \int e^{P (x)} g (x) d x .$
If

$y^{'} + p (t) y = g (t)$

is a differential equation such that $y (t_{0}) = y_{0},$ and $p (t)$ and $g (t)$ are continuous on the open interval $I = (α, β),$ then there exists a unique function $y = ϕ (t)$ satisfying the differential equation and the initial condition on $I .$

🔗

Reading Questions 1.5.7 Reading Questions

🔗

1.

🔗

Explain in your own words what a first-order linear differential equation is.

🔗

2.

🔗

What important rule from differential calculus do we use when solving a first-order differential equation?

🔗

Exercises 1.5.8 Exercises

🔗

Finding General Solutions.

🔗

Find the general solution for each equation in Exercise Group 1.5.8.1–10.

🔗

1.

🔗

\frac{d y}{d x} + 5 y = 0

🔗

2.

🔗

x^{'} - 7 x = 0

🔗

3.

🔗

y^{'} + 2 x y = 0

🔗

4.

🔗

\frac{d x}{d t} - t^{2} x = 0

🔗

5.

🔗

y^{'} + \frac{2}{x} y = 0

🔗

6.

🔗

\frac{d y}{d x} - 5 y = e^{x}

🔗

7.

🔗

y^{'} - 5 y = 10 x

🔗

8.

🔗

\frac{d x}{d t} - 5 x = \sin 2 t

🔗

9.

🔗

y^{'} - \frac{3}{x^{2}} y = \frac{1}{x^{2}}

🔗

10.

🔗

\frac{d y}{d x} + \frac{2}{x} y = \frac{\sin x}{x^{2}}

🔗

Solving Initial Value Problems.

🔗

Solve the initial value problems in Exercise Group 1.5.8.11–20.

🔗

11.

🔗

\frac{d y}{d x} + 5 y = 0,

y (0) = 2

🔗

12.

🔗

x^{'} - 7 x = 0,

x (0) = 1

🔗

13.

🔗

y^{'} + 2 x y = 0,

y (0) = 3

🔗

14.

🔗

\frac{d x}{d t} - t^{2} x = 0,

x (0) = - 1

🔗

15.

🔗

y^{'} + \frac{2}{x} y = 0,

y (1) = - 3

🔗

16.

🔗

y^{'} = - \frac{2 y}{x + 1} + 2 x,

y (0) = 2

Hint.

y = \frac{3 x^{4} + 8 x^{3} + 6 x^{2} + 12}{6 (x + 1)^{2}}

🔗

17.

🔗

y^{'} = - \frac{2 y}{1 + x} + e^{x},

y (0) = 1

Hint.

y = (e^{x} + 1) / (x + 1)^{2}

🔗

18.

🔗

\frac{d x}{d t} - 5 x = \sin 2 t,

x (0) = π

🔗

19.

🔗

y^{'} = y \tan x + \frac{e^{x}}{\cos x},

y (0) = 1

Hint.

y = e^{x} \sec x

🔗

20.

🔗

y^{'} = - \frac{y}{x + 2} + \frac{\cos x}{x + 2},

y (0) = 1

Hint.

y = \sin x / (x + 2)

🔗

21.

🔗

A 600-liter tank initially contains 200 liters of water containing 10 kilograms of salt. Supposed that water containing 0.1 kilograms of salt flows into the top of the tank at a rate of 10 liters per minute. The water in the tank is kept well mixed, and 5 liters per minute are removed from the bottom of the tank. How much salt is in the tank when the tank is full?

Hint.

x (t)

is the amount of salt in the tank at time

t,

we know that

x (0) = 10 .

The volume of the tank is

V = 200 + 5 t .

We can model the amount of salt in the tank at time

t

with a differential equation,

\begin{aligned} \frac{d x}{d t} & = rate in - rate out \\ = 10 (0.1) - 5 \frac{x}{V} \\ = 1 - 5 \frac{x}{200 + 5 t} \\ = 1 - \frac{x}{40 + t} . \end{aligned}

The resulting equation

\frac{d x}{d t} + \frac{1}{40 + t} x = 1

is a first order linear differential equation. An integrating factor for this equation is given by

μ (t) = \exp (\int \frac{1}{40 + t} d t) = 40 + t .

Multiplying both sides of the differential equation by

μ (t),

we have

\frac{d}{d t} [(40 + t) x] = (40 + t) \frac{d x}{d t} + x = (40 + t) (\frac{d x}{d t} + \frac{1}{40 + t} x) = 40 + t .

Integrating both sides of this equation, we obtain

(40 + t) x = 40 t + \frac{t^{2}}{2} + C .

Using the intial condition

x (0) = 10,

we can determine that

C = 400

x (t) = \frac{t^{2} + 80 t + 800}{2 t + 80} .

The tank is full at time

t = 400 / 5 = 80,

and the tank contains

x (80) = 170 / 3 \approx 56.67

kilograms of salt when the tank is full.

🔗

22.

🔗

A manager in a communications company contributes $2400 per year into her retirement fund by making many small deposits throughout the year. The fund grows at a rate of 3.5% per year compounded continuously. After 35 years, she retires and begins and begins withdrawing from the retirement fund at a rate of $3500 per month. If she does not make any deposits after she retires, how long will her retirement fund last? [Hint: Divide the problem into two smaller problems—one that deals with the situation before retirement and one that deals with the problem after retirement.]

🔗

23.

🔗

Lake Baikal, located in southern Siberia, is only the seventh largest lake in the world in terms of surface area; however, it is the world’s deepest lake. The lake has a depth of 1,642 meters, and the bottom lies 1,186.5 meters below sea level. Lake Baikal has a volume of 23,600 cubic kilometers and contains 20% of the world’s unfrozen fresh water. In contrast, Lake Superior, the largest of the Great Lakes, has a volume of only 12,100 cubic kilometers. Although 544 rivers flow into Lake Baikal, there is only a single outlet, the Angara River. The outflow of the lake is pretty constant at 60.4 cubic kilometers per year.

🔗

Pollution is an increasing concern in Lake Baikal. One of the major polluters has been Baykalsk Pulp and Paper Mill. The mill was constructed in 1966 on the shoreline of Lake Baikal and regularly discharged waste into the lake. The plant was closed in November 2008 due to unprofitability, but production resumed in January 2010. The mill underwent a final bankruptcy in September 2013, but the longterm fate of the mill is still undecided.

Levy, Clifford J. (November 8, 2010). “Last Gasp for Factory Bequeathed by Soviets.” The New York Times. Retrieved March 14, 2014 from www.nytimes.com/2010/11/09/world/europe/09baikal.html.

🔗

Suppose that we wish to understand how the pollution level changes in Lake Baikal over a period of years. Hypothetically, let us assume that the Baykalsk Pulp and Paper Mill has been responsible for most of the pollution in Lake Baikal for the last several decades. Suppose that at

t = 0

years the mill ceases operation and there are no longer any pollutants discharged into the lake from the mill although there are still other sources of pollution. Let us assume that the lake is currently 6 times more polluted than these other sources of contaminants. We wish to know how long it will take for the pollution level to reduce to half of its current level of lake. Lake Baikal’s waters are well-mixed and well-oxygenated in spite of its great depth, so we can model this situation as a simple mixing problem.

The volume of water in Lake Baikal is $V = 23, 600 {km}^{3},$ and

$r_{in} = r_{out} = r = 60 {km}^{3} / year$

be the rates of inflow from the numerous rivers that feed the lake and outflow to the Angara River. Assume that $C = C_{in}$ is the pollutant concentration flowing into Lake Baikal and $C_{out}$ is the concentration of the outflow into the Angara River (measured in metric tons per cubic km). If $x (t)$ is the amount of solute at time $t$ in Lake Baikal, then

$x_{0} = x (0) = 6 C V$

is the initial amount of solute in the lake. Estimate $Δ x$ during the time interval $[t, t + Δ t],$ where $Δ t > 0$ is small.
From your estimate of $Δ x$ in part (1), write an initial value problem that describes the amount of pollutants in the lake at time $t .$
The equation that you found in part (2) is a first-order linear equation. Solve this equation.
Using part (3), predict how many years it will take to reduce the pollution in Lake Baikal to half of its current level.

🔗

24.

🔗

Variation of Parameters. Consider the following method of solving the general linear equation of the first order,

\begin{matrix} (1.5.8) & y^{'} + p (t) y = g (t) . \end{matrix}

If $g (t)$ is identically zero, show that the solution is

$y = A \exp [- \int p (t) d t],$

where $A$ is a constant.
If $g (t)$ is not identically zero, assume that the solution is of the form

$\begin{matrix} (1.5.9) & y = A (t) \exp [- \int p (t) d t], \end{matrix}$

where $A$ is now a function of $t .$ By substituting for $y$ in the given differential equation (1.5.8), show that $A (t)$ must satisfy the condition

$\begin{matrix} (1.5.10) & A^{'} (t) = g (t) \exp [\int p (t) d t] . \end{matrix}$
Find $A (t)$ from Equation (1.5.10). Then substitute for $A (t)$ in equation (1.5.9) and determine $y .$ Verify that the solution obtain in this manner agrees with the solution given in the proof of Theorem Theorem 1.5.7. That is, show that this solution is equivalent to the solution

$y = \frac{1}{μ (t)} (\int μ (t) g (t) d t + C),$

where

$μ (t) = \exp (\int p (t) d t) .$

🔗

T his technique is know as variation of parameters, which we will revisit when we study second order linear differential equations.

🔗

25.

🔗

Bernoulli’s equation is an equation of the form

y^{'} = a (t) y + f (t) y^{n},

🔗

where

n \neq 0

or 1. Bernoulli’s equation is nonlinear and cannot be solved by the techniques that we have used to solve first order linear equations.

Polking p. 63

Using the substitution $z = y^{1 - n},$ show that we can transform Bernoulli’s equation into the linear equation

$z^{'} = (1 - n) a (t) z + (1 - n) f (t) .$
Solve the equation $x y^{'} + y = x^{4} y^{3} .$

🔗

26.

🔗

The first-order nonlinear differential equation

\begin{matrix} (1.5.11) & y^{'} = p (t) + q (t) y + r (t) y^{2} \end{matrix}

🔗

is known as the Ricatti equation and has some useful applications in control theory. If one solution

y_{1} (t)

of the Ricatti equation is known, then a more general solution containing an arbitrary constant can be found by substituting

y = y_{1} (t) + 1 / v (t)

into equation (1.5.11) to find a first-order linear equation in

v

and

t,

which we can solve to find a general solution to the Ricatti equation.

Show that this first-order linear equation is $v^{'} + [q (t) + 2 r (t) y_{1} (t)] v = - r (t) .$
Find the solution to the Ricatti equation

$y^{'} = - \frac{1}{t^{2}} - \frac{1}{t} y + y^{2}$

given the particular solution $y_{1} = 1 / t .$
Find the solution to the Ricatti equation

$y^{'} = \cos t - y \tan t + y^{2} \sec t$

given the particular solution $y_{1} = \sin t .$
Find the solution to the Ricatti equation

$y^{'} = 2 - 3 y + y^{2}$

given the particular solution $y_{1} \equiv 2 .$

Hint.

If $y = y_{1} + 1 / v,$ then $y^{'} = y_{1}^{'} - v^{'} / v^{2} .$ Substituting into our original equation, we obtain

$y^{'} = y_{1}^{'} - \frac{v^{'}}{v^{2}} = p + q y_{1} + r y_{1}^{2} - \frac{v^{'}}{v^{2}} .$

On the other hand,

$\begin{aligned} y^{'} & = p + q (y_{1} + \frac{1}{v}) + r {(y_{1} + \frac{1}{v})}^{2} \\ = p + q y_{1} + \frac{q}{v} + r y_{1}^{2} + \frac{2 r y_{1}}{v} + \frac{r}{v^{2}} \\ = y_{1}^{'} + \frac{q}{v} + \frac{2 r y_{1}}{v} + \frac{r}{v^{2}} . \end{aligned}$

Therefore,

$- \frac{v^{'}}{v^{2}} = \frac{q}{v} + \frac{2 r y_{1}}{v} + \frac{r}{v^{2}},$

which is just the first-order linear equation

$v^{'} + [q (t) + 2 r (t) y_{1} (t)] v = - r (t) .$
$y = t + \frac{1}{C - t}$
$y (t) = \frac{1}{C \cos t - \sin t} + \sin t$
$y (t) = 2 + \frac{1}{C e^{t} - 1}$

🔗

27.

🔗

Suppose that we have a population that not only grows logistically but also requires a minimum threshold population to survive. For example, the case of the North Pacific right whale, a species now very much on the endangered list. If the population drops too low, whales might not be able to find suitable mates and the species will eventually go extinct. In other words, the population will die out if it drops below a certain threshold. We can model this with the following equation,

\begin{matrix} (1.5.12) & \frac{d P}{d t} = k (1 - \frac{P}{N}) (P - a N), \end{matrix}

🔗

where

P

is the population of the whales at time

t

and

N

is the carrying capacity. The constants

k

and

a

are positive with

a < 1 .

Find the equilibrium solutions of this equation.
Since equation (1.5.12) is autonomous, we can find a solution using separation of variables. Find this solution.
Equation (1.5.12) is also a Ricatti equation (1.5.11). Since we know an equilibrium solution from part (1), we can use the method of the previous problem to find a general solution to (1.5.12). Find the general solution using the fact that we have a Ricatti equation and show that your solution agrees with the solution that you found in part (2).

Prev Top Next