(41C/CV/CX) Epidemiologic Programs for Computers & Calculators

[Thomas Klemm]

REED

This is the program REED for the HP-41C:

Code:

01*LBL "REED"
02 "REED-FROST"
03 ├" EPIDEMICS"
04 AVIEW
05 "<DETERMINISTIC>"
06 AVIEW
07 ADV
08*LBL A
09 1 E-1
10 STO 04
11 CLX
12 STO 05
13 STO 06
14 CF 00
15 FIX 3
16 1
17 TONE 9
18 "P? ="
19 PROMPT
20 -
21 STO 00
22 ARCL L
23 FS? 55
24 PRA
25 CF 29
26 FIX 0
27 TONE 9
28 "C0? ="
29 PROMPT
30 STO 02
31 ARCL X
32 FS? 55
33 PRA
34 TONE 9
35 "S0? ="
36 PROMPT
37 STO 03
38 ARCL X
39 FS? 55
40 PRA
41 TONE 9
42 "S. INFLUX? ="
43 PROMPT
44 STO 01
45 ARCL X
46 FS? 55
47 PRA
48*LBL 00
49 TONE 9
50 "INF. TIME? ="
51 PROMPT
52 ARCL X
53 FS? 55
54 PRA
55 ADV
56 2
57 X=Y?
58 GTO 01
59 CLX
60 1
61 X#Y?
62 GTO 00
63 SF 00
64*LBL 01
65 ISG 04
66 GTO 02
67 TONE 9
68 STOP
69*LBL 02
70 RCL 00
71 RCL 02
72 FC? 00
73 XEQ 03
74 Y^X
75 1
76 -
77 CHS
78 RCL 03
79 *
80 RND
81 "T="
82 ARCL 04
83 ├"  C="
84 ARCL X
85 RCL 06
86 RCL 02
87 FS? 00
88 ST+ 05
89 STO 06
90 RDN
91 FC? 00
92 ST+ 05
93 RDN
94 ST- 03
95 STO 02
96 RCL 01
97 ST+ 03
98 ├"  S="
99 ARCL 03
100 X#0?
101 GTO 04
102 ├"  I="
103 ARCL 05
104*LBL 04
105 FC? 55
106 PROMPT
107 FS? 55
108 PRA
109 RCL 02
110 X=0?
111 GTO 05
112 GTO 01
113*LBL 03
114 RCL 06
115 +
116 RTN
117*LBL 05
118 TONE 9
119 "END OF EPIDEMIC"
120 AVIEW
121 ADV
122 RTN
123*LBL B
124 CF 00
125 SF 29
126 FIX 2
127 END

And here the same program for the HP-42S:

Code:

