SIR模型： import scipy.integrate import numpy as np import matplotlib.pyplot as plt # model def SIR_model(y, t, beta, gamma): S, I, R = y dS_dt = -beta * S * I dI_dt = beta * S *…