How to compute the Lotka-Volterra equation: system of coupled non-linear ordinary differential equations describing prey-predator population dynamics in closed ecosystems. Formulas: dx/dt = αx - βxy (prey) and dy/dt = δxy - γy (predators), where x(t) prey count, y(t) predator count, α prey birth rate without predators, β predation rate hunting efficiency encounters, δ predator conversion rate prey-to-new-predators, γ predator natural mortality without prey. Equilibrium: x* = γ/δ, y* = α/β (non-linear stable center). Oscillation period approximate T = 2π/√(αγ). Conserved quantity V = δx - γ ln(x) + βy - α ln(y) first integral along closed trajectories. Non-linear systems produce closed oscillatory cycles around equilibrium, characteristic delay response predator-prey dynamics. History: Alfred James Lotka (USA 1880-1949) Johns Hopkins University Baltimore MD biophysicist 1925 paper "Elements of Physical Biology" Williams Wilkins Baltimore introducing autocatalytic equations reaction kinetics applied populations;

Vito Volterra (Italy 1860-1940) Università La Sapienza Roma mathematician professor mathematics 1881-1932 founder CNR Consiglio Nazionale delle Ricerche Roma 1923, son-in-law Umberto D Ancona Adriatico zoologist 1925 statistical observations Adriatic fisheries Trieste Venezia Pescara Bari showing during WWI 1914-1918 when fisheries reduced predator fish fraction (sharks tuna) increased relative to prey (sardines), 1926 paper "Variazioni e fluttuazioni del numero d individui in specie animali conviventi" Memorie Reale Accademia Nazionale Lincei Roma 2:31-113 derived independent same Lotka equations. Essential tool for: population ecology Italian Adriatic sardine-tuna Trieste Pescara Ancona, lynx-hare Canada Hudson Bay Company records 1845-1935 famous classic dataset, wolf-deer Yellowstone USA reintroduction 1995 Wolf Restoration Project National Park Service NPS, ISPRA Istituto Superiore Protezione Ricerca Ambientale Roma; fisheries management NMFS National Marine Fisheries Service NOAA Silver Spring MD, NOAA Fisheries Magnuson-Stevens Act 1976, FAO Food Agriculture Organization Roma 1945; ecology research USA NSF Long-Term Ecological Research LTER 26 sites founded 1980; Yellowstone National Park 1872 first national park world; USA universities ecology: UC Davis CA Department Wildlife Fish Conservation Biology; UC Berkeley CA Environmental Science Policy Management ESPM; Cornell Ithaca NY Department Natural Resources;

Yale New Haven CT School Forestry Environmental Studies; Stanford CA Hopkins Marine Station Biology; University Wisconsin-Madison Aldo Leopold Center Wildlife Ecology; UC Santa Cruz CA Department Ecology Evolution; University Washington WA Center Quantitative Science. NMFS NOAA Silver Spring MD fisheries management Magnuson-Stevens Act regulation US fisheries Pacific Atlantic Caribbean Gulf, MMPA Marine Mammal Protection Act 1972, ESA Endangered Species Act 1973.

Lotka-Volterra formula

Differential equations system: dx/dt = αx - βxy (prey population), dy/dt = δxy - γy (predator population). 2D non-linear ODE first order system. α positive: prey grow exponentially without predators. β positive: predation reduces prey proportional to product x*y interaction term. δ positive: predators reproduce by capturing prey, conversion factor. γ positive: predators die naturally without prey.

Non-trivial equilibrium: x* = γ/δ, y* = α/β. Internal critical point non-linear stable but NEUTRAL center (Hopf-stable for ODE Hamiltonian systems): small perturbations around equilibrium produce closed periodic cycles NOT convergence. Trivial equilibrium (0,0) unstable.

Hamiltonian conserved quantity: V(x,y) = δx - γ ln(x) + βy - α ln(y) constant along closed trajectories. Hamiltonian system canonical form, energy conservation analogous. Level curves V=constant in (x,y) plane represent closed cyclic trajectories, no point attractor convergence.

