188bet±¸ÓÃÍøÖ·

Applications of Systems of Differential Equations

Epidemiology: The Spread of Disease

Theoretical Introduction

Note: Before attempting this laboratory, you should first complete theÌýPredator-Prey Laboratory.

Let's consider modeling a new situation, namely, the way in which aÌýdiseaseÌýspreads through a population. The population could be made up of any life form, but let's keep things personal, and say that we're working with a human population. Also, let's talk about a rather familiar disease—influenza.

Choosing Variables

Before we can go very far in creating our model, we need to choose some variables to represent the quantities of interest. We'll need to talk about how many people have caught influenza already, how many people haven't caught it yet, and how many people have already recovered from the illness. Of course, we can think of each of these quantities as being a function of time. The following list of variable names should cover what we need:

• S(t)Ìýrepresents the number of people who have not yet caught the disease. The letterÌýSÌýis used because these people areÌýsusceptibleÌýto the illness.

• I(t)Ìýrepresents the number of people who currently have the disease. The letterÌýIÌýis used because these people areÌýinfectedÌýby the illness.

• R(t)Ìýrepresents the number of people who haveÌýrecoveredÌýfrom the illness. In the model we are creating, we consider patients that have already had the disease as having developed anÌýimmunityÌýto re-infection. This means that recovery from the illness doesÌýnotÌýreturn you to the susceptible group.Ìý(Note: the model would work in just the same way if the patients actually died from the illness. After all, once they die they are no longer susceptible or infected. Because of this fact, when some researchers work with this model they say that theÌýRÌýrepresents people that have beenÌýremovedÌýfrom the disease-interacting population, either by death or recovery. But let's be positive—they recovered!)

• tÌý°ù±ð±è°ù±ð²õ±ð²Ô³Ù²õÌýtime.

Finding the Rates of Change

Now let's think about how these sub-groups interact with one another. We'll consider the rates of change, (derivatives,) of each of our functions in turn.

³§â€²(t): the rate of change of the susceptible population. In the context of the spread of disease, this group is reduced by one person each time one of them becomes infected. So, the rate at whichÌýSÌýis changing is the same as the rate at which infection is taking place. The underlying question, then, is "When does infection take place?"

Infection occurs when the disease is passed from an infected person to a susceptible person. This occurs duringÌýsomeÌýof the times when susceptibles come in contact with infecteds. Now the number of susceptible-infected contacts will be proportional to the product ofÌýSÌý²¹²Ô»åÌýI, and of these contacts a proportion will catch the disease. It makes sense then to come up with the following formula for the rate of change of the susceptible population:

³§â€²Ìý=Ìý-Ìýa S I

·É³ó±ð°ù±ðÌýaÌýis some positive constant of proportionality which takes into account the total number of contacts between the two groups, and how many of these contacts actually lead to infection. The negative sign is present because when infection occurs, the susceptible population isÌýreduced.

±õ′(t): the rate of change of the infected population. We have already discussed why this population grows. They are constantly recruiting members from the susceptible group. We can even name the rate at which growth occurs due to this transformation—it'sÌýa S I. After all, the rate at which people are leaving the susceptible population is exactly the same as the rate at which they are joining the infected population.ÌýThat's where they are going!

But we have another factor to consider! There are also peopleÌýrecoveringÌýfrom the disease. It seems reasonable to assume that the number of people moving from the infected group to the recovered group is some kind of fixed proportion of the infected groups size, e.g. we can see the sense in a statement such as "On any given day, an average of 15% of our influenza patients enter recovery and are released from our care." So, the rate at which the infected group isÌýreducedÌýby recoveries should be a quantity that looks like -b I, ·É³ó±ð°ù±ðÌýbÌýis the constant of proportionality representing the recovery rate. (For the example above,ÌýbÌý=Ìý0.15.)

Putting our infection and recovery rates together, we get a quantity representing the rate of change of the infected population:

±õ′Ìý=Ìýa S IÌý-Ìýb I

Whether the infected group is growing or shrinking in size would be determined by the relative sizes of these two terms.

¸é′(t): the rate of change of the recovered population. We already know exactly what rate this group is growing at. Their members come exclusively from the recoveries we mentioned above, at a rate ofÌýb I. The formula for the rate of change of the recovered population is therefore:

¸é′Ìý=ÌýÌýb I

Note that this quantity is positive, since the recovered group is growing all of the time as people overcome the illness.

Putting it All TogetherÌý Ìý

So far, our three rates of change have led to three differential equations. If we throw in three initial conditions to keep them company we should have a nicely formed initial value problem.

The initial conditions are obvious. We'll use the sizes of the three groups at whatever time we're choosing to call time zero—the beginning of our analysis. Let's call these sizesÌýS_o,ÌýI_o, andÌýR_oÌýrespectively. The initial value problem becomes:

Make sure you understand why each part of the model makes sense. It's now time to begin our computer analysis, soÌýÌýand see how we might implement these ideas inÌýMathematica.

If you're lost, impatient, want an overview of this laboratory assignment, or maybe even all three, you can click on the compass button on the left to go to the table of contents for this laboratory assignment.