Mathematical derivation of the Lotka–Volterra equations
- Assumptions (minimal)
- Two species: prey x(t) and predator y(t).
- Prey reproduce exponentially with per-capita rate α when predators absent.
- Predation rate proportional to encounters ~ x·y (law of mass action). Each encounter removes prey at rate β and yields predator growth proportional to prey consumed with factor δ.
- Predator dies at per-capita rate γ in absence of prey.
- Well-mixed population, continuous time, parameters α, β, γ, δ > 0.
- Construct per-capita rates → differential equations
- Prey growth without predators: dx/dt = α x.
- Loss by predation: encounters ~ xy, so subtract β x y.
- Predator death without prey: −γ y.
- Predator gain from feeding: +δ x y.
Combine to get the canonical system dx/dt = α x − β x y dy/dt = −γ y + δ x y
- Fixed points
- Solve dx/dt = dy/dt = 0 → (x,y) = (0,0) and (x,y) = (γ/δ, α/β).
- Linear stability (Jacobian) J(x,y) = [[α − β y, −β x], [δ y, δ x − γ]].
- At (0,0): eigenvalues λ1 = α, λ2 = −γ → saddle (unstable).
- At (γ/δ, α/β): Jacobian has purely imaginary eigenvalues ±i√(αγ) → center (neutral stability, closed orbits for the ideal model).
-
First integral (conserved quantity) Divide equations: (dx/dt)/(x) − (β/δ)(dy/dt)/(y) = 0 and integrate to obtain constant of motion V(x,y) = δ x − γ ln x + β y − α ln y = const. Level curves V(x,y)=C are closed orbits → neutral periodic oscillations (no damping).
-
Remarks (concise)
- Model is minimal and idealized; adding density-dependence (e.g., logistic prey term), functional responses, or stochasticity changes stability and realism.
- The system is Hamiltonian-like with the conserved first integral; orbits are neutrally periodic (not asymptot
Leave a Reply