00 { 302-Byte Prgm }
01▸LBL "REED"
02 "REED-FROST"
03 ├" EPIDEMICS"
04 AVIEW
05 "<DETERMINISTIC>"
06 AVIEW
07 ADV
08▸LBL A
09 1ᴇ-1
10 STO 04
11 CLX
12 STO 05
13 STO 06
14 CF 00
15 FIX 03
16 1
17 TONE 9
18 "P? ="
19 PROMPT
20 -
21 STO 00
22 ARCL ST L
23 FS? 55
24 PRA
25 CF 29
26 FIX 00
27 TONE 9
28 "C0? ="
29 PROMPT
30 STO 02
31 ARCL ST X
32 FS? 55
33 PRA
34 TONE 9
35 "S0? ="
36 PROMPT
37 STO 03
38 ARCL ST X
39 FS? 55
40 PRA
41 TONE 9
42 "S. INFLUX? ="
43 PROMPT
44 STO 01
45 ARCL ST X
46 FS? 55
47 PRA
48▸LBL 00
49 TONE 9
50 "INF. TIME? ="
51 PROMPT
52 ARCL ST X
53 FS? 55
54 PRA
55 ADV
56 2
57 X=Y?
58 GTO 01
59 CLX
60 1
61 X≠Y?
62 GTO 00
63 SF 00
64▸LBL 01
65 ISG 04
66 GTO 02
67 TONE 9
68 STOP
69▸LBL 02
70 RCL 00
71 RCL 02
72 FC? 00
73 XEQ 03
74 Y↑X
75 1
76 -
77 +/-
78 RCL 03
79 ×
80 RND
81 "T="
82 ARCL 04
83 ├"  C="
84 ARCL ST X
85 RCL 06
86 RCL 02
87 FS? 00
88 STO+ 05
89 STO 06
90 R↓
91 FC? 00
92 STO+ 05
93 R↓
94 STO- 03
95 STO 02
96 RCL 01
97 STO+ 03
98 ├"  S="
99 ARCL 03
100 X≠0?
101 GTO 04
102 ├"  I="
103 ARCL 05
104▸LBL 04
105 FC? 55
106 PROMPT
107 FS? 55
108 PRA
109 RCL 02
110 X=0?
111 GTO 05
112 GTO 01
113▸LBL 03
114 RCL 06
115 +
116 RTN
117▸LBL 05
118 TONE 9
119 "END OF EPIDEMIC"
120 AVIEW
121 ADV
122 RTN
123▸LBL B
124 CF 00
125 SF 29
126 FIX 02
127 END

This allows us to reproduce FIGURE 3:

Quote:Three examples of printed output from program "REED" which is based on the deterministic formulation. The three epidemics displayed differed in the magnitude of the p values: 0.02, 0.03, and 0.04, respectively.

---

REED-FROST EPIDEMICS
<DETERMINISTIC>

P? =0.020
C0? =1
S0? =100
S. INFLUX? =0
INF. TIME? =1

T=1 C=2 S=98 I=1
T=2 C=4 S=94 I=3
T=3 C=7 S=87 I=7
T=4 C=11 S=76 I=14
T=5 C=15 S=61 I=25
T=6 C=16 S=45 I=40
T=7 C=12 S=33 I=56
T=8 C=7 S=26 I=68
T=9 C=3 S=23 I=75
T=10 C=1 S=22 I=78
T=11 C=0 S=22 I=79
END OF EPIDEMIC

---

REED-FROST EPIDEMICS
<DETERMINISTIC>

P? =0.030
C0? =1
S0? =100
S. INFLUX? =0
INF. TIME? =1

T=1 C=3 S=97 I=1
T=2 C=8 S=89 I=4
T=3 C=19 S=70 I=12
T=4 C=31 S=39 I=31
T=5 C=24 S=15 I=62
T=6 C=8 S=7 I=86
T=7 C=2 S=5 I=94
T=8 C=0 S=5 I=96
END OF EPIDEMIC

---

REED-FROST EPIDEMICS
<DETERMINISTIC>

P? =0.040
C0? =1
S0? =100
S. INFLUX? =0
INF. TIME? =1

T=1 C=4 S=96 I=1
T=2 C=14 S=82 I=5
T=3 C=36 S=46 I=19
T=4 C=35 S=11 I=55
T=5 C=8 S=3 I=90
T=6 C=1 S=2 I=98
T=7 C=0 S=2 I=99
END OF EPIDEMIC

---

FROST

This is the program FROST for the HP-41C:

Code:

