Predator-Prey Models | Varsity Tutors

SECTION 1

Historical Context & Motivation

For centuries, farmers and naturalists noticed a puzzling pattern: when rabbit populations boomed, fox populations soon followed — but then the rabbits crashed, and the foxes crashed shortly after. These cycles seemed to repeat endlessly. Before the twentieth century, nobody had a mathematical way to explain why predator and prey populations rise and fall together in such a predictable rhythm.

The breakthrough came when two mathematicians, working independently and on different continents, each discovered the same pair of equations that could model this behavior. Their work launched an entire field of mathematical ecology — using math to understand how species interact. The story of these discoveries is tied to real-world data, from fish catches in the Adriatic Sea to fur-trapping records in Canada.

1900s

Hudson's Bay Fur Records

The Hudson's Bay Company in Canada kept detailed records of lynx and snowshoe hare pelts traded each year. These records revealed striking population cycles — roughly every 10 years, hare numbers peaked, followed by a peak in lynx numbers.

1925

Alfred Lotka's Chemical Analogy

American mathematician Alfred J. Lotka published a pair of differential equations originally inspired by chemical reactions, but he recognized they perfectly described predator-prey dynamics.

1926

Vito Volterra's Fish Market Insight

Italian mathematician Vito Volterra independently derived the same equations after his son-in-law, a marine biologist, asked why predatory fish populations increased in the Adriatic Sea during World War I, when fishing slowed down.

1930s

Ecological Experiments

Biologist Georgii Gause tested predator-prey models in the laboratory using microscopic organisms (Paramecium and Didinium), confirming the oscillatory behavior the equations predicted.

2000s

Modern Computational Models

With powerful computers, ecologists now simulate complex predator-prey systems involving multiple species, spatial movement, and environmental change — all built on the foundation that Lotka and Volterra laid.

The central question these pioneers tackled was deceptively simple: Can we write a pair of equations that explains why predator and prey populations oscillate over time? The answer turned out to be yes, and the resulting model — known as the Lotka–Volterra equations — remains one of the most important tools in ecology and differential equations.

SECTION 2

Core Principles & Definitions

Before diving into equations, you need to understand the key ideas that make predator-prey models work. Every Lotka–Volterra model is built on a set of simple biological assumptions — each one easy to understand on its own, but together they produce complex and fascinating behavior.

Exponential Prey Growth

Without predators, the prey population (like rabbits) grows exponentially — the more rabbits there are, the faster the population increases. This is described by a positive growth rate constant.

Predator Decline Without Food

Without prey, the predator population (like foxes) declines exponentially because they starve. This is described by a negative death rate constant.

Encounters Drive Interaction

The rate at which predators catch prey depends on how often they encounter each other. More prey and more predators means more encounters. Mathematically, this is modeled by multiplying the two population sizes together.

Coupled System

The prey equation depends on the predator population, and the predator equation depends on the prey population. These equations are coupled — you cannot solve one without the other. This linkage creates the oscillations.

Closed Ecosystem

The basic model assumes a closed ecosystem — no migration in or out, no other species, and no environmental changes. This simplification lets us focus on the core predator-prey interaction.

✦ KEY TAKEAWAY

Think of predator-prey dynamics like a thermostat in your house. When the temperature drops (prey population rises), the heater kicks on (predators thrive and multiply). When the temperature rises too high (too many predators eat the prey), the heater shuts off (predators decline from lack of food). The system naturally oscillates around a balance point, never quite settling at a fixed temperature. The key insight is that predators and prey are locked in a feedback loop — each population's fate depends on the other.

SECTION 3

Visualizing Population Oscillations

The most striking feature of predator-prey models is the way the two populations cycle over time. The diagram below shows a typical output of the Lotka–Volterra equations: the prey population (shown in cyan) rises first, followed by the predator population (shown in pink). When predators become too numerous, prey decline sharply, and predators soon follow.

Notice how the prey curve (cyan) peaks before the predator curve (pink). The predator population lags behind because predators need time to reproduce after the food supply increases. This time delay is what creates the oscillation.

The diagram reveals the signature pattern of predator-prey dynamics. First, prey are abundant and predators are few, so the prey population grows rapidly. As predators feast on the plentiful prey, the predator population grows. Eventually predators become so numerous that prey are caught faster than they can reproduce, causing the prey to decline. With less food available, predators then decline too — and the whole cycle starts again. This quarter-cycle lag between the two curves is the hallmark of a predator-prey system.

SECTION 4

The Lotka–Volterra Equations

Now that you have a visual intuition, let's translate those ideas into equations. The Lotka–Volterra model is a system of two differential equations — one for the prey population and one for the predator population. Each equation describes the rate of change of a population over time. If you've seen expressions like dx/dt before, you're ready. If not, just think of dx/dt as "how fast x is changing at this moment."

