The Birth of Epidemic Modeling: Daniel Bernoulli and the Smallpox Inoculation Controversy

Author

Jong-Hoon Kim

Published

April 14, 2025

1 Introduction

Mathematical modeling of infectious diseases is now a cornerstone of public health decision-making. During COVID-19, we all became familiar with concepts like “flattening the curve” and reproduction numbers. But did you know that the first mathematical model for infectious disease control was developed nearly 300 years ago? In 1760, a Swiss mathematician named Daniel Bernoulli created a revolutionary approach to evaluate a controversial medical procedure: smallpox inoculation (). His work represents the birth of what we now call mathematical epidemiology ().

2 Who Was Daniel Bernoulli?

Daniel Bernoulli (1700-1782) came from a family of brilliant mathematicians. While his relatives focused on pure mathematics, Daniel pursued medicine alongside his mathematical studies (). This unique combination placed him perfectly to bridge these two worlds when a major public health controversy erupted in 18th century Europe.

Daniel Bernoulli

3 The Smallpox Crisis

Today, we celebrate smallpox as the first human disease to be completely eradicated through vaccination. But in Bernoulli’s time, smallpox was a devastating illness that:

  • Killed approximately 400,000 Europeans annually
  • Caused death in 1 out of every 7-10 infected individuals
  • Left survivors permanently scarred and sometimes blind
  • Affected nearly everyone, with most contracting it during childhood

4 The Controversial Solution: Variolation

In the early 18th century, a precursor to vaccination called “variolation” (or inoculation) was introduced to Europe (). This procedure involved:

  1. Taking a small amount of pus from a person with a mild case of smallpox
  2. Inserting it under the skin of a healthy person
  3. Inducing a usually milder form of the disease
  4. Conferring lifelong immunity if the person survived

This technique reduced the mortality rate from smallpox from about 15% to approximately 1-2%. However, it was extremely controversial because:

  • It deliberately infected healthy people with a deadly disease
  • Some people died from the procedure itself
  • Inoculated individuals could spread smallpox to others
  • Religious and ethical objections were raised against “interfering with God’s will”

5 Bernoulli’s Breakthrough: The First Disease Model

Bernoulli recognized that this controversy needed more than opinions—it needed data and mathematical analysis. In 1760, he presented his work “Essai d’une nouvelle analyse de la mortalité causée par la petite vérole” (Essay on a new analysis of the mortality caused by smallpox) to the Royal Academy of Sciences in Paris ().

His approach was revolutionary because:

  1. It was the first mathematical model specifically designed to evaluate a public health intervention
  2. It separated the risk of dying from smallpox from other causes of death
  3. It calculated potential years of life gained through inoculation
  4. It used real demographic data from the city of London

6 The Model Explained Simply

Bernoulli’s approach can be understood through a simple analogy:

Imagine 1,000 children born in the same year. Without inoculation, these children face two risks as they age: - The risk of contracting and possibly dying from smallpox - The risk of dying from any other cause

Bernoulli used available data to estimate: - How many people would get smallpox each year - How many would die from it - How many would die from other causes

He then compared this natural scenario with one where everyone was inoculated against smallpox. Although inoculation carried its own small risk of death, Bernoulli’s calculations showed that widespread inoculation would:

  • Increase life expectancy by about 3 years (significant for that time period)
  • Save thousands of lives across a population

7 Mathematical Innovation: The Life Table Approach

Bernoulli’s work was mathematically innovative because he:

  1. Created what we now call a “dynamic life table” - tracking population changes over time
  2. Separated specific causes of death (what epidemiologists now call “competing risks”)
  3. Introduced differential equations to model population changes due to disease ()

In modern mathematical notation, his core equation looked something like:

dp(a)da=μ(a)p(a)λps(a)

Where: - p(a) represents the proportion of people surviving to age a - μ(a) is the force of mortality from causes other than smallpox - λ is the rate of smallpox infection - ps(a) is the proportion of susceptible people at age a

While this may look intimidating, the concept is straightforward: Bernoulli tracked how many people would die from different causes under different scenarios ().

8 Impact and Legacy

Bernoulli’s mathematical analysis had profound effects:

  1. Quantified benefits: He showed that inoculation would increase life expectancy by approximately 3 years - remarkable in an era when life expectancy was around 30 years ().

  2. Risk comparison: He demonstrated that the 1-2% risk of death from inoculation was far smaller than the lifetime risk of dying from natural smallpox infection.

  3. Policy influence: His work helped shift public opinion and policy toward accepting inoculation as a worthwhile public health measure.

  4. Methodological innovation: He established mathematical modeling as a valid approach to public health questions ().

9 The Birth of Mathematical Epidemiology

