1.4 Linear equations and the integrating factor

One of the most important types of equations we will learn how to solve are the so-called linear equations. In fact, the majority of the course is about linear equations. In this lecture we focus on the first order linear equation. A first order equation is linear if we can put it into the form:

$\label{eq:1}y' + p(x)y = f(x).$

(1.4.1)

Here the word “linear” means linear in $y$ and $y'$ ; no higher powers nor functions of $y$ or $y'$ appear. The dependence on $x$ can be more complicated.

Solutions of linear equations have nice properties. For example, the solution exists wherever $p(x)$ and $f(x)$ are defined, and has the same regularity (read: it is just as nice). But most importantly for us right now, there is a method for solving linear first order equations. The trick is to rewrite the left hand side of (1.4.1) as a derivative of a product of $y$ with another function. To this end we find a function $r(x)$ such that

$r(x)y' + r(x)p(x)y = \frac{d}{dx}\left[r(x)y\right] \nonumber$

This is the left hand side of (1.4.1) multiplied by $r(x)$ . So if we multiply (1.4.1) by $r(x)$ , we obtain

$\frac{d}{dx}\left[r(x)y\right] = r(x)f(x) \nonumber$

Now we integrate both sides. The right hand side does not depend on $y$ and the left hand side is written as a derivative of a function. Afterwards, we solve for $y$ . The function $r(x)$ is called the integrating factor and the method is called the integrating factor method.

We are looking for a function $r(x)$ , such that if we differentiate it, we get the same function back multiplied by $p(x)$ . That seems like a job for the exponential function! Let

$r(x) = e^{\int p(x)dx} \nonumber$

We compute:

$\begin{align}\begin{aligned} y' + p(x)y &= f(x), \\ e^{\int p(x)dx}y' + e^{\int p(x)dx}p(x)y &= e^{\int p(x)dx}f(x), \\ \frac{d}{dx}\left[e^{\int p(x)dx}y\right] &= e^{\int p(x)dx}f(x), \\ e^{\int p(x)dx}y &= \int e^{\int p(x)dx}f(x)dx + C, \\ y &= e^{-\int p(x)dx}\left( \int e^{\int p(x)dx}f(x)dx + C \right).\end{aligned}\end{align} \nonumber$

Of course, to get a closed form formula for $y$ , we need to be able to find a closed form formula for the integrals appearing above.

Example 1

Solve

$y' + 2xy = e^{x-x^2}, \quad y(0) = -1 \nonumber$

Show solution

First note that $p(x) = 2x$ and $f(x) = e^{x-x^2}$ . The integrating factor is $r(x) = e^{\int p(x)dx} = e^{x^2}$ . We multiply both sides of the equation by $r(x)$ to get

x^{2}

The integrating-factor exponent: the integral of 2x dx ✓ Computed · mojocas 0.1.0 ✓ Agrees with the textFor y' + 2xy = e^(x - x^2) the text takes p(x) = 2x and states the integrating factor r(x) = e^(integral of p) = e^(x^2). The exponent is this integral.

$\begin{align}\begin{aligned} e^{x^2}y' + 2xe^{x^2}y &= e^{x-x^2}e^{x^2}, \\ \frac{d}{dx}\left [ e^{x^2}y\right] &= e^{x}.\end{aligned}\end{align} \nonumber$

We integrate

$\begin{align} \begin{aligned} e^{x^2} y &= e^x +C, \\ y &= e^{x-x^2} + C e^{-x^2}.\end{aligned}\end{align} \nonumber$

Next, we solve for the initial condition

-1 = y(0) = 1 + C

, so

C=-2

. The solution is

$y = e^{x-x^2} - 2 e^{-x^2}. \nonumber$

Note that we do not care which antiderivative we take when computing $e^{\int p(x) dx}$ . You can always add a constant of integration, but those constants will not matter in the end.

Since we cannot always evaluate the integrals in closed form, it is useful to know how to write the solution in definite integral form. A definite integral is something that you can plug into a computer or a calculator. Suppose we are given

$y' + p(x)y = f(x), \quad y(x_0) = y_0 \nonumber$

. Look at the solution and write the integrals as definite integrals.

$\label{eq:17}y(x) = e^{\int-_{x_0}^x p(s)ds} \left(\int_{x_0}^x e^{\int_{x_0}^t p(s)ds} f(t) dt + y_0 \right)$

(1.4.2)

You should be careful to properly use dummy variables here. If you now plug such a formula into a computer or a calculator, it will be happy to give you numerical answers.

Your Turn