PREY EQUATION

dx/dt = αx − βxy

x = prey population, y = predator population, α (alpha) = prey birth rate, β (beta) = rate at which predators catch prey. The term αx represents natural growth, and −βxy represents deaths from predation.

PREDATOR EQUATION

dy/dt = δxy − γy

δ (delta) = rate at which predators convert food into offspring, γ (gamma) = predator death rate. The term δxy represents predator births (which depend on eating prey), and −γy represents natural deaths.

Notice the xy terms that appear in both equations. This product is the mathematical translation of "encounters depend on both populations." If there are 100 rabbits and 10 foxes, the encounter rate is proportional to 100 × 10 = 1000. If rabbits double to 200, encounters double to 2000. The xy term acts like a bridge that couples the two equations together.

EQUILIBRIUM POINTS

x* = γ/δ, y* = α/β

Setting both dx/dt = 0 and dy/dt = 0 gives us the equilibrium — the population values where neither species is growing nor shrinking. There is a trivial equilibrium at (0, 0) and a non-trivial one at (γ/δ, α/β).

💡 Why Can't We Solve These Exactly?

Because of the xy terms, these equations are nonlinear — meaning we generally cannot write a neat formula for x(t) and y(t). Instead, we use numerical methods (computers stepping through tiny time intervals) or analyze the behavior qualitatively using phase-plane diagrams, which you'll see in the next section.

SECTION 5

The Phase-Plane Portrait

Instead of plotting each population against time separately, ecologists and mathematicians often plot prey vs. predator directly. This is called a phase-plane diagram. In this view, the x-axis shows the prey population and the y-axis shows the predator population. As time passes, the system traces out a closed loop — a cycle that the populations travel around endlessly.

Each closed loop represents a possible orbit of the system, depending on initial populations. The green dot marks the equilibrium point. Notice that the system travels counterclockwise around the loop: prey increase first (move right), then predators increase (move up), then prey decrease (move left), then predators decrease (move down).

The phase-plane diagram reveals something remarkable about the basic Lotka–Volterra model: the orbits are closed loops, meaning the populations cycle forever without ever settling at the equilibrium. Different starting conditions produce different loops — start with more prey, and you get a wider orbit. In the real world, random events and other factors cause orbits to shift over time, but this idealized view captures the essential cycling behavior perfectly.

SECTION 6

Worked Example: Rabbits and Foxes

Let's apply the Lotka–Volterra equations to a specific scenario. Suppose we have a meadow with rabbits (prey) and foxes (predators), and we know the following parameter values:

Finding Equilibrium and Interpreting Rates of Change

Step 1 — State the Given ParametersWe are told: α = 0.4 (rabbits reproduce at rate 0.4 per month), β = 0.01 (predation rate), δ = 0.005 (conversion rate of prey into predator offspring), γ = 0.3 (foxes die at rate 0.3 per month). Initial populations: x₀ = 100 rabbits, y₀ = 20 foxes.

α = 0.4, β = 0.01, δ = 0.005, γ = 0.3, x₀ = 100, y₀ = 20

Step 2 — Find the Non-Trivial EquilibriumAt equilibrium, both dx/dt = 0 and dy/dt = 0. Using our formulas: x* = γ/δ = 0.3/0.005 = 60. And y* = α/β = 0.4/0.01 = 40. This means the system would be perfectly balanced with 60 rabbits and 40 foxes.

Equilibrium: x* = 60 rabbits, y* = 40 foxes

Step 3 — Calculate the Initial Rate of Change for PreyAt time t = 0, we plug x₀ = 100 and y₀ = 20 into the prey equation: dx/dt = αx − βxy = 0.4(100) − 0.01(100)(20) = 40 − 20 = 20. This means at this moment, the rabbit population is increasing at a rate of 20 rabbits per month.

dx/dt = 20 rabbits per month (prey increasing)

Step 4 — Calculate the Initial Rate of Change for PredatorsSimilarly, we plug into the predator equation: dy/dt = δxy − γy = 0.005(100)(20) − 0.3(20) = 10 − 6 = 4. The fox population is increasing at 4 foxes per month.

dy/dt = 4 foxes per month (predators increasing)

Step 5 — Interpret the ResultsSince we started with more prey than the equilibrium (100 > 60) and fewer predators than the equilibrium (20 < 40), both populations are growing. The system is in the lower-right portion of the phase-plane loop, moving up and to the right. As foxes multiply, they'll eventually eat rabbits faster than the rabbits can reproduce, and the cycle will continue.

Both populations grow initially. The system is heading toward peak prey, then peak predators will follow.

SECTION 7

Strengths and Limitations of the Model

