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

#?N為需要傳播的人群總數，這里可以修改為變化的函數哦
N?=?10000
#?β為傳染率系數，也就是傳染的概率，復雜網絡中一般為固定的數值，社交網絡或者其他網絡中應該是變化的函數哦
beta?=?0.25
#?gamma為恢復率系數，也就是用戶對于某個主題的接收概率，應該針對不同的用戶興趣進行相關的函數調整。
gamma?=?0
#?I_0為感染者的初始人數，可以設置為社區中用戶的人數，多個社區可以進行同步傳播
I_0?=?1
#?S_0為易感者的初始人數
S_0?=?N?-?I_0
#?感染者每天接觸人數
P?=?1
#?T為傳播時間
T?=?150

#?INI為初始狀態下的數組
INI?=?（S_0I_0）


def?funcSI（inivalue_）:
????Y?=?np.zeros（2）
????X?=