源代码参考 defaults.py，自然状态下(不封城)
r = 20.
beta = 0.048
h = 0.048
theta = 0.1
gamma = 0.1
封城后，r = 3.
人口迁移引力模型中，迁移系数 K = 0.04，迁移指数 alpha = 2
封城后，城市间无人口迁移。
以上数据中，r = 20 以及人口迁移引力模型中的 K = 0.04 源于猜测，K 值猜测过程详见 README.md 文件，其余数据来源在张同学的参数化报告中介绍 supplement/张_模型参数化.docx
城市布局 5 * 5 矩阵，人口均为 1000 万人，城市间距 100 km，初始时向中心城市 city22 投放一位潜伏者。分别模拟以下 4 种封城时机：
| 序号 | 封城时间 | 结果文件目录 |
|---|---|---|
| 0 | 5 | reports/sim_0 |
| 1 | 10 | reports/sim_1 |
| 2 | 15 | reports/sim_2 |
| 3 | 10000(不封城) | reports/sim_3 |
模拟时间区间 [0, 360] (一年)，从第 0 天开始，模拟时间步长与采样步长均为 0.1 天，模拟使用 RK4 算法。
相应结果文件目录中包含：
- 国家整体变化曲线
country.svg - 各城市变化曲线
cityxx.svg - 数据报告，包含三项补充信息：感染最高点时间，感染最高点，易感人群下降至初始的 50% 耗时
report.csv - 模拟动画
animation.mp4，描述各城市感染人群随时间变化情况。圆圈面积与感染人数成正比。
程序和说明文件中，人口的单位均为 万人 ，距离单位均为 km，时间单位均为 天。