[Predator & Prey][Non-linear Dynamics]Feb 13

The Calculus of Starvation

Coupled dynamics of population oscillations and the mathematics of survival.

The Calculus of Starvation

Act 1: The Wealth Paradox (c. 1798)

While the Industrial Revolution injected unprecedented economic energy into Britain, the benefits were remarkably lopsided. England was producing more than ever—thanks to the Agricultural Revolution boosting domestic crop yields and new steam engines supercharging textile mills—yet the average person felt little improvement.

Economists call this phenomenon "Engels' Pause": a 50-year period when the nation grew wealthy, yet workers remained poor. The wealth flowed elsewhere—primarily into building more factories.

Does this remind you of something? More and more money being poured into making Data-Centers? The "AI pause", contemporary talk in 2026.

As cities crowded and smokestacks multiplied, a darker anxiety began to fester. The concern transcended mere poverty; it touched on physics itself. Many feared this growth was a temporary illusion, and that nature would inevitably slam on the brakes.

One of the first thinkers to mathematize this anxiety was Thomas Malthus.[1]

Thomas Robert Malthus

Thomas Robert Malthus (1766–1834). Wellcome Collection. Public domain.
Source: Wikimedia Commons

Malthus argued that when resources become abundant, populations increase faster than living standards improve. This creates a self-limiting cycle: rising populations lead to resource scarcity, reduced per-capita availability, and eventual population decline—only to repeat when conditions improve.

Mathematical Thought Experiment: Let us formalize these intuitions. The following derivations are exercises in understanding, not endorsements of any particular worldview.

Act 2: Exponential Growth Under Unlimited Resources

Consider a simplified scenario with two key assumptions:

Assumption A1: Resources are unlimited (R=R = \infty). A population NN converts resources into offspring with constant efficiency. The rate of population change is proportional to the current population—the more individuals present, the more births occur:

dNdt=rN(1)\frac{dN}{dt} = rN \tag{1}

where rr is a positive constant (intrinsic growth rate).

To solve this, we integrate:

N0Nt1NdN=0trdtln(Nt)ln(N0)=rtln(N(t)N0)=rtN(t)N0=ertN(t)=N0ert\begin{align} \int_{N_0}^{N_t} \frac{1}{N} \, dN &= \int_{0}^{t} r \, dt \tag{2} \\ \ln(N_t) - \ln(N_0) &= rt \tag{3} \\ \ln\left( \frac{N(t)}{N_0} \right) &= rt \tag{4} \\ \frac{N(t)}{N_0} &= e^{rt} \tag{5} \\ N(t) &= N_0 e^{rt} \tag{6} \end{align}

where N0N_0 is the initial population.

Assumption A2: Only existing organisms can produce offspring (no spontaneous generation).

Exponential Function plot

Exponential Function used to model population growth based on assumption A1 and A2.

Exponential Function plot

Same exponential model compared to observed world population (1750–2026):

Caveats: The exponential model provides limited predictive value for real-world populations. While the fit captures the general upward trend, it diverges significantly in later periods. Yet it validates our core inference: any population with abundant resources and unchecked reproduction follows exponential trajectories. Understanding this baseline is essential before introducing constraints.

Act 3: Carrying Capacity—The Saturation Limit

The exponential model captures initial growth patterns but ultimately fails because it ignores environmental constraints. If Malthus was the pessimist who saw the cliff edge ahead, Pierre François Verhulst was the engineer attempting to predict when the upward trajectory would bend back down.[2]

Thought Experiment: A bacterial colony in a petri dish does not conquer the universe. Bacteria multiply exponentially until they exhaust nutrients, then population stabilizes or crashes. Verhulst recognized that the growth term (rNrN) wasn't wrong—merely incomplete. The missing piece was a negative feedback mechanism reflecting environmental resistance.

Verhulst introduced the Logistic Equation—a modified growth model incorporating a saturation term.[3] This equation is named after its characteristic S-shaped (sigmoid) curve.

dNdt=rN(1NK)(7)\frac{dN}{dt} = rN \left( 1 - \frac{N}{K} \right) \tag{7}

where:

Physical Interpretation:

When NKN \ll K (low density), the feedback ≈ 1, and growth is nearly exponential—the system has room to expand.

When NKN \to K (high density), the feedback → 0, and growth halts. The population approaches KK asymptotically.

This single variable, KK (carrying capacity), fundamentally changed how we understand natural populations. It revealed that populations are not controlled by abstract "fate," but by concrete density-dependent resource limitations—a principle confirmed across ecology.[4]

Act 4: Interactive Exploration—The Logistic Growth Simulator

Adjust the parameters below to observe how feedback reshapes population dynamics. Notice when the logistic and exponential curves diverge, and why.

[Logistic Growth Model]

dN/dt = rN(1 − N/K) = 0.50 × N × (1 − N/1000)
0.1 (slow)2.0 (fast)
1005000
1500
5s50s
Physical Interpretation:
  • G-term (rN): Growth engine — "Multiply!"
  • B-term (1 − N/K): Brake — slows growth as N approaches K
  • When N ≪ K: Growth is nearly exponential (brake ≈ 1)
  • When N → K: Growth halts (brake → 0)

Reflection Question: At what population relative to KK does exponential growth become visibly constrained?

Test your understanding with the following questions:

[Diagnostic Check: Logistic Growth]

1. If you increase the rate of growth (r) in the logistic model, does the population grow beyond the carrying capacity (K)?
2. For r=0.5, K = 100, N0 = 1, up to how many seconds do the exponential growth and logistic functions behave similarly?

References

  1. Malthus, T. R. (1798). An Essay on the Principle of Population - The foundational text on population dynamics and resource scarcity. Available via Internet Archive.
  2. Verhulst, P.-F. (1845). Recherches mathématiques sur la loi d'accroissement de la population - Original derivation of the logistic equation. Nouveaux Mémoires de l'Académie Royale des Sciences et Belles-Lettres de Bruxelles.
  3. Wikipedia: Logistic function - Overview of the logistic curve and its applications in biology, demography, and epidemiology.
  4. Gotelli, N. J. (2001). A Primer of Ecology (3rd ed.) - Comprehensive treatment of carrying capacity and density-dependent population regulation in ecological systems.
  5. United Nations: World Population Data - Official global population statistics and demographic trends used for curve fitting examples.
The Price of Speed: Runge-Kutta vs EulerALL CHAPTERSThe Couple-Dance