Problem 38

Question

The spread of a non-fatal disease through a population of fixed size \(M\) can be modeled as follows. The rate that healthy people are infected, in people per day, is proportional to the product of the numbers of healthy and infected people. The constant of proportionality is \(0.01 / M .\) The rate of recovery, in people per day, is 0.009 times the number of people infected. Construct a differential equation that models the spread of the disease. Assuming that initially only a small number of people are infected, plot a graph of the number of infected people against time. What fraction of the population is infected in the long run?

Step-by-Step Solution

Verified
Answer
In the long run, 10% of the population is infected.
1Step 1: Define Variables
Let's define our variables. Let \( I(t) \) be the number of infected people at time \( t \). Let \( S(t) \) be the number of healthy people. Since we have a fixed population size \( M \), we have that \( S(t) + I(t) = M \).
2Step 2: Determine Infection and Recovery Rates
The infection rate is proportional to the product of the number of healthy and infected people. Therefore, the infection rate can be written as \( \frac{0.01}{M} \cdot S(t) \cdot I(t) = 0.01 I(t) \left( M - I(t) \right) / M \). The recovery rate is \( 0.009 I(t) \).
3Step 3: Formulate the Differential Equation
The rate of change of infected people \( \frac{dI}{dt} \) is the infection rate minus the recovery rate. Hence, our differential equation is: \[ \frac{dI}{dt} = 0.01 I(t) \frac{(M - I(t))}{M} - 0.009 I(t). \]
4Step 4: Solve the Differential Equation
Due to the complexity, solving the differential equation analytically can be challenging without simplification assumptions or numerical methods. We can explore its behavior by examining equilibrium points or using computational software for numerical solutions.
5Step 5: Analyze Long-Term Behavior
At equilibrium, the number of infected people does not change, meaning \( \frac{dI}{dt} = 0 \). Setting our equation to 0, we solve for \( I(t) \) to find potential equilibrium fractions: \[ 0.01 I(t) \frac{(M - I(t))}{M} = 0.009 I(t) \] leads to the solutions **I(t) = 0** or **I(t) = 0.1 M**. Eventually, a fraction of 0.1 of the population is infected in the long run.

Key Concepts

Disease ModelingPopulation DynamicsEquilibrium Analysis
Disease Modeling
Disease modeling is a crucial method for understanding how diseases spread through populations. It helps predict the future course of an outbreak and assesses the potential impact of interventions. In the provided problem, we use a differential equation to model the spread of a non-fatal disease, which assumes a constant population size.

In this scenario, the infection rate is proportional to the number of healthy and infected individuals. This is a common assumption in disease modeling, recognizing that more interactions between healthy and infected individuals increase the chances of transmission. The constant of proportionality is given as \( \frac{0.01}{M} \), which means the rate changes directly with population size.

Moreover, the model includes a recovery rate, expressed as a constant rate of 0.009 times the number of infected individuals. This helps us understand how quickly people can recover from the disease. By putting these elements together, we form a differential equation that provides insights into the spread of the disease over time. Such models are fundamental in epidemiology, helping to design control strategies through understanding the dynamics of disease transmission.
Population Dynamics
Population dynamics involve the study of short and long-term changes in the size and composition of populations, and the processes driving these changes. In our model, the population is composed of healthy and infected individuals, with the total population size \( M \) assumed constant.

The dynamics in this model are driven by:
  • Infection rates: Related to how healthy individuals come into contact with infected individuals, leading to new infections.
  • Recovery rates: Representing how infected individuals return to a healthy state.
In our problem, the interaction between these rates alters the number of infected individuals over time, governing the population dynamics of the disease. The differential equation \( \frac{dI}{dt} = 0.01 I(t) \frac{(M - I(t))}{M} - 0.009 I(t) \) captures these dynamics precisely.

Understanding these dynamics is key in predicting how quickly a disease will spread or die out, and the peak number of individuals that might be infected at any given time. It also informs decisions on health strategies and allocation of limited resources, such as vaccines or other medical interventions.
Equilibrium Analysis
Equilibrium analysis in differential equations examines points at which the system does not change over time, meaning that the rate of change, \( \frac{dI}{dt} \), is zero. At these points, the population dynamics have reached a stable state.

For the given differential equation, setting \( \frac{dI}{dt} = 0 \) allows us to find the equilibrium states. Solving our specific equation, we find two solutions: \( I(t) = 0 \) and \( I(t) = 0.1 M \). At \( I(t) = 0 \), the disease has been eradicated and at \( I(t) = 0.1 M \), 10% of the population remains infected in the long run.

Such analysis is essential because it informs us about the potential outcomes of the disease spreading through the population. Knowing the stable states, health authorities can prepare and respond better to an outbreak. They also indicate the effectiveness of public health interventions in possibly shifting these equilibrium points towards more favorable outcomes.