感染症流行過程の数理モデルをPythonで解く - Qiita

注意！この記事はなんとなくやってみた以上の情報を含みません。 感染症数理モデル 昨今のCOVID19事情を受けて色々感染症の数理モデルの話は出ていると思うわけですが、それをPythonのodeintで解いて、そのままmatplotlibで可視化してみよう、という事です。 くどいようですがこの記事のなんかを適当に実装した結果からCOVID19に関して何かのお気持ちを表明しては駄目だぞ。 で、解くのは一番簡単なSIRモデルとしまして、 \begin{align} \frac{\mathrm{d} S}{\mathrm{d} t} & = - \beta I S\\ \frac{\mathrm{d} I}{\mathrm{d} t} & = \beta I S - \gamma I\\ \frac{\mathrm{d} R}{\mathrm{d} t} & = \gamma I \end{ali