01*LBL "FROST"
02 "REED-FROST"
03 ├" EPIDEMICS"
04 AVIEW
05 "<STOCHASTIC>"
06 ├"<TROUGH>"
07 AVIEW
08 ADV
09*LBL A
10 1 E-1
11 STO 13
12 CLX
13 STO 07
14 STO 09
15 STO 10
16 STO 14
17 CF 00
18 CF 01
19 CF 02
20 CF 03
21 FIX 3
22 2
23 TONE 9
24 "P? ="
25 PROMPT
26 ARCL X
27 FS? 55
28 PRA
29 CF 29
30 FIX 0
31 /
32 2
33 -
34 RND
35 STO 04
36 TONE 9
37 "C0? ="
38 PROMPT
39 STO 06
40 ARCL X
41 FS? 55
42 PRA
43 TONE 9
44 "S0? ="
45 PROMPT
46 STO 05
47 ARCL X
48 FS? 55
49 PRA
50 TONE 9
51 "S. INFLUX? ="
52 PROMPT
53 STO 12
54 ARCL X
55 X#0?
56 SF 03
57 FS? 55
58 PRA
59 TONE 9
60*LBL 04
61 "INF. TIME? ="
62 PROMPT
63 ARCL X
64 FS? 55
65 PRA
66 1
67 X=Y?
68 GTO 05
69 CLX
70 2
71 X#Y?
72 GTO 04
73 SF 02
74*LBL 05
75 TONE 9
76 "SEED<=9999?="
77 PROMPT
78 ARCL X
79 FS? 55
80 PRA
81 ADV
82 1 E4
83 /
84 STO 08
85*LBL 06
86 FS? 49
87 OFF
88 "<"
89 FS? 55
90 ACA
91 RCL 07
92 RCL 06
93 RCL 05
94 RCL 04
95 STO 00
96 +
97 STO 01
98 +
99 STO 02
100 FC? 03
101 +
102 STO 03
103 ISG 13
104 GTO 08
105 BEEP
106 STOP
107*LBL 08
108 FC? 02
109 GTO 07
110 RCL 14
111 ST+ 02
112 ST+ 03
113*LBL 07
114 .003
115 STO 11
116 RCL 08
117 9821
118 *
119 .211327
120 +
121 FRC
122 STO 08
123 RCL 03
124 *
125 INT
126 1
127 +
128*LBL 09
129 RCL IND 11
130 X<>Y
131 X<=Y?
132 GTO IND 11
133 ISG 11
134 GTO 09
135 ">"
136 FS? 55
137 ACA
138 GTO 13
139*LBL 00
140 "*"
141 SF 01
142 1
143 GTO 10
144*LBL 01
145 "S"
146 1
147 ST+ 10
148 GTO 10
149*LBL 02
150 "C"
151 SF 00
152 1
153 GTO 10
154*LBL 03
155 "I"
156 1
157*LBL 10
158 ST- IND 11
159 ISG 11
160 GTO 10
161 FS? 55
162 ACA
163*LBL 11
164 FC?C 01
165 GTO 07
166*LBL 13
167 FC?C 00
168 GTO 12
169 RCL 10
170 ST+ 09
171*LBL 12
172 0
173 STO 10
174 RCL 01
175 RCL 02
176 X>Y?
177 GTO 07
178 RCL 09
179 RCL 14
180 RCL 06
181 FC? 02
182 ST+ 07
183 STO 14
184 RDN
185 FS? 02
186 ST+ 07
187 RDN
188 ST- 05
189 STO 06
190 RCL 12
191 ST+ 05
192 0
193 STO 09
194 BEEP
195 FS? 55
196 PRBUF
197 ADV
198 "T="
199 ARCL 13
200 ├"  C="
201 ARCL 06
202 ├"  S="
203 ARCL 05
204 FS? 03
205 GTO 14
206 ├"  I="
207 ARCL 07
208*LBL 14
209 AVIEW
210 FC? 55
211 STOP
212 CLD
213 ADV
214 RCL 06
215 X>0?
216 RCL 05
217 X>0?
218 GTO 06
219 TONE 9
220 "END OF EPIDEMIC"
221 AVIEW
222 ADV
223 RTN
224*LBL B
225 CF 00
226 CF 01
227 CF 02
228 CF 03
229 SF 29
230 FIX 3
231 END

It uses the printer's ACA and PRBUF commands, which are not supported by the HP-42S.
I haven't tested it yet, but I hope I didn't make any typos when transcribing from the original PDF.

These programs might have been useful four years ago.
We now also have other tools for modeling epidemics:




As always, thank you for your contributions.