1、首先是SIR传染病模型，也就是三个微分方程，如图。其中中的S表示易感者，I表示感染者，R表示移出者。a表示易感者的感染概率，b表示感染者的恢复概率,也就是移出率。
2、然后是编一个模型的代码，保存为sir.m。a和b的值可根据自己的情况设定。
3、再编写一个运行的m文件，可命名为sirrun.m，代码如下。[0,40]是t的范围，[0.04 0.95 0.01]分别是I、S和R的初始值。
4、运行的结果如下 ，其中红绿线表示S，蓝线表示I，线表示R。
5、运行完也会显示结果数据，由于数据较长，这里只截图一部分，如下：