Linear oscillation period: T = 2π/√(αγ) approximation for small oscillations around equilibrium (Jacobian linearization). Period depends only on αγ NOT β δ. For α=1, γ=1.5: T = 2π/√1.5 = 5.13 time units. Non-linear cycles period larger than linear for large amplitudes (anharmonicity).

Phase shift predators vs prey: predator cycles lag π/2 (quarter period) compared to prey cycles. Prey peak precedes predator peak, prey decline precedes predator decline, fundamental time lag. Maximum prey when dy/dt changes sign (predators start growing), maximum predators when dx/dt changes sign (prey start declining).

History and context

Vito Volterra (Italy 1860-1940) Ancona-Roma mathematician born Anchora Marche Jewish origin cultured family father merchant, professor mathematics University Pisa 1883 age 23, University Torino 1893, University La Sapienza Roma 1900-1932 (until resignation refusing fascism oath Mussolini 1931). Senator Kingdom Italy 1905 appointment King Vittorio Emanuele III. Founder CNR Consiglio Nazionale delle Ricerche Roma 1923 first president, Italian National Research Council. Contributed functional analysis theories (1887 paper "Sur les fonctions qui dépendent d'autres fonctions"), integral equations (1896 Volterra equations), Volterra series non-linear dynamic systems expansions, mathematical biology populations 1925-1930 after son-in-law Umberto D Ancona zoologist Adriatic statistical analysis fisheries data WWI 1914-1918. 1925-1926 paper Lincei Roma "Variazioni e fluttuazioni del numero d'individui in specie animali conviventi" Memorie Reale Accademia Nazionale Lincei 2:31-113 derived prey-predator equations same Alfred Lotka independent (Lotka 1925 'Elements of Physical Biology' Williams Wilkins Baltimore first chronological publication). Heritage: Volterra without life post-fascism Italian university, Jewish anti-fascist Mussolini 1938 racial laws persecution died Roma 1940 age 80; CNR successor Guglielmo Marconi 1928, Italy post-war scientific research re-foundation 1945.

Alfred James Lotka (USA 1880-1949) Lemberg Galicia Austria-Hungary born 1880 emigrated USA 1909 Cornell PhD biology 1909, Johns Hopkins University Baltimore Maryland 1922-1949 statistician actuary insurance companies Metropolitan Life Insurance Co NYC 1924-1947 chief, biophysicist mathematical biologist 1925 'Elements of Physical Biology' Williams Wilkins Baltimore first book mathematical biology. Studied autocatalytic reactions chemical kinetics applied to population dynamics produce same prey-predator equations. Lotka curve population growth, demographic transition Lotka theorem stable populations. Umberto D Ancona (Italy 1896-1964) Vienna-Roma zoologist marine biologist Trieste son-in-law Volterra 1925 doctoral thesis University Trieste statistical analysis Adriatic fisheries pre-WWI vs WWI 1914-1918 vs post-WWI, observed during war fisheries reduction predator fraction (sharks tuna) increased 12-15% compared to prey (sardines anchovies), pre-war 5-7% predators, post-WWI return 5-7%. Hypothesis: reduced fishing during war causes natural prey-predator dynamics undisturbed predator population recovery. USA datasets famous predator-prey: Hudson Bay Company Canada records 1845-1935 90 years canadian lynx-snowshoe hare cycles 9.6-year period oscillations Mackenzie River Yukon Northwest Territories fur trade ledgers Charles Elton Oxford 1924 statistical analysis classical reference;