The Lotka–Volterra model is elegant and powerful, but like all mathematical models, it simplifies reality. Understanding what it gets right and what it misses is essential for applying it wisely. The table below compares the strengths and limitations side by side.

Comparison of model strengths versus real-world limitations

Strengths	Limitations
Captures the fundamental oscillatory behavior seen in real predator-prey systems	Assumes unlimited prey growth without predators — ignores carrying capacity and limited resources
Requires only four parameters, making it easy to set up and study	Predicts perfectly periodic orbits that never decay — real populations eventually settle or go extinct
Equilibrium analysis provides clear insight into balance between species	Ignores other species, disease, migration, seasonal effects, and random events
Serves as a foundation for all more advanced ecological models	The encounter term (xy) assumes predators are never full — a predator always eats more if it meets more prey

✦ KEY TAKEAWAY

Think of the Lotka–Volterra model like a simplified flight simulator. It won't teach you about turbulence, crosswinds, or engine failure, but it perfectly demonstrates the core principle of how lift and gravity interact to keep a plane in the air. In the same way, this model strips away all the noise of real ecosystems to reveal the fundamental mechanism that makes predator and prey populations cycle. Once you understand this core, you can add complexity layer by layer.

SECTION 8

Connection to Advanced Models

The basic Lotka–Volterra model is a starting point, not an endpoint. Over the decades, mathematicians and ecologists have developed more realistic extensions. These advanced models fix the limitations we just discussed while preserving the core insight of coupled differential equations. Here's how the basic model compares to some important extensions.

How the basic model connects to more advanced population dynamics

Feature	Basic Lotka–Volterra	Advanced Extensions
Prey growth	Exponential (unlimited)	Logistic (with carrying capacity K)
Predator feeding	Linear (always eats more)	Functional response (predators get "full")
Orbit behavior	Closed loops (forever)	Spirals toward stable equilibrium or limit cycles
Number of species	Exactly 2	3+ species food webs
Randomness	None (deterministic)	Stochastic noise added for realism

One of the most common improvements is adding logistic growth to the prey equation. Instead of αx, you write αx(1 − x/K), where K is the carrying capacity — the maximum prey population the environment can sustain. This single change makes the orbits spiral inward toward a stable equilibrium, which is much more realistic. If you continue your study of differential equations in college, you'll encounter these richer models in courses on nonlinear dynamics and mathematical biology.

SECTION 9

Practice Problems

PROBLEM 1 — CONCEPTUAL

In the Lotka–Volterra model, why does the predator population peak after the prey population, rather than at the same time? Explain using biological reasoning.

PROBLEM 2 — BASIC CALCULATION

Given α = 0.5, β = 0.02, δ = 0.01, γ = 0.4, find the non-trivial equilibrium populations x* and y*.

PROBLEM 3 — INTERMEDIATE

Using the same parameters from Problem 2 (α = 0.5, β = 0.02, δ = 0.01, γ = 0.4), suppose the initial populations are x₀ = 80 prey and y₀ = 10 predators. Calculate dx/dt and dy/dt at t = 0. Are both populations increasing, both decreasing, or is one rising while the other falls?

PROBLEM 4 — APPLIED

A wildlife reserve monitors a wolf-deer system. Currently there are 500 deer and 50 wolves, with α = 0.6, β = 0.005, δ = 0.002, and γ = 0.5. A conservation manager proposes introducing 20 more wolves to "control" the deer. Without calculating trajectories, use the equilibrium values to argue whether this will permanently reduce the deer population.

PROBLEM 5 — CRITICAL THINKING

The basic Lotka–Volterra prey equation is dx/dt = αx − βxy. Suppose you modify it to dx/dt = αx(1 − x/K) − βxy, where K is the carrying capacity. Explain in words what the new term (1 − x/K) does when x is small compared to K, when x is close to K, and when x exceeds K. How would this change the shape of the orbits in the phase plane?

SUMMARY

Lesson Summary

The Lotka–Volterra predator-prey model is a system of two coupled differential equations that describes how predator and prey populations change over time. The prey equation (dx/dt = αx − βxy) captures exponential growth reduced by predation, while the predator equation (dy/dt = δxy − γy) captures predator growth from eating prey minus natural death. The xy interaction term links the two equations, creating a feedback loop that produces oscillating populations.

The equilibrium point (x* = γ/δ, y* = α/β) represents the population balance, and the phase-plane diagram reveals that orbits form closed loops around this point — a signature of perpetual oscillation. While the basic model has limitations (no carrying capacity, no randomness), it captures the essential mechanism that drives predator-prey cycles and serves as the foundation for all advanced ecological models. Developed independently by Lotka (1925) and Volterra (1926), these equations remain a cornerstone of mathematical biology.