Write the solution of the following problem as a definite integral, but try to simplify as far as you can. You will not be able to find the solution in closed form.

$y' + y = e^{x^2-x}, \quad y(0) = 10 \nonumber$

Example 2

Let us discuss a common simple application of linear equations. This type of problem is used often in real life. For example, linear equations are used in figuring out the concentration of chemicals in bodies of water (rivers and lakes).

Figure

\PageIndex{1}

A $100$ liter tank contains $10$ kilograms of salt dissolved in $60$ liters of water. Solution of water and salt (brine) with concentration of $0.1$ kilograms per liter is flowing in at the rate of $5$ liters a minute. The solution in the tank is well stirred and flows out at a rate of $3$ liters a minute. How much salt is in the tank when the tank is full?

Show solution

Let us come up with the equation. Let $x$ denote the kilograms of salt in the tank, let $t$ denote the time in minutes. For a small change $\Delta t$ in time, the change in $x$ (denoted $\Delta x$ ) is approximately

$\Delta x \approx \text{(rate in x concentration in)}\Delta t - \text{(rate out x concentration out)}\Delta t. \nonumber$

Dividing through by $\Delta t$ and taking the limit $\Delta t \rightarrow 0$ we see that

$\frac{dx}{dt} = \text{(rate in x concentration in)} - \text{(rate out x concentration out)} \nonumber$

In our example, we have

$\begin{align}\begin{aligned} \text{rate in} &= 5, \\ \text{concentration in} &= 0.1, \\ \text{rate out} &= 3, \\ \text{concentration out} &= \frac{x}{\text{volume}} = \frac{x}{60 + (5 - 3)t}.\end{aligned}\end{align} \nonumber$

Our equation is, therefore,

$\frac{dx}{dt} = \text{(5 x 0.1)} - \left ( 3 \frac{x}{60 + 2t} \right) \nonumber$

Or in the form (1.4.1)

$\frac{dx}{dt} + \frac{3}{60 + 2t}x = 0.5 \nonumber$

Let us solve. The integrating factor is

$r(t) = \text{exp} \left ( \int \frac{3}{60 + 2t} dt \right ) = \text{exp} \left ( \frac{3}{2} \ln (60 + 2t) \right ) = {(60 + 2t)}^{3/2} \nonumber$

We multiply both sides of the equation to get

$\begin{align}\begin{aligned} (60 + 2t)^{3/2} \frac{dx}{dt} + (60 + 2t)^{3/2} \frac{3}{60 + 2t}x &= 0.5(60 + 2t)^{3/2}, \\ \frac{d}{dt} \left [ (60 +2t)^{3/2}x \right ] &= 0.5(60 + 2t)^{3/2}, \\ (60 + 2t)^{3/2} x&= \int 0.5(60 + 2t)^{3/2}dt + C, \\ x &= (60 + 2t)^{-3/2} \int \frac {(60 + 2t)^{3/2}}{2}dt + C(60 + 2t)^{-3/2}, \\ x &= (60 + 2t)^{-3/2} \frac{1}{10} (60 + 2t)^{5/2} + C(60 + 2t)^{-3/2}, \\ x &= \frac {(60 + 2t)}{10} + C(60 + 2t)^{-3/2}.\end{aligned}\end{align} \nonumber$

We need to find $C$ . We know that at $t = 0$ , $x = 10$ . So

$10 = x(0) = \frac {60}{10} + C{(60)}^{-3/2} = 6 + C{(60)}^{-3/2} \nonumber$

$C = 4({60}^{3/2}) \approx 1859.03. \nonumber$

We are interested in $x$ when the tank is full. So we note that the tank is full when $60 + 2t = 100$ , or when $t = 20$ . So

$\begin{align}\begin{aligned} x(20) &= \frac{60 + 40}{10} + C{(60 + 40)}^{-3/2}\\ & \approx 10 + 1859.03{(100)}^{-3/2} \approx 11.86.\end{aligned}\end{align} \nonumber$

See Figure $\PageIndex{2}$ for the graph of $x$ over $t$ .

Figure

\PageIndex{2}

: Graph of the solution

x

kilograms of salt in the tank at time

t

The concentration at the end is approximately $0.1186$ ^kg/_liter and we started with $\frac{1}{6}$ or $0.167$ ^kg/_liter.

Adapted from Differential Equations for Engineers by Jiří Lebl (https://www.jirka.org/diffyqs/), © Jiří Lebl, licensed under CC BY-SA 4.0. Changes were made. License: CC-BY-SA-4.0.