Model predator-prey population dynamics using the Lotka-Volterra equations. Simulate oscillating populations, explore equilibria, and visualize classic ecological cycles.
The Lotka-Volterra equations are the foundational mathematical model of predator-prey dynamics in ecology. Published independently by Alfred Lotka (1925) and Vito Volterra (1926), these coupled differential equations describe how two interacting populations—a prey species and its predator—oscillate in perpetual, out-of-phase cycles. When prey are abundant, predators thrive and multiply; their growing numbers then deplete the prey; with fewer prey, predators decline; reduced predation allows prey to recover; and the cycle repeats.
The classic equations are: dN/dt = αN - βNP (prey growth minus predation) and dP/dt = δNP - γP (predator reproduction from consumption minus natural death). The four parameters—prey growth rate (α), predation rate (β), predator reproduction efficiency (δ), and predator death rate (γ)—determine the amplitude, period, and equilibrium of the oscillations.
This calculator numerically solves the Lotka-Volterra system using Euler's method, letting you explore how parameter changes affect population dynamics. Model classic examples like lynx-hare cycles, wolf-deer interactions, or abstract ecological scenarios. Visualize the time series, identify equilibrium populations, and understand the sensitivity of ecosystems to parameter perturbations.
The Lotka-Volterra equations are a cornerstone of mathematical ecology. Understanding predator-prey dynamics is essential for conservation biology, wildlife management, pest control, and understanding ecosystem stability. This calculator makes the math accessible and visual. Keep these notes focused on your operational context. Tie the context to the calculator’s intended domain. Use this clarification to avoid ambiguous interpretation.
dN/dt = αN − βNP (prey), dP/dt = δNP − γP (predator). Equilibrium: N* = γ/δ, P* = α/β. Period ≈ 2π/√(αγ). Parameters: α = prey intrinsic growth rate, β = predation rate coefficient, δ = predator reproduction efficiency per prey consumed, γ = predator natural death rate. Solved via Euler method: N(t+dt) = N(t) + dt × (αN − βNP).
Result: Oscillating populations: prey peaks ~120, predator peaks ~25
With prey growth α=0.1, predation β=0.01, predator efficiency δ=0.005, and predator death γ=0.1, the equilibrium is N*=γ/δ=20, P*=α/β=10. Starting at N=100, P=20 (far from equilibrium), both populations oscillate with ~60-step period. Prey peaks about 2.5× higher than predators due to higher turnover.
The classic Lotka-Volterra system has a remarkable mathematical property: it conserves a quantity V = δN - γ ln(N) + βP - α ln(P). This conservation means orbits in the N-P phase plane are closed loops—populations cycle forever without damping or amplification. This is analogous to conservation of energy in physics and means the model is "structurally unstable": any modification (adding carrying capacity, functional responses, or stochasticity) changes the qualitative behavior.
The period of oscillation is approximately 2π/√(αγ), depending on the prey growth rate and predator death rate. The amplitude depends on the starting conditions—orbits further from equilibrium have larger oscillations. Unlike a simple pendulum, the oscillations are not sinusoidal; prey peaks are sharper and predator peaks are broader.
Ecologists have developed numerous extensions to the basic model. The Rosenzweig-MacArthur model adds logistic prey growth (carrying capacity K) and a Type II functional response (predator satiation). This produces damped oscillations converging to a stable equilibrium at low K, but paradoxically destabilizes as K increases (the "paradox of enrichment"—adding resources can cause population crashes).
The competitive Lotka-Volterra equations model two competing species rather than predator-prey, leading to outcomes including competitive exclusion (one species wins), stable coexistence, or bistability (outcome depends on initial conditions). These models underpin the competitive exclusion principle—two species competing for the same niche cannot stably coexist.
Modern conservation biology uses Lotka-Volterra-derived models extensively. Wolf reintroduction in Yellowstone, deer population management, marine predator-prey dynamics (sharks and prey fish), and invasive species control all rely on understanding predator-prey dynamics. The Lotka-Volterra framework helps predict: How will removing a predator affect prey populations? What happens when prey habitat is reduced? How does introducing a new predator reshape the ecosystem?
They model the simplest predator-prey interaction: prey grow exponentially in the absence of predators, predation is proportional to encounters between prey and predators (proportional to both populations), predators reproduce in proportion to prey consumed, and predators die at a constant per-capita rate. This produces characteristic oscillating population cycles.
The non-trivial equilibrium (where both species coexist) is: N* = γ/δ (prey equilibrium = predator death rate / predator efficiency) and P* = α/β (predator equilibrium = prey growth rate / predation rate). Counterintuitively, the prey equilibrium depends only on predator parameters, and vice versa.
The oscillation arises from the time lag in predator response to prey changes. When prey are abundant, predators reproduce—but by the time predator numbers peak, prey have been depleted. Predators then starve and decline, but by the time they've declined enough, prey have already begun recovering. This delayed feedback creates persistent cycles.
The basic Lotka-Volterra model is highly simplified—it assumes unlimited prey resources, no spatial structure, no age structure, no other species, and perfectly linear functional responses. Real ecosystems are far more complex. However, the model captures the essential mechanism of predator-prey cycles seen in nature (e.g., Canadian lynx-hare cycles) and serves as the foundation for more sophisticated ecological models.
The most famous real-world example of Lotka-Volterra dynamics is the snowshoe hare and Canadian lynx cycle, documented through Hudson's Bay Company fur trading records spanning 200+ years. The hare population peaks every ~10 years, with the lynx population peaking 1-2 years later—exactly the phase-shifted oscillation the Lotka-Volterra equations predict.
If a parameter changes (e.g., prey growth rate increases due to habitat improvement), the equilibrium shifts and populations transition to new oscillation patterns. If predator death rate increases (e.g., disease), prey equilibrium rises (N* = γ/δ) and the system oscillates around a new point. The model predicts that pesticides targeting prey paradoxically reduce predator populations more.