Yellowstone Wolf Restoration Project National Park Service NPS 1995 reintroduction grey wolves Canis lupus 31 wolves McMurchie Mackenzie Pack Druid Pack from Canadian Rockies, Yellowstone elk Cervus elaphus population reduced from 19000 1995 to 5000 2024 trophic cascade beaver willow regeneration. USA: NSF NOAA NMFS National Marine Fisheries Service Silver Spring MD fisheries management Magnuson-Stevens Act 1976; ESA Endangered Species Act 1973; MMPA Marine Mammal Protection Act 1972; CITES Convention International Trade Endangered Species 1973; USFS US Forest Service Washington DC; USFWS US Fish Wildlife Service Falls Church VA. USA Yellowstone National Park established 1872 first national park world, Wyoming Montana Idaho 8983 km² area; National Park Service NPS Washington DC founded 1916 426 parks managed. USA universities ecology: UC Davis CA Department Wildlife Fish Conservation Biology; UC Berkeley CA Environmental Science Policy Management ESPM; Cornell Ithaca NY Department Natural Resources; Yale New Haven CT School Forestry Environmental Studies; Stanford CA Hopkins Marine Station Biology; University Wisconsin-Madison Aldo Leopold Center Wildlife Ecology; UC Santa Cruz CA Department Ecology Evolution; University Washington WA Center Quantitative Science. Famous mathematicians USA biology: Robert May (Australia-USA-UK 1936-2020) Princeton-Oxford ecology mathematics theoretical Lord May Oxford; James D.

Murray (UK 1931-) Oxford 'Mathematical Biology' textbook 2002 standard 3rd edition Springer; Lee Segel (USA 1932-2005) UC Berkeley NY mathematical biology; Joel E. Cohen (USA 1944-) Rockefeller University NY ecology; Daniel Botkin (USA) Yale ecology theoretical. Italian famous mathematicians: Volterra Vito 1860-1940; Tullio Levi-Civita 1873-1941 Padua tensor calculus general relativity Einstein contributions; Federigo Enriques 1871-1946 Bologna projective geometry algebraic geometry Italian school; Beniamino Segre 1903-1977 Bologna algebraic geometry; Ennio De Giorgi 1928-1996 Pisa SNS calculus variations Hilbert problem 19 1957 solution. Software simulation: SciPy Python scipy.integrate.odeint differential equations, R deSolve package, MATLAB ode45 ode23 ode15s, Mathematica NDSolve, Berkeley Madonna; XPPaut bifurcation analysis Bard Ermentrout Pittsburgh PA; AUTO continuation software. Lotka-Volterra extensions: competitive Lotka-Volterra two-species competition same resources; logistic Lotka-Volterra carrying capacity correction realistic; n-species generalization community ecology Hubbell unified neutral theory 2001; SIR epidemiology Kermack-McKendrick 1927 disease modeling.

Lotka-Volterra parameters examples table

• Lotka-Volterra classic textbook: alpha=1, beta=0.1, delta=0.075, gamma=1.5, x*=20, y*=10
• Lynx-hare Canada Hudson Bay: alpha=0.5, beta=0.025, delta=0.005, gamma=0.5, T=12.6 years
• Wolf-elk Yellowstone post-1995: alpha=0.3, beta=0.001, delta=0.0002, gamma=0.2, x*=1000 elk, y*=300 wolves
• Sardine-tuna Adriatic Volterra 1925: alpha=2, beta=0.05, delta=0.02, gamma=1.5, x*=75, y*=40
• Fox-rabbit Britain: alpha=0.6, beta=0.04, delta=0.01, gamma=0.4, x*=40, y*=15
• Equilibrium prey x*: gamma/delta
• Equilibrium predator y*: alpha/beta
• Linear period T: 2*pi/sqrt(alpha*gamma)
• Phase shift: π/2 quarter period predators lag
• Conserved quantity V: delta*x - gamma*ln(x) + beta*y - alpha*ln(y)
• Yellowstone elk 1995: 19000 individuals pre-reintroduction wolves
• Yellowstone elk 2024: 5000 individuals post-reintroduction
• Yellowstone wolves 2024: 100 individuals 9 packs
• Hudson Bay records 1845-1935: 90 years canadian lynx-snowshoe hare
• Lynx-hare period: 9.6 years average cycle
• Trophic cascade Yellowstone: wolves -> elk -> willow -> beaver positive
• Adriatic WWI 1914-1918: predators 12-15% post-war return 5-7%
• Mediterranean sardine: Sardina pilchardus stock
• Mediterranean anchovies: Engraulis encrasicolus stock
• Mediterranean bluefin tuna: Thunnus thynnus largest tuna species
• Mazara del Vallo: largest Mediterranean fish market Sicily
• Pescara Adriatic: Italian fisheries port
• Trapani Sicily: sardine industry historical
• Italian fish markets: Mazara, Trapani, Pescara, Anzio, Catania
• Volterra CNR foundation: 1923 first president Roma
• Lotka 1925 publication: Elements of Physical Biology Williams Wilkins Baltimore
• Volterra 1926 publication: Lincei Roma 2:31-113
• D Ancona thesis 1925: University Trieste Adriatic fisheries WWI analysis
• FAO foundation: 1945 Roma viale delle Terme di Caracalla

