Category Archives: Logistic Growth

Verhulst…

I have to admit, there was a time when I was in love with an equation.

The class of nonlinear differential equations are fairly notoriously pains in the butt, but for me, they represented something other than that pain. They reflected the artful beauty found in integral calculus. Sergey Lvin, was a soft spoken mathematics professor at the University of Maine who took pride in his teaching. I recall asking him once why my differential equations class was so difficult while derivatives were really automatic and straightforward. He smiled and said to me that “Differentiation was a science and integration was an art.” What he meant was that one could follow the rules by rote and take the derivative of any function imaginable, and yet to find the integral of a function was often a journey through considering various solution tactics. And to make matters worse, differential equations are implicit relationships between derivatives and the state quantities. Rather than \frac{dx}{dt}=f(t) we have that \frac{dx}{dt}=f(x(t),t) and that makes all the difference in the world.

However, utilizing the tools you’ve learned with the experience you’ve gained to approach problems, often from oblique lines of inquiry, is the joy of mathematics. I enjoyed learning methods, building my tools belt to handle any situation, and this is in part why I remember liking the Verhulst equation so much:

(1)   \begin{align*}\frac{dN}{dt}&=rN\left(1-\frac{N}{K}\right)\end{align*}

Pierre Verhulst was a mathematician of some note from the early 19th century. In a collection of papers starting in 1938 (of which I must get my hands on) he proposes that for a population of N(t) individuals at time t that there is some intrinsic growth rate to the population, r, and a natural carrying capacity for it within its environment, K. It took him 7 years to call the solution to this differential equation the logistic equation and some how the wires crossed and that is what I see the differential equation called. If you look close then you can show that for N(t)>0 and N(t)\approx0 we can approximate the logistic growth equation with

(2)   \begin{align*}\frac{dN}{dt}=rN\end{align*}

which is the only differential equation that I know the solution to N(t)=N(0)e^{rt} so mathematically r matches that biological meaning of growth rate. Furthermore, for r>0 the logistic growth equation results in a growing population (dN/dt>0) for all N(t)\in(0,K) gives credence to the interpretation that K is the carrying capacity. This was amazing to me. Here was a fairly complex (for me at the time) differential equation where the parameters just made sense. To boot it was an adventure to solve. Recall that

(3)   \begin{align*}\frac{dN}{dt}=rN\left(1-\frac{N}{K}\right)\leftrightarrow\frac{dN}{N\left(1-\frac{N}{K}\right)}=rdt,\end{align*}

but now the left hand-side is just a mess.  This is where the art comes in.

We note that \frac{1}{N\left(1-\frac{N}{K}\right)}=\frac{K}{N(K-N)}=\frac{1}{N}+\frac{1}{K-N} (when did you learn partial fraction decomposition? I cannot recall when I was supposed to have and Eisso Atzema – my professor for some course requiring said knowledge – liked to tease me for my inability to recall said information. Is it he who I first gained math related snarkiness from?). The integral of each of these parts is trivial at this point, let me be pendantic for anyone not cruising along:

(4)   \begin{align*} \frac{dN}{dt}&=rN\left(1-\frac{N}{K}\right),\\ \frac{dN}{N\left(1-\frac{N}{K}\right)}&=rdt,\\ \frac{K}{N(K-N)}dN&=rdt,\\ \left(\frac{1}{N}+\frac{1}{K-N}\right)dN&=rdt,\\ \int\left(\frac{1}{N}+\frac{1}{K-N}\right)dN&=\int rdt,\\ \int\frac{dN}{N}+\int\frac{dN}{K-N}&=rt+c. \end{align*}

Now use u-substitution:

(5)   \begin{align*} u=K-N&\leftrightarrow du=-dN,\\ \ln|N(t)|-\int\frac{du}{u}&=rt+c,\\ \ln|N(t)|-\ln|u(t)|&=rt+c,\\ \ln\left|\frac{N(t)}{u(t)}\right|&=rt+c,\\ \ln\left|\frac{N(t)}{K-N(t)}\right|&=rt+c,\\ \left|\frac{N(t)}{K-N(t)}\right|&=e^{rt+c}=e^{rt}e^c,\\ \frac{N(t)}{K-N(t)}&=\pm e^ce^{rt},\\ C&=\pm e^c,\\ N(t)&=(K-N(t))Ce^{rt},\\ N(t)\left(1+Ce^{rt}\right)&=KCe^{rt},\\ N(t)&=\frac{KCe^{rt}}{1+Ce^{rt}},\\ N_0=N(0)&=\frac{KC}{1+C},\\ N_0&=C(K-N_0),\\ \frac{N_0}{K-N_0}&=C,\\ N(t)&=\frac{K\frac{N_0}{K-N_0}e^{rt}}{1+\frac{N_0}{K-N_0}e^{rt}},\\ N(t)&=\frac{KN_0e^{rt}}{K+N_0(1-e^{rt})}. \end{align*}

There.  Breathe. I’ll have you know that out of those 21 steps that 15 of them are strictly algebra (just a hair under 75%).

The point, to bring us back to it, is that I was shown an equation where everything has real life meaning and the solution was a fun journey through using multiple tools (separation of variables – a change of differential notation, partial fraction decomposition, linearity of integral operator, integration by substitution, knowledge of transcendentals like \ln and exponentiation, and fancy \pm obfuscation)and a lot of my knowledge. I could sketch the solution and understand what it meant. I could consider the limiting behavior explicitly or qualitatively.  I loved this thing. I loved it. Loved.  Past tense. But the funny thing was that I would not fall out of love for a number of years and in the interim my admiration of this equation would only grow, but that is for another time.