#?-*-?coding:?utf-8?-*-
“““
Created?on?Tue?Apr?14?22:00:16?2020

@author:?TSON
“““

import?scipy.integrate?as?spi
import?numpy?as?np
import?matplotlib.pyplot?as?plt

#?N為人群總數
N?=?100000
#?β為傳染率系數
beta?=?0.19
#?gamma為恢復率系數
gamma?=?0.1
#?Te為疾病潛伏期
Te?=?5
#?I_0為感染者的初始人數
I_0?=?1
#?E_0為潛伏者的初始人數
E_0?=?0
#?R_0為治愈者的初始人數
R_0?=?0
#?S_0為易感者的初始人數
S_0?=?N?-?I_0?-?E_0?-?R_0
#?T為傳播時間
T?=?150
#?INI為初始狀態下的數組
INI?=?（S_0E_0I_0R_0）

def?funcSEIR（inivalue_）:
?Y?=?np.zeros（4）
?X?=?inivalue
?#?易感個體變化
?Y[0]?=?-?（beta?*?X[0]?*（?X[2]+X[1]））?/?N
?#?潛伏個體變化
?Y[1]?=?（beta?*?X[0]?*（?X[2]+X[1]））?/?N?-?X[1]?/?Te
?#?感染個體變化
?Y[2]?=?X[