Separable differential equations

$\newenvironment {prompt}{}{} \newcommand {\ungraded }[0]{} \newcommand {\todo }[0]{} \newcommand {\oiint }[0]{{\large \bigcirc }\kern -1.56em\iint } \newcommand {\mooculus }[0]{\textsf {\textbf {MOOC}\textnormal {\textsf {ULUS}}}} \newcommand {\npnoround }[0]{\nprounddigits {-1}} \newcommand {\npnoroundexp }[0]{\nproundexpdigits {-1}} \newcommand {\npunitcommand }[1]{\ensuremath {\mathrm {#1}}} \newcommand {\RR }[0]{\mathbb R} \newcommand {\R }[0]{\mathbb R} \newcommand {\N }[0]{\mathbb N} \newcommand {\Z }[0]{\mathbb Z} \newcommand {\sagemath }[0]{\textsf {SageMath}} \newcommand {\d }[0]{\,d} \newcommand {\l }[0]{\ell } \newcommand {\ddx }[0]{\frac {d}{\d x}} \newcommand {\zeroOverZero }[0]{\ensuremath {\boldsymbol {\tfrac {0}{0}}}} \newcommand {\inftyOverInfty }[0]{\ensuremath {\boldsymbol {\tfrac {\infty }{\infty }}}} \newcommand {\zeroOverInfty }[0]{\ensuremath {\boldsymbol {\tfrac {0}{\infty }}}} \newcommand {\zeroTimesInfty }[0]{\ensuremath {\small \boldsymbol {0\cdot \infty }}} \newcommand {\inftyMinusInfty }[0]{\ensuremath {\small \boldsymbol {\infty -\infty }}} \newcommand {\oneToInfty }[0]{\ensuremath {\boldsymbol {1^\infty }}} \newcommand {\zeroToZero }[0]{\ensuremath {\boldsymbol {0^0}}} \newcommand {\inftyToZero }[0]{\ensuremath {\boldsymbol {\infty ^0}}} \newcommand {\numOverZero }[0]{\ensuremath {\boldsymbol {\tfrac {\#}{0}}}} \newcommand {\dfn }[0]{\textbf } \newcommand {\unit }[0]{\mathop {}\!\mathrm } \newcommand {\eval }[1]{\bigg [ #1 \bigg ]} \newcommand {\seq }[1]{\left ( #1 \right )} \newcommand {\epsilon }[0]{\varepsilon } \newcommand {\phi }[0]{\varphi } \newcommand {\iff }[0]{\Leftrightarrow } \DeclareMathOperator {\arccot }{arccot} \DeclareMathOperator {\arcsec }{arcsec} \DeclareMathOperator {\arccsc }{arccsc} \DeclareMathOperator {\si }{Si} \DeclareMathOperator {\scal }{scal} \DeclareMathOperator {\sign }{sign} \newcommand {\arrowvec }[1]{{\overset {\rightharpoonup }{#1}}} \newcommand {\vec }[1]{{\overset {\boldsymbol {\rightharpoonup }}{\mathbf {#1}}}\hspace {0in}} \newcommand {\point }[1]{\left (#1\right )} \newcommand {\pt }[1]{\mathbf {#1}} \newcommand {\Lim }[2]{\lim _{\point {#1} \to \point {#2}}} \DeclareMathOperator {\proj }{\mathbf {proj}} \newcommand {\veci }[0]{{\boldsymbol {\hat {\imath }}}} \newcommand {\vecj }[0]{{\boldsymbol {\hat {\jmath }}}} \newcommand {\veck }[0]{{\boldsymbol {\hat {k}}}} \newcommand {\vecl }[0]{\vec {\boldsymbol {\l }}} \newcommand {\uvec }[1]{\mathbf {\hat {#1}}} \newcommand {\utan }[0]{\mathbf {\hat {t}}} \newcommand {\unormal }[0]{\mathbf {\hat {n}}} \newcommand {\ubinormal }[0]{\mathbf {\hat {b}}} \newcommand {\dotp }[0]{\bullet } \newcommand {\cross }[0]{\boldsymbol \times } \newcommand {\grad }[0]{\boldsymbol \nabla } \newcommand {\divergence }[0]{\grad \dotp } \newcommand {\curl }[0]{\grad \cross } \newcommand {\lto }[0]{\mathop {\longrightarrow \,}\limits } \newcommand {\bar }[0]{\overline } \newcommand {\surfaceColor }[0]{violet} \newcommand {\surfaceColorTwo }[0]{redyellow} \newcommand {\sliceColor }[0]{greenyellow} \newcommand {\vector }[1]{\left \langle #1\right \rangle } \newcommand {\sectionOutcomes }[0]{} \newcommand {\HyperFirstAtBeginDocument }[0]{\AtBeginDocument }$

Separable differential equations are those in which the dependent and independent variables can be separated on opposite sides of the equation.

In this section we will see that the phrase

Divide and Conquer

is literally true (with the mathematical definition of “divide”) in the context of differential equations. Rather than talk about math, let’s just show you what we’re getting at.

Solve the differential equation $y'(x) = x\cdot y(x)^2$ where $y(0)= 5$ .

We’re going to divide and conquer, literally. Divide both sides of $y'(x) = x\cdot y(x)^2$ by $y(x)^2$ to find $\frac {y'(x)}{y(x)^2} = \answer [given]{x}.$ Now integrate both sides with respect to $x$ $\begin{align*} \int \frac{y'(x)}{y(x)^2} \d x &= \int x \d x,\\ \frac{-1}{y(x)} &= \answer[given]{\frac{x^2}{2}}+C \end{align*}$ and solve for $y(x)$ to find $y(x) = \answer [given]{\frac {-1}{\frac {x^2}{2}+C}}.$ Since $y(0) =5$ we write $\begin{align*} 5&= y(0)\\ &=\frac{-1}{C} \end{align*}$ So $C= -1/5$ . Our final answer is $y(x) = \frac {-1}{\frac {x^2}{2} \answer [given]{-\frac {1}{5}}}.$

Now that we’ve seen an example, let’s discuss this in general.

A separable differential equation is a differential equation which can be written in the form $G(y) y' = H(x)$ In other words, the independent variable $x$ and the function $y$ can be placed on separate sides of the equals sign.

Calculus has lots of different notation. Each of these expressions means exactly the same thing:

$y' = r\cdot y(1-y)$
$y'(x) = r\cdot y(x)(1-y(x))$
$\dd [y]{x} = r \cdot y(1-y)$
$\dd {x} y(x) = r\cdot y(x) (1-y(x))$

We will work with all the notation above.

Which of the following are separable differential equations? Select all that apply.

$y' = y/x$ $y' = \sqrt {x^2+y^2}$ $\cos (y(x))y'(x) = \frac {x}{1+y(x)}$ $\dd [y]{x} = y+5$ $y'(x) = \sin (x+y(x))$ $\dd [y]{x} + x^2 y= x$ $\dd {x} y(x) = \sin (x) y(x) - \cos (x)$ $\dd {x} y(x) = -5 y(x)^{1/3}$

To be a separable differential equation, we must collect $y$ and $y'$ on one side, and $x$ on the other. Typically, you will do this with multiplication or division.

The reason we care about separable differential equations is that:

Separable differential equations help model many real-world contexts.
Separable differential equations are solvable by humans.

The basic ideas is if $G(y)y' = H(x)$ then we can integrate both sides, writing

$\begin{align*} G(y)y' &= H(x)\\ \int G(y) y' \d x &= \int H(x) \d x\\ \int G(y) \d y &= \int H(x) \d x \end{align*}$ If we can symbolically compute these integrals, then we can solve for $y$ . It is now time to work some examples.

Proportional reasoning

In this example, we will show the power of knowing that one quantity is proportional to another.

A snowball has a radius of $2$ inches. After $1.5$ hours, it has a radius of $1$ inch. Assume

the radius changes at rate proportional to the snowball’s surface area, and
the surface area is proportional to the square of the snowball’s radius.

When is the radius of the snowball $0.1$ inches?

We’ll do this by setting up a differential equation. Let $r$ be the radius of the snowball. From the statement of the problem we have $\dd [r]{t} = k_1 \cdot \text {(surface area)}$ and $\text {(surface area)} = k_2 \cdot r^2,$ hence $\dd [r]{t} = k_1 \cdot k_2 \cdot r^2.$ Letting $K = k_1\cdot k_2$ we now can write $\dd [r]{t} = \answer [given]{K}\cdot r^2.$ We now have a differential equation. Let’s divide and conquer. Write with me $\answer [given]{\frac {1}{r^2}}\dd [r]{t} = K$ integrating both sides with respect to $t$ we find $\begin{align*} \int \answer[given]{\frac{1}{r^2}}\dd[r]{t}\d t &= \int K \d t\\ \int \answer[given]{\frac{1}{r^2}}\d r &= \int K \d t\\ \answer[given]{\frac{-1}{r}} &= K\cdot t+C\\ \frac{1}{r} &= -K\cdot t- C\\ r &=\answer[given]{\frac{1}{-K\cdot t- C}}. \end{align*}$ We know that $r(0) = \answer [given]{2}$ , so let’s put that information to use $\begin{align*} \answer[given]{2} &= r(0)\\ \answer[given]{2} &= \frac{1}{-K\cdot 0- C}\\ \answer[given]{2} &= \frac{1}{- C}\\ C &= \answer[given]{-1/2}. \end{align*}$ We also know that $r(1.5) = \answer [given]{1}$ . Again write with me, noting that $1.5 = 3/2$ $\begin{align*} \answer[given]{1} &= r(3/2)\\ \answer[given]{1} &= \frac{1}{-K\cdot 3/2 +1/2}\\ -K\cdot 3/2 +1/2 &=1\\ -K\cdot 3 +1 &=2\\ -K\cdot 3 &=1\\ K &=\answer[given]{-1/3}. \end{align*}$ Hence $r(t) = \answer [given]{\frac {1}{t/3+1/2}}.$ To find when the radius is $0.1$ inches, solve $0.1 = \frac {1}{t/3+1/2}$ and we find $t$ is $\answer [given]{28.5}$ hours.

Some things to note about our last example:

It was solved without using the formula for the surface area of a sphere. For us, it was sufficient to know that the surface area of a sphere is proportional to the square of the radius.
Our solution, $r(t) = \frac {1}{t/3+1/2}$ can be plotted in the slope field determined by $r'=(-1/3)r^2$ :

Which of the following are equilibrium solutions to the differential equation $r'=(-1/3)r^2$ ?

$r = 0$ $r = 1$ $r = t$ $t = 0$ $t = 1$

Exponential and logistic growth

In the science fiction television series Star Trek, a tribble is an alien species that is furry, spherical (radius $\approx 5$ inches), that essentially does nothing but eat and reproduce.

Mr. Spock claims that there are $1771561$ tribbles aboard the space station K7, “assuming one tribble, multiplying with an average litter of ten, producing a new generation every twelve hours over a period of three days.” Explain Mr. Spock’s computation.

Perhaps Mr. Spock used an exponential differential equation to model this. Letting $P$ be the population of tribbles at time $t$ (in hours), it makes sense that the rate that the population is growing is proportional to the size of the population; this is the definition of an exponential model. Write with me $P' = k P.$ Ah, this is a separable differential equation. Let’s solve it: $\begin{align*} P' &= k P\\ \frac{P'}{P} &= k\\ \int \frac{P'}{P} \d t &= \int k\d t\\ \answer[given]{\ln(P)} &= k\cdot t + C. \end{align*}$ However we also know that $P(0) = \answer [given]{1}$ . Write with me $\begin{align*} \ln(P(0)) &= k\cdot 0 + C\\ \ln(\answer[given]{1}) &= C\\ \answer[given]{0} &=C. \end{align*}$ After $12$ hours, we should have $\answer [given]{11}$ tribbles. So $P(12) = \answer [given]{11},$ and we may write $\begin{align*} \ln(P(12)) &= k\cdot \answer[given]{12}\\ \ln(\answer[given]{11}) &= k\cdot \answer[given]{12}\\ \answer[given]{\frac{\ln(11)}{12}} &= k. \end{align*}$ Hence we can model the population growth with $\begin{align*} \ln(P(t)) &= \frac{\ln(11)}{12}\cdot t\\ P(t) &= e^{\frac{\ln(11)}{12} \cdot t}\\ P(t) &= \answer[given]{11}^{t/12}. \end{align*}$ Checking this with Mr. Spock’s computation, we see that $\begin{align*} P(72) &= 11^{72/12}\\ &= 11^6\\ &=1771561. \end{align*}$ just as Mr. Spock stated.

If we graph the solution to the differential equation representing the population of tribbles that we found above

we see that some time after $72$ hours, the growth rate of the tribbles explodes. In reality this would be an ecological disaster. Exponential growth can be scary.

Now we will see a model of population growth with environmental limitations. Suppose that the birthrate of the tribbles is limited by how much food and space is available. This gets us to the idea of carrying capacity.

The carrying capacity of a biological system is the maximum population that can be sustained indefinitely with the given resources.

Let’s work an example involving this concept:

Suppose that the carrying capacity of tribbles aboard a space station is $1000$ tribbles, we start with $1$ tribble, and that growth rate (per tribble) is $10$ tribbles every twelve hours. Model this population with a differential equation.

First, note that the population of tribbles is proportional to the product $P$ and $1000-P$ . This makes sense as when $P$ is small, the growth rate should be small; and when $1000-P$ is small, the growth rate should also be small. Write with me $\dd [P]{t} = k \cdot P (1000-P)$ Ah! This is a separable differential equation. Let’s solve it: $\begin{align*} \dd[P]{t} &= k \cdot P (1000-P)\\ \frac{1}{P (1000-P)}\dd[P]{t} &= k \\ \int \frac{1}{P (1000-P)}\dd[P]{t}\d t &= \int k\d t \\ \int \answer[given]{\frac{1}{P (1000-P)}} \d P &= \int k\d t. \end{align*}$ Use partial fraction decomposition to compute the antiderivative: $\frac {1}{1000}\ln \left (\answer [given]{\frac {P}{1000-P}}\right ) = k\cdot t + C.$ However we also know that $P(0) = \answer [given]{1}$ . Write with me $\begin{align*} \frac{1}{1000}\ln\left(\frac{P(0)}{1000-P(0)}\right) &= k\cdot 0 + C\\ \frac{1}{1000}\ln\left(\frac{1}{999}\right) &= C\\ \answer[given]{\frac{-\ln(999)}{1000}} &= C. \end{align*}$ After $12$ hours, we should have $\answer [given]{11}$ tribbles. So $P(12) = \answer [given]{11},$ and we may write $\begin{align*} \frac{1}{1000}\ln\left(\frac{P(12)}{1000-P(12)}\right) &= k\cdot 12-\frac{\ln(999)}{1000}\\ \ln\left(\frac{11}{989}\right) &= k\cdot 12000-\ln(999) \\ \frac{\ln(10989)-\ln(989)}{12000} &= k. \end{align*}$ Hence we can model the population growth with $\frac {1}{1000}\ln \left (\frac {P(t)}{1000-P(t)}\right ) = \frac {\ln (10989)-\ln (989)}{12000}\cdot t - \frac {\ln (999)}{1000}$ $\begin{align*} \ln\left(\frac{P(t)}{1000-P(t)}\right) &= \frac{\ln(10989)-\ln(989)}{12}\cdot t - \ln(999)\\ \frac{P(t)}{1000-P(t)} &= \left(\frac{10989}{989}\right)^{t/12}\cdot \frac{1}{999} \end{align*}$ Solving for $P(t)$ we find $P(t) = \frac {\left (\frac {10989}{989}\right )^{t/12}\cdot \frac {1000}{999}}{1+\left (\frac {10989}{989}\right )^{t/12}\cdot \frac {1}{999}}.$ If we graph this solution, we see a nice “S-curve”

In our last example we used the differential equation $P' = k\cdot P(A-P)$ This is called the logistic differential equation where $k$ and $A$ are constants. Let’s examine a slope field for this model with some reasonable values of $k$ and $A$ :

Which of the following are equilibrium solutions to the differential equation above?

$P = 0$ $P = 1$ $P = 2$ $P = t$

Predator-prey model

In nature there are animals that are predators and animals that are prey. In the early 20th Century, Lotka and Volterra suggested the following model to help us understand populations of predators and prey. Let

$\begin{align*} x(t) &= \text{population of predators at time $t$}\\ y(t) &= \text{population of prey at time $t$}. \end{align*}$ Lotka and Volterra made the following assumptions:

If there are no prey, then the predators starve at some rate proportional to the number of predators. If there are prey, then the predator’s population grows at a rate proportional to the product of the population of the predators and the prey. In the language of calculus we write $\dd [x]{t} = -p x + q x y.$
If there are no predators, then the prey’s population grows at some rate proportional to the number of prey. If there are predators, then the prey’s population is reduced at a rate proportional to the product of the population of the predators and the prey. In the language of calculus we write $\dd [y]{t} = r y - s xy.$

Taking these two equations, and dividing them we find

$\begin{align*} \dd[y]{x} &= \frac{\dd[y]{t}}{\dd[x]{t}}\\ &=\frac{r y - s xy}{-p x + q x y}\\ &=\frac{y(r-sx)}{x(qy-p)}. \end{align*}$ In the previous examples, it worked out that we could solve for $y$ as a function of $x$ . As we will see with the predator-prey model, this is not always the case.

Find the solution to the differential equation $y' = \frac {y(4-2x)}{x(y-3)}$ that passes through the point $(1,e)$ .

To start, note that this is a separable differential equation. Write with me: $\begin{align*} y' &= \frac{y(4-2x)}{x(y-3)}\\ \frac{(y-3)y'}{3}&=\frac{(4-2x)}{x}\\ \int \frac{(y-3)y'}{y}\d x&=\int \frac{(4-2x)}{x}\d x\\ \int y'-\frac{3y'}{y}\d x&=\int \frac{4}{x}-2\d x\\ \answer[given]{y-3\ln(y)} &=4\ln(x)-2x+C. \end{align*}$ Remembering that the solution must pass through point $(1,e)$ , we write $\begin{align*} e-3\ln(e) &=4\ln(1)-2+C\\ \answer[given]{e-3} &=-2+C\\ \answer[given]{e-1} &=C. \end{align*}$ Hence we have the following implicit equation relating $x$ and $y$ : $y-3\ln (y) =4\ln (x)-2x+e-1.$ Let’s look at this solution along with the slope field:

From this graph, we can see the cyclic nature of the predator-prey relationship, where the population of the predators is on the $x$ -axis, and the population of the prey is on the $y$ -axis.