Differential Equations · Lotka–Volterra equations

Lotka–Volterra equations

The Lotka–Volterra equations model a prey and a predator population: the prey grow on their own and are eaten, and the predators grow by eating them and otherwise die out. Their numbers rise and fall in cycles around a balance point.

Lotka–Volterra equations in this visualization

Here xx is the number of prey and yy the number of predators, both in thousands. Each changes at a rate per head: the prey grow at α\alpha and are eaten in proportion to the predators, and the predators die at γ\gamma and grow in proportion to the prey. The phase plane follows the point (x,y)(x, y) from your start, and the time graph below draws both numbers year by year.

The model without fishing or a ceiling
x′=x(α−βy),y′=y(δx−γ)x' = x(\alpha - \beta y), \qquad y' = y(\delta x - \gamma)

The prey's rate α−βy\alpha - \beta y is zero on the line y=α/βy = \alpha/\beta, and the predators' rate δx−γ\delta x - \gamma is zero on the line x=γ/δx = \gamma/\delta. These dashed zero-growth lines cross at the balance point (x∗,y∗)=(γ/δ,α/β)(x^*, y^*) = (\gamma/\delta, \alpha/\beta) and cut the plane into four regions. Below y∗y^* the prey increase, and above it they decrease. To the right of x∗x^* the predators increase, and to the left they decrease. So the path goes round counterclockwise. The predators peak when they stop growing, which is when the prey have fallen back to x∗x^*, after the prey's own peak. That delay is the lag in the readout.

Going round does not by itself bring the path back; a spiral goes round too. What closes it is a quantity that stays the same along every solution:V=δx−γln⁡x+βy−αln⁡y,V′=(δx−γ)(α−βy)+(βy−α)(δx−γ)=0.\begin{aligned} V &= \delta x - \gamma \ln x + \beta y - \alpha \ln y, \\ V' &= (\delta x - \gamma)(\alpha - \beta y) + (\beta y - \alpha)(\delta x - \gamma) = 0. \end{aligned}Each level of VV above its minimum is a closed curve around the balance point, so the path follows it back to its start and repeats. With the ceiling KK, crowding slows the prey's growth and VV is no longer constant; in the ceiling highlight the path spirals in toward the balance point instead of closing.

However big the swings, every full cycle averages exactly the balance point. Dividing the predators' equation by yy gives (ln⁡y)′=δx−γ(\ln y)' = \delta x - \gamma. Over one cycle yy comes back to its start, so this integrates to 0, and the average of xx is γ/δ=x∗\gamma/\delta = x^*. The prey's equation gives the average of yy, α/β=y∗\alpha/\beta = y^*, the same way. Larger swings change the extremes and the period, not the averages.

Fishing both populations continuously at the same rate per head, ε\varepsilon, turns α\alpha into α−ε\alpha - \varepsilon and γ\gamma into γ+ε\gamma + \varepsilon. Without the ceiling the balance point moves to x∗=(γ+ε)/δx^* = (\gamma + \varepsilon)/\delta, y∗=(α−ε)/βy^* = (\alpha - \varepsilon)/\beta: more prey on average, fewer predators. The prey's average is the number at which the predators just break even, so it rises when fishing makes the predators die faster.

Things to notice

  • In the equations a population can come very close to 0 but never reaches it. A real population that small can die out.
  • The coefficients and units are illustrative, not fitted to any real species.

Controls

  • The start cannot be on an axis: a population that starts at 0 stays at 0.

Related