Bernoulli’s approach contained the seeds of modern infectious disease modeling:

  • It used mathematics to predict disease outcomes
  • It quantified the impact of interventions
  • It separated competing causes of death
  • It estimated population-level effects from individual risks

The field that Bernoulli pioneered would later be expanded by other mathematical innovators like Ronald Ross (, ) and the team of Kermack and McKendrick (), eventually leading to the sophisticated models we use today ().

10 Modern Applications of Bernoulli’s Approach

The principles Bernoulli established are still used to answer questions like:

  • How many lives will a new vaccine save? ()
  • Is a screening program worth its costs?
  • Which age groups should receive priority for limited medical resources?
  • What is the optimal strategy for controlling an epidemic? ()

Modern epidemic models have incorporated additional mathematical sophistication, particularly around the concept of the basic reproduction number (R₀) and the construction of next-generation matrices (; ), but they still follow Bernoulli’s fundamental approach of using mathematics to guide public health decisions.

11 Conclusion

Daniel Bernoulli’s 1760 smallpox model represents a pivotal moment in the history of public health (). By applying mathematical reasoning to a contentious medical issue, he demonstrated that data and analysis could guide decisions more effectively than opinions alone.

As we continue to face infectious disease challenges in the modern world, from emerging pathogens to vaccine hesitancy, Bernoulli’s approach remains relevant (). His work reminds us that mathematical modeling is not just an academic exercise but a powerful tool for saving lives and informing policy.

The next time you hear about disease models predicting the course of an epidemic or evaluating the impact of interventions, remember that this approach began with a mathematician applying his skills to a pressing public health problem nearly three centuries ago ().

12 References

Anderson, Roy M., and Robert M. May. 1991. “Infectious Diseases of Humans: Dynamics and Control.” Oxford University Press.
Bernoulli, Daniel. 1766. “Essai d’une Nouvelle Analyse de La Mortalité Causée Par La Petite vérole.” Mémoires de Mathématiques Et de Physique, Académie Royale Des Sciences, 1–45.
Blower, Sally, and Daniel Bernoulli. 2004. “An Attempt at a New Analysis of the Mortality Caused by Smallpox and of the Advantages of Inoculation to Prevent It.” Reviews in Medical Virology 14 (5): 275–88.
Diekmann, Odo, J. A. P. Heesterbeek, and Michael G. Roberts. 2010. “The Construction of Next-Generation Matrices for Compartmental Epidemic Models.” Journal of the Royal Society Interface 7 (47): 873–85.
Dietz, Klaus, and J. A. P. Heesterbeek. 2002. “Daniel Bernoulli’s Epidemiological Model Revisited.” Mathematical Biosciences 180 (1-2): 1–21.
Ferguson, Neil M., Derek A. T. Cummings, Christophe Fraser, James C. Cajka, Philip C. Cooley, and Donald S. Burke. 2006. “Strategies for Mitigating an Influenza Pandemic.” Nature 442 (7101): 448–52.
Hald, Anders. 2003. History of Probability and Statistics and Their Applications Before 1750. Hoboken, NJ: John Wiley & Sons.
Heesterbeek, Hans, and Mick G. Roberts. 2015. “How Mathematical Epidemiology Became a Field of Biology: A Commentary on Anderson and May (1981) ’the Population Dynamics of Microparasites and Their Invertebrate Hosts’.” Philosophical Transactions of the Royal Society B 370 (1666): 20140307.
Heesterbeek, J. A. P., and Klaus Dietz. 1996. “Mathematical Epidemiology of Infectious Diseases: Model Building, Analysis and Interpretation.” Mathematical Biosciences 146 (1): 77–79.
Hethcote, Herbert W. 2000. “The Mathematics of Infectious Diseases.” SIAM Review 42 (4): 599–653.
Keeling, Matt J., and Pejman Rohani. 2008. “Modeling Infectious Diseases in Humans and Animals.” Princeton University Press.
Kermack, William O., and Anderson G. McKendrick. 1927. “A Contribution to the Mathematical Theory of Epidemics.” Proceedings of the Royal Society of London. Series A 115 (772): 700–721.
Riley, Steven, Christophe Fraser, Christl A. Donnelly, Azra C. Ghani, Laith J. Abu-Raddad, Anthony J. Hedley, Gabriel M. Leung, et al. 2003. “Transmission Dynamics of the Etiological Agent of SARS in Hong Kong: Impact of Public Health Interventions.” Science 300 (5627): 1961–66.
Ross, Ronald. 1908. “Report on the Prevention of Malaria in Mauritius.”
———. 1911. The Prevention of Malaria. London: John Murray.
Wallinga, Jacco, and Marc Lipsitch. 2007. “How Generation Intervals Shape the Relationship Between Growth Rates and Reproductive Numbers.” Proceedings of the Royal Society B: Biological Sciences 274 (1609): 599–604.