# 1.4: Linear equations and the integrating factor

- Page ID
- 333

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\)

\( \newcommand{\id}{\mathrm{id}}\)

\( \newcommand{\Span}{\mathrm{span}}\)

\( \newcommand{\kernel}{\mathrm{null}\,}\)

\( \newcommand{\range}{\mathrm{range}\,}\)

\( \newcommand{\RealPart}{\mathrm{Re}}\)

\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

\( \newcommand{\Argument}{\mathrm{Arg}}\)

\( \newcommand{\norm}[1]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)

\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)

\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vectorC}[1]{\textbf{#1}} \)

\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

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). \]

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 \(\eqref{eq: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 \(\eqref{eq:1}\) multiplied by \(r(x)\). So if we multiply \(\eqref{eq: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 \]

Solve

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

**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

\[\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.

Try it! Add a constant of integration to the integral in the integrating factor and show that the solution you get in the end is the same as what we got above. An advice: Do not try to remember the formula itself, that is way too hard. It is easier to remember the process and repeat it.

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) \]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.

Check that \(y(x_0) = y_0\) in formula \(\eqref{eq:17}\).

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 \]

Before we move on, we should note some interesting properties of linear equations. First, for the linear initial value problem \(y' + p(x) y = f(x)\), \(y(x_0) = y_0\), there is always an explicit formula \(\eqref{eq:17}\) for the solution. Second, it follows from the formula \(\eqref{eq:17}\) that if \(p(x)\) and \(f(x)\) are continuous on some interval \((a,b)\), then the solution \(y(x)\) exists and is differentiable on \((a,b)\). Compare with the simple nonlinear example we have seen previously, \(y'=y^2\), and compare to Theorem 1.2.1.

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).

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?

**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 \(\eqref{eq: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 \]

or

\[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\).

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