(*^ ::[ Information = "This is a Mathematica Notebook file. It contains ASCII text, and can be transferred by email, ftp, or other text-file transfer utility. It should be read or edited using a copy of Mathematica or MathReader. If you received this as email, use your mail application or copy/paste to save everything from the line containing (*^ down to the line containing ^*) into a plain text file. On some systems you may have to give the file a name ending with ".ma" to allow Mathematica to recognize it as a Notebook. The line below identifies what version of Mathematica created this file, but it can be opened using any other version as well."; FrontEndVersion = "Macintosh Mathematica Notebook Front End Version 2.2"; MacintoshStandardFontEncoding; fontset = title, inactive, noPageBreakBelow, nohscroll, preserveAspect, groupLikeTitle, center, M7, e8, 24, "B Univers 65 Bold"; fontset = subtitle, inactive, noPageBreakBelow, nohscroll, preserveAspect, groupLikeTitle, center, M7, e6, 18, "B Garamond Bold"; fontset = subsubtitle, inactive, noPageBreakBelow, nohscroll, preserveAspect, groupLikeTitle, center, M7, e6, 14, "I Garamond LightItalic"; fontset = section, inactive, noPageBreakBelow, nohscroll, preserveAspect, groupLikeSection, grayBox, M22, a20, 18, "B Univers 65 Bold"; fontset = subsection, inactive, noPageBreakBelow, nohscroll, preserveAspect, groupLikeSection, blackBox, M19, a15, 14, "AGaramond Semibold"; fontset = subsubsection, inactive, noPageBreakBelow, nohscroll, preserveAspect, groupLikeSection, whiteBox, M18, bold, a12, 12, "Garamond"; fontset = text, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, 12, "AGaramond"; fontset = smalltext, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, 10, "AGaramond"; fontset = input, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeInput, M42, N23, bold, L-5, 12, "Courier"; fontset = output, output, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeOutput, M42, N23, L-5, 12, "Courier"; fontset = message, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeOutput, M42, N23, R65535, L-5, 12, "Courier"; fontset = print, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeOutput, M42, N23, L-5, 12, "Courier"; fontset = info, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeOutput, M42, N23, B65535, L-5, 12, "Courier"; fontset = postscript, PostScript, formatAsPostScript, output, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeGraphics, M7, l34, w282, h287, 12, "Courier"; fontset = name, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, italic, 10, "Geneva"; fontset = header, inactive, noKeepOnOnePage, preserveAspect, M7, 12, "Garamond"; fontset = leftheader, inactive, L2, 12, "Garamond"; fontset = footer, inactive, noKeepOnOnePage, preserveAspect, center, M7, 12, "Garamond"; fontset = leftfooter, inactive, L2, 12, "Garamond"; fontset = help, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, 10, "Times"; fontset = clipboard, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, 12, "Times"; fontset = completions, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, 12, "Times"; fontset = special1, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M18, N55, 12, "Garamond"; fontset = special2, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, 12, "Times"; fontset = special3, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, 12, "Times"; fontset = special4, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, 12, "Times"; fontset = special5, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, 12, "Times"; showRuler; currentKernel; ] :[font = text; inactive; preserveAspect; leftWrapOffset = 17; leftNameWrapOffset = 1] Mathematics 162 Laboratory 8 Week of March 22, 1993 Name: _____________________________ Lab Partner: ___________________________ Consulted with: ____________________________________________________________ :[font = smalltext; inactive; preserveAspect; right] © Lafayette College, 1994. :[font = title; inactive; preserveAspect] A Model for an Epidemic :[font = text; inactive; preserveAspect] Suppose a contagious disease has appeared in a population. The disease is not fatal, but people who have recovered are still susceptible to the disease. We will use sequences to build a mathematical model for the spread of the disease in the population. :[font = section; inactive; preserveAspect] Part 0: The Model :[font = text; inactive; preserveAspect] Define the following sequences :[font = text; inactive; preserveAspect] In = the proportion of the population that has the disease on day n; Fn = the proportion of the population that is newly infected on day n; Cn = the proportion of the population that recovers on day n. ;[s] 15:0,1;1,0;2,2;66,1;67,2;69,1;70,0;71,2;137,1;138,2;140,1;141,0;142,2;199,1;200,2;202,-1; 3:3,17,10,AGaramond,66,9,0,0,0;6,14,10,AGaramond,2,12,0,0,0;6,14,10,AGaramond,0,12,0,0,0; :[font = text; inactive; preserveAspect] From the definitions, it follows that :[font = text; inactive; preserveAspect] In = InÊÐ 1 + Fn Ð Cn. ;[s] 12:0,1;31,2;32,3;35,1;36,2;37,0;41,1;45,2;46,3;51,1;52,2;53,3;55,-1; 4:1,17,10,AGaramond,64,9,0,0,0;4,14,10,AGaramond,2,12,0,0,0;4,17,10,AGaramond,66,9,0,0,0;3,14,10,AGaramond,0,12,0,0,0; :[font = text; inactive; preserveAspect] This is called a recursive relationship between these sequences, because the values of I, F, and C on day n depend on the value of I the preceding day. We will make two more assumptions that can be written mathematically as recursive relationships. First, we will assume that the new cases of the disease on a day are jointly proportional to the number of people infected and the number who are uninfected the preceding day. (Jointly proportional means proportional to the product.) Mathematically, ;[s] 17:0,0;17,1;26,0;87,1;88,0;90,1;91,0;97,1;98,0;106,1;107,0;131,1;132,0;320,1;327,0;429,1;436,0;503,-1; 2:9,14,10,AGaramond,0,12,0,0,0;8,14,10,AGaramond,2,12,0,0,0; :[font = text; inactive; preserveAspect] Fn = a InÊÐ 1(1 Ð InÊÐ 1). ;[s] 13:0,1;30,2;31,0;34,1;35,0;36,1;37,2;38,3;42,0;47,1;48,2;49,3;53,0;56,-1; 4:4,14,10,AGaramond,0,12,0,0,0;4,14,10,AGaramond,2,12,0,0,0;3,11,7,AGaramond,66,9,0,0,0;2,11,7,AGaramond,64,9,0,0,0; :[font = text; inactive; preserveAspect] (The constant a is positive.) Our second assumption is that the recoveries on a day are proportional to the number of people infected the previous day. We will write this ;[s] 3:0,0;14,1;15,0;173,-1; 2:2,14,10,AGaramond,0,12,0,0,0;1,14,10,AGaramond,2,12,0,0,0; :[font = text; inactive; preserveAspect] Cn = (1 Ð b) InÊÐ 1. ;[s] 9:0,1;29,2;30,0;38,1;39,0;41,1;42,2;43,3;47,0;49,-1; 4:3,14,10,AGaramond,0,12,0,0,0;3,14,10,AGaramond,2,12,0,0,0;2,11,7,AGaramond,66,9,0,0,0;1,11,7,AGaramond,64,9,0,0,0; :[font = text; inactive; preserveAspect] b is the proportion of infected people who will not be well the next day. Substituting these two assumptions into our initial equation gives us a recursive definition for I : ;[s] 6:0,1;1,0;48,1;51,0;172,1;173,0;176,-1; 2:3,14,10,AGaramond,0,12,0,0,0;3,14,10,AGaramond,2,12,0,0,0; :[font = text; inactive; preserveAspect] In = InÊÐ 1 + a InÊÐ 1(1 Ð InÊÐ 1) Ð (1 Ð b) InÊÐ 1 = (a + b) InÊÐ 1Ð a InÊÐ 12 ;[s] 34:0,1;8,2;9,3;12,1;13,2;14,0;18,1;22,3;23,1;24,2;25,0;29,3;34,1;35,2;36,0;40,3;50,1;51,3;53,1;54,2;55,0;59,3;63,1;64,3;67,1;68,3;70,1;71,2;72,0;76,3;79,1;82,2;83,0;87,4;89,-1; 5:6,11,7,AGaramond,64,9,0,0,0;11,14,10,AGaramond,2,12,0,0,0;7,11,7,AGaramond,66,9,0,0,0;9,14,10,AGaramond,0,12,0,0,0;1,11,7,AGaramond,32,9,0,0,0; :[font = text; inactive; preserveAspect] In class, we showed that if we made the definitions k = a + b and dn = a/k In, then this definition is equivalent to a simpler one for the sequence {dn}: ;[s] 20:0,0;52,1;53,0;56,1;57,0;60,1;61,0;66,1;67,2;68,0;71,1;72,0;73,1;74,0;75,1;76,2;77,0;149,1;150,2;151,0;154,-1; 3:9,14,10,AGaramond,0,12,0,0,0;8,14,10,AGaramond,2,12,0,0,0;3,11,7,AGaramond,66,9,0,0,0; :[font = text; inactive; preserveAspect] dn = k dn Ð 1(1 Ð dn - 1). ;[s] 11:0,0;28,1;29,2;32,0;35,1;36,3;40,2;45,0;46,1;47,3;51,2;54,-1; 4:3,14,10,AGaramond,2,12,0,0,0;3,11,7,AGaramond,66,9,0,0,0;3,14,10,AGaramond,0,12,0,0,0;2,11,7,AGaramond,64,9,0,0,0; :[font = text; inactive; preserveAspect] Notice that by definition, 0 < k ² a, so 0 < dn ² In ² 1. ;[s] 11:0,0;31,1;32,0;35,1;36,0;45,1;46,2;47,0;50,1;51,2;52,0;58,-1; 3:5,14,10,AGaramond,0,12,0,0,0;4,14,10,AGaramond,2,12,0,0,0;2,11,7,AGaramond,66,9,0,0,0; :[font = section; inactive; preserveAspect] Part 1: Controlling the Epidemic :[font = text; inactive; preserveAspect] Explain in a couple of sentences why the following statement is true: If the disease is mild, that is, it is not very contagious and it is easy to recover from, then the value of k will be small. ;[s] 3:0,0;180,1;181,0;197,-1; 2:2,14,10,AGaramond,0,12,0,0,0;1,14,10,AGaramond,2,12,0,0,0; :[font = text; inactive; preserveAspect] :[font = text; inactive; preserveAspect] Let's begin by assuming the disease is mild and plot the growth of {dn}. Suppose that k = 0.9 and d1 = 0.01. (If you wish, assume that a = 0.9 and b = 0 for concreteness; then dn = In.) The following command will define the appropriate sequence (and tell Mathematica to remember the values of terms it has already computed) ;[s] 22:0,0;68,1;69,3;70,0;87,1;88,0;99,1;100,2;101,0;137,1;138,0;149,1;150,0;179,1;180,3;181,0;184,1;185,3;186,1;187,0;259,1;270,0;328,-1; 4:9,14,10,AGaramond,0,12,0,0,0;9,14,10,AGaramond,2,12,0,0,0;1,11,7,AGaramond,64,9,0,0,0;3,11,7,AGaramond,66,9,0,0,0; :[font = input; preserveAspect] k = 0.9; d[1] = 0.01; d[n_] := d[n] = k d[n-1](1 - d[n-1]) :[font = text; inactive; preserveAspect] We can plot the graph of d as we learned in the lab before break, using the command ;[s] 3:0,0;25,1;26,0;84,-1; 2:2,14,10,AGaramond,0,12,0,0,0;1,14,10,AGaramond,2,12,0,0,0; :[font = input; preserveAspect] ListPlot[Table[{n, d[n]}, {n,1,20}]] :[font = text; inactive; preserveAspect] What happens to the number of people infected in the first twenty days of the epidemic? What will happen to the proportion of people infected over the long run? (You may plot more terms if you wish, to support your answer.) :[font = text; inactive; preserveAspect] :[font = text; inactive; preserveAspect] Now consider the case when the disease is a little less mild; take k = 1.1. In order to define a new sequence to examine this case, first we must clear the old definition with the command ;[s] 3:0,0;67,1;68,0;189,-1; 2:2,14,10,AGaramond,0,12,0,0,0;1,14,10,AGaramond,2,12,0,0,0; :[font = input; preserveAspect] Clear[d,k] :[font = text; inactive; preserveAspect] Now modify the original definition to treat this case (assuming d1 is still 0.01). ;[s] 4:0,0;64,1;65,2;66,0;83,-1; 3:2,14,10,AGaramond,0,12,0,0,0;1,14,10,AGaramond,2,12,0,0,0;1,11,7,AGaramond,64,9,0,0,0; :[font = text; inactive; preserveAspect] Now what happens to the proportion of people infected over the first twenty days? Over the long run? :[font = text; inactive; preserveAspect] :[font = section; inactive; preserveAspect] Part 2: Analysis of the Model using Limits :[font = text; inactive; preserveAspect] Show that if the sequence {dn} has a limit L, then the limit must satisfy the equation: ;[s] 6:0,0;27,1;28,2;29,0;43,1;44,0;88,-1; 3:3,14,10,AGaramond,0,12,0,0,0;2,14,10,AGaramond,2,12,0,0,0;1,11,7,AGaramond,66,9,0,0,0; :[font = text; inactive; preserveAspect] L = k L (1 Ð L). ;[s] 6:0,1;36,0;39,1;42,0;48,1;49,0;52,-1; 2:3,14,10,AGaramond,0,12,0,0,0;3,14,10,AGaramond,2,12,0,0,0; :[font = text; inactive; preserveAspect] Hint: Take the limit of the recursive definition as n goes to infinity. ;[s] 4:0,1;5,0;53,1;54,0;73,-1; 2:2,14,10,AGaramond,0,12,0,0,0;2,14,10,AGaramond,2,12,0,0,0; :[font = text; inactive; preserveAspect] :[font = text; inactive; preserveAspect] Use the equation you just derived to find the limit of the sequence in each of the two cases considered in Part 1. You may use Mathematica, but write your answers below. ;[s] 3:0,0;128,1;139,0;171,-1; 2:2,14,10,AGaramond,0,12,0,0,0;1,14,10,AGaramond,2,12,0,0,0; :[font = text; inactive; preserveAspect] :[font = text; inactive; preserveAspect] Public health officials often have some control over the parameters a and b, and hence over k. For instance, they may be able to decrease b by giving proper treatment and they might be able to decrease a through quarantines, or education. For what values of k will the disease tend to be eliminated from the population? ;[s] 13:0,0;68,1;69,0;74,1;75,0;92,1;93,0;139,1;140,0;203,1;204,0;260,1;261,0;322,-1; 2:7,14,10,AGaramond,0,12,0,0,0;6,14,10,AGaramond,2,12,0,0,0; :[font = section; inactive; preserveAspect] Part 3: Dependence on parameters :[font = text; inactive; preserveAspect] We have seen that given two constants: the value of k, and the initial value of d, we can find a sequence to model an epidemic. We call these constants parameters. In this part we will examine how the evolution of an epidemic depends on these parameters. ;[s] 7:0,0;53,1;54,0;81,1;82,0;154,1;164,0;257,-1; 2:4,14,10,AGaramond,0,12,0,0,0;3,14,10,AGaramond,2,12,0,0,0; :[font = text; inactive; preserveAspect] Assuming d1 is still 0.01, investigate the long term effect of an epidemic on a population with each of the following k values: 2.8, 3.84, 4.0. In each case, write a few sentences describing your conclusion. ;[s] 6:0,0;9,1;10,2;11,0;118,1;119,0;209,-1; 3:3,14,10,AGaramond,0,12,0,0,0;2,14,10,AGaramond,2,12,0,0,0;1,11,7,AGaramond,64,9,0,0,0; :[font = text; inactive; preserveAspect] :[font = text; inactive; preserveAspect] Suppose that in the case when k = 4.0, the sequence had a limit. What would it be? ;[s] 3:0,0;30,1;31,0;84,-1; 2:2,14,10,AGaramond,0,12,0,0,0;1,14,10,AGaramond,2,12,0,0,0; :[font = text; inactive; preserveAspect] :[font = text; inactive; preserveAspect] Finally, consider the effect of changing the initial value d1. Define a sequence using k = 4.0 and d1 = 0.75, and plot its graph. Repeat this using k = 4.0 and d1 = 0.7500001. Compare the long term behavior of these two sequences. ;[s] 14:0,0;59,1;60,2;61,0;88,1;89,0;100,1;101,2;102,0;150,1;151,0;162,1;163,2;164,0;234,-1; 3:6,14,10,AGaramond,0,12,0,0,0;5,14,10,AGaramond,2,12,0,0,0;3,11,7,AGaramond,64,9,0,0,0; :[font = text; inactive; preserveAspect] :[font = text; inactive; preserveAspect] Suppose a public health official said "The disease is not at all mild, the a value is about 3.1 and the b value is 0.9, but we appear to have the epidemic under control. The proportion of the population that is infected has been staying near 58% for a few days now." What would you reply? Should the official expect the epidemic to remain under control for long? What can you tell her concerning the future of this epidemic? What must be done to eradicate the disease? ;[s] 5:0,0;75,1;76,0;104,1;105,0;473,-1; 2:3,14,10,AGaramond,0,12,0,0,0;2,14,10,AGaramond,2,12,0,0,0; ^*)