# CoVID-19 SEIR model implemented in HP 35S RPN
C001 LBL C
C002 INPUT A #S0
C003 INPUT B #E0
C004 INPUT C #I0
C005 INPUT D #R0
C006 INPUT E #alpha
C007 INPUT F #beta
C008 INPUT G #gamma
C009 INPUT H #rho
C010 INPUT I #dt
C011 INPUT J #tmax
# get nsteps from tmax over dt
C012 RCL J
C013 RCL I
C014 /
C015 STO K #nsteps
# calculate dS
C016 RCL H
C017 RCL F
C018 *
C019 RCL A
C020 RCL C
C021 *
C022 RCL I
C023 *
C024 *
C025 +/-
C026 STO L #dS
# calculate dE
C027 RCL H
C028 RCL F
C029 *
C030 RCL A
C031 RCL C
C032 *
C033 *
C034 RCL E
C035 RCL B
C036 *
C037 -
C038 RCL I
C039 *
C040 STO M #dE
# calculate dI
C041 RCL E
C042 RCL B
C043 *
C044 RCL G
C045 RCL C
C046 *
C047 -
C048 RCL I
C049 *
C050 STO N #dI
# calculate dR
C051 RCL G
C052 RCL C
C053 *
C054 RCL I
C055 *
C056 STO O #dR
# calculate new S, E, I, and R
C057 RCL A
C058 RCL L
C059 +
C060 STO A #new S
C061 RCL B
C062 RCL M
C063 +
C064 STO B #new E
C065 RCL C 
C066 RCL N 
C067 +
C068 STO C #new I
C069 RCL D
C070 RCL O
C071 +
C072 STO D #new R
C073 DSE K
C074 GTO C016
C075 RTN