コロナのSIRモデル（Rの変化を検証）

SIRモデルによる感染者数の変化

date of the change

1/γ(averate recovery dates)

total dates

R:再生産率, R_0:Rの初期値, R_1:Rの新しい値.

initial value of infected=

途中でRが変わった時の感染の推移です．SIRモデルを使っています．

a = 44 (slider from 1 to 200, step 1)
R_1 = 0.5 (slider from 0.1 to 5)
R_0 = 2.5 (slider from 0.1 to 5, step 0.1)
g(t) = (R_0 + R_1) / 2 + ((R_0 - R_1) / 2 * sgn(a - t))
γ' = 10.1 (slider from 1 to 20)
γ: 1 / γ'
y'(t, x, y, z) = (γ * x)
x'(t, x, y, z) = (((g(t) * γ) * z) * x) - (γ * x)
z'(t, x, y, z) = ((((-g(t)) * γ) * x) * z)
n = 200 (slider from 100 to 600, step 50)
f(x) = If[0 ≤ x ≤ n, 1]
ratio: R_1 / R_0
i0: 1.0E-4
InputBox1 = Textfield(i0)
infected, recovered, susceptible = NSolveODE({x', y', z'}, 0, {i0, 0, 1 - i0}, n)
text1: "{R_1}/{R_0}=" + (LaTeX[ratio]) + ""