Chapter 4.2, Problem 106E

Differential equations can be used to model disease epidemics. In the next set of problems, we examine the change of size of two sub-populations of people living in a city: individuals who are infected and individuals who are susceptible to infection. S represents the size of the susceptible population, and I represents the size of the infected population. We assume that if a susceptible person interacts with an infected person, there is a probability c that the susceptible person will become infected. Each infected person recovers from the infection at a rate r and becomes susceptible again. We consider the case of influenza, where we assume that no one dies from the disease, so we assume that the total population size of the two sub-populations is a constant number, N. The differential equations that model these population sizes are

S' = rI − cSI and

I' =cSI − rI.

Here c represents the contact rate and r is the recovery rate.

106. [T] Use computational software or a calculator to compute the solution to the initial-value problem y'= ty, y(0) = 2 using Euler’s Method with the given step size h. Find the solution at t = 1. For a hint, here is “pseudo-code” for how to write a computer program to perform Euler’s Method for y' = f(t, y). y(0) = 2: Create function f(t, y) Define parameters y( 1) =y₀, t(0) = 0, step size N, and total number of steps, N Write a for loop: for k = 1 to N fn = f(t(k), y(k))

y(k+ 1) = y(k) + h*fn

t(k+1) = t(k)+ h