Detailed numerical examples

Example 1 - Lotka-Volterra classic textbook. Parameters: alpha=1.0 (prey grow fast), beta=0.1 (moderate predation), delta=0.075 (moderate conversion efficiency), gamma=1.5 (moderate predator mortality). Equilibrium: x* = 1.5/0.075 = 20 prey, y* = 1.0/0.1 = 10 predators. Linear period T = 2*pi/sqrt(1*1.5) = 5.13 time units. Initial x0=10, y0=5: point below prey equilibrium above predator equilibrium. Clockwise cycles in (x,y) plane, closed trajectories, no convergence.

Example 2 - Lynx-hare Canada Hudson Bay 9.6-year period. Hudson Bay Company Canada records 1845-1935 90 years canadian lynx-snowshoe hare cycles classical dataset Mackenzie River Yukon Northwest Territories fur trade ledgers Charles Elton Oxford 1924 statistical analysis. Estimated parameters: alpha=0.5, beta=0.025, delta=0.005, gamma=0.5. Equilibrium x* = 0.5/0.005 = 100 hares, y* = 0.5/0.025 = 20 lynxes. Period T = 2*pi/sqrt(0.5*0.5) = 12.6 years vs observed 9.6 years (Lotka-Volterra approximation).

Example 3 - Wolf-elk Yellowstone post-1995. Yellowstone Wolf Restoration Project National Park Service NPS 1995 reintroduction 31 grey wolves Canis lupus from Canadian Rockies. Parameters: alpha=0.3 (elk reproduction), beta=0.001 (predation per wolf-elk encounter), delta=0.0002 (conversion), gamma=0.2 (wolves natural mortality). Equilibrium x* = 0.2/0.0002 = 1000 elk, y* = 0.3/0.001 = 300 wolves. Trophic cascade observed wolves -> elk reduction 19000 1995 to 5000 2024 -> willow recovery -> beaver populations increase positive ecological cascade USA NPS.

Example 4 - Sardine-tuna Adriatic Volterra 1925. Italy observations Adriatic Trieste Pescara Ancona Bari WWI 1914-1918 fisheries reduction. Historical parameters Volterra 1925: alpha=2 (sardine reproduction high), beta=0.05 (moderate predation), delta=0.02 (tuna conversion efficiency moderate), gamma=1.5 (tuna mortality). Equilibrium x* = 1.5/0.02 = 75 (sardine units), y* = 2/0.05 = 40 (tuna units). WWI fishery reduction 1914-1918 increased predator fraction 12-15% from pre-war 5-7%, evidence Lotka-Volterra dynamics naturally restoring predator population.

Example 5 - Fox-rabbit Britain. UK British Isles wildlife Vulpes vulpes red fox Oryctolagus cuniculus rabbit. Estimated parameters: alpha=0.6 (rabbits high reproduction), beta=0.04, delta=0.01, gamma=0.4. Equilibrium x* = 0.4/0.01 = 40 rabbits/km², y* = 0.6/0.04 = 15 foxes/km². Period T = 2*pi/sqrt(0.6*0.4) = 12.8 years. Observed UK studies typical 9-12 year cycles. USA equivalent: coyote-rabbit, owls-mice, white-tailed deer-cougars Western states.

Limits and considerations

