k = 10
b = 0.06
q = 0
INIT infectuous_undetected = 0
INIT susceptible = 999999
INIT exposed = 1
N0 = INIT(susceptible) + INIT(exposed)
u = 0.1
INIT susceptible_quarantined = 0
sus_to_quar = q * k * (1 - b) * infectuous_undetected * susceptible / N0
quar_to_sus = u * susceptible_quarantined
sus_to_expos = k * b * (1 - q) * infectuous_undetected * susceptible / N0
sus_to_expos_quar = q * k* b * infectuous_undetected * susceptible / N0
p = 0.2
expos_to_undet = p * exposed
INIT exposed_quarantined = 0
expos_quar_to_inf_quar = p * exposed_quarantined
w = 0.0375
undet_to_iso = w * infectuous_undetected
m = 0.1
undet_to_death = m * infectuous_undetected
v = 0.0625
undet_to_im = v * infectuous_undetected
INIT infectuous_quarantined = 0
quar_to_death = m * infectuous_quarantined
inf_quar_to_iso = w * infectuous_quarantined
inf_quar_to_rec = v * infectuous_quarantined
INIT infectuous_isolated = 0
iso_to_immune = v * infectuous_isolated
iso_to_death = m * infectuous_isolated
INIT SARS_death = 0
INIT recovered_immune = 0
R = k * b * (1 - q) /(v + m + w)
infecteds_total = exposed + exposed_quarantined + infectuous_undetected + infectuous_quarantined + infectuous_isolated