函数A 如下:
def SIR_model(y, t, B, r0):
S, I, R = y
dS_dt = -B * S * I
dI_dt = B * S * I - r0 * I
dR_dt = r0* I
O.append(I)
r.append(R)
s.append(S)
return ([dS_dt, dI_dt, dR_dt ])
这里的B指的是函数B中的B
S,I,R有初值,B的初值也能求出来,
S0=0.999999551536256
I0=0.000000448463744
R0 = 0.0
函数B如下:
B = B0 - C * (i+r)
B0=0.36,C=1250,i和r是从函数A中获得的
我不知道是不是该用类来写,如果用类来定义的话该怎么定义呢?谢谢!