1. Idealized model: assumes constant environment, no resource competition, no spatial heterogeneity, no carrying capacity (logistic correction needed realistic).
2. Constant parameters: real ecology variables seasonal, climate change, human impact, harvest pressure modify dynamics.
3. Only 2 species: real ecosystems multi-species interactions food web complex, n-species generalization required community ecology.
4. Neutral cycles: NOT observed real systems, perturbations cause damped or amplified oscillations, logistic Lotka-Volterra preferred.
5. Type I functional response (linear xy): Holling Type II saturation real predators, Type III sigmoidal switching.
6. Continuous deterministic: real populations stochastic discrete individuals, demographic stochasticity important small populations extinction risk.
7. Adriatic Volterra 1925 data: limited fisheries records WWI, modern statistical methods may interpret differently.
8. Yellowstone simplification: trophic cascade complex, multiple predator species (bears coyotes), elk migration patterns.
9. Numerical integration: Euler method simple but inaccurate, Runge-Kutta RK4 preferred long-term simulations.
10. V conservation quantity: theoretical exact, numerical drift accumulates float64 long simulations.

Common errors

1. Confusing stable with unstable equilibrium: Lotka-Volterra NEUTRAL center not stable attractor.
2. Neglecting predator-prey phase shift π/2: critical for understanding prediction.
3. Pretending cycle amplitude independent of initial conditions: depends on V conserved quantity, larger amplitude larger V.
4. Confusing Lotka 1925 with Volterra 1926: Lotka first chronological publication, Volterra independent derivation.
5. Ignoring logistic correction realistic: Lotka-Volterra ideal, real systems have carrying capacity.
6. Pretending linear T = 2π/√(αγ) exact: only for small oscillations linearization.
7. JavaScript Math.log natural ln (NOT log10): check conserved quantity formula.
8. Numerical Euler integration unstable: small dt needed convergence.
9. Confusing x prey y predator indices: convention varies, check derivation carefully.
10. Neglecting extinction risk: stochastic populations small risk extinction not in deterministic ODE.

Numbered applications

1. USA Yellowstone Wolf Restoration: National Park Service NPS 1995 reintroduction 31 grey wolves Canis lupus, trophic cascade observation 1995-2024, elk 19000->5000.
2. Hudson Bay Company Canada: records 1845-1935 90 years canadian lynx-snowshoe hare cycles classical dataset Charles Elton Oxford analysis 1924.
3. USA universities ecology: UC Davis Wildlife Conservation; UC Berkeley ESPM; Cornell Ithaca Natural Resources; Yale Forestry; Stanford Hopkins Marine; Wisconsin-Madison Aldo Leopold Center.
4. NMFS NOAA fisheries USA: National Marine Fisheries Service Silver Spring MD Magnuson-Stevens Act 1976 regulation US fisheries Pacific Atlantic Caribbean Gulf.
5. Conservation acts USA: ESA Endangered Species Act 1973; MMPA Marine Mammal Protection Act 1972; CITES Convention International Trade Endangered Species 1973.
6. Italy Volterra Adriatic legacy: 1925 observations sardine-tuna Trieste Pescara Ancona; D Ancona Università Trieste 1925 thesis WWI fisheries.
7. FAO Roma: Food Agriculture Organization founding 1945 viale delle Terme di Caracalla, fisheries management worldwide.
8. Software simulation: SciPy Python scipy.integrate.odeint differential equations; R deSolve package; MATLAB ode45 ode23 ode15s; XPPaut bifurcation analysis Bard Ermentrout Pittsburgh PA.
9. Famous USA mathematical biologists: Robert May Princeton-Oxford ecology Lord May; James D Murray Oxford 'Mathematical Biology' textbook standard; Joel E Cohen Rockefeller University NY ecology.
10. Lotka-Volterra extensions: competitive same resources two-species; logistic carrying capacity; n-species community ecology Hubbell unified neutral theory 2001; SIR epidemiology Kermack-McKendrick 1927.

For support

For technical questions: info@calcoi.com. Educational mathematical ecology tool. Applications Lotka A 1925 'Elements of Physical Biology' Williams Wilkins Baltimore, Volterra V 1926 'Variazioni e fluttuazioni del numero d individui in specie animali conviventi' Memorie Reale Accademia Nazionale Lincei Roma 2:31-113, D Ancona U 1925 University Trieste thesis Adriatic fisheries WWI, Murray JD 2002 'Mathematical Biology' 3rd edition Springer.

Updated July 2026