20/04/04 11:37:34 ZFu90Xbq.net
SEIR MODEL
dS(t)/dt = mu*(N-S) - b*S(t)*I(t)/N - nu*S(t)
dE(t)/dt = b*S(t)I(t)/N - (mu+sig)*E(t)
dI(t)/dt = sig*E(t) - (mu+g)*I(t)
dR(t)/dt = g*I(t) - mu*R + nu*S(t)
mu:自然死亡率 b:感染率(S->I)
nu:ワクチン有効率(S->R) sig:発症率(E->I),g:回復率(I->R)
の微分方程式の数値解を使ってシミュレーション
対策しない（外出を控えず、マスクもしない）方が患者や死者は増えるけど早く収束するな。
contacｔ_rate と trannsmission_probabilityを変化させてグラフにしてみた。
URLﾘﾝｸ(i.imgur.com)