3
所以我試圖解決三個ODE的系統,並開發了下面的代碼來解決它們使用ODEint。但是當我運行時,ODEint有一個問題要求我的函數用於方程組。ODEINT無法調用功能
from scipy.integrate import odeint
#initial parameters
X0 = 75.
M0 = 150.
G0 = 105.
N0 = 80.
T0 = 7.
u10 = 0.03
u20 = 0.0078947020
u30 = 0.0010708464
Rxg = 1.92
Rxm = 3.84
Rxma = 5.76
#define system of equations
def sugar(t,y,u10,u20,u30,Rxg,Rxm,Rxma):
ddt = [-u10*(X0+(Rxg*(G0-y[0]))+(Rxm*(M0-y[1]))+(Rxma*(N0-y[2]))),
-u20*(X0+(Rxg*(G0-y[0]))+(Rxm*(M0-y[1]))+(Rxma*(N0-y[2]))),
-u30*(X0+(Rxg*(G0-y[0]))+(Rxm*(M0-y[1]))+(Rxma*(N0-y[2])))]
return(ddt)
#time parameter
tinit = 0.
tend = 10.
h = 2.
t = arange(tinit,tend+h,h)
#initial guess
y0 = [G0,M0,N0]
#run odeint
y = odeint(sugar,y0,t,args = (u10,u20,u30,Rxg,Rxm,Rxma))
print 't =',t
print 'y =',y
而且在我的輸出獲得:
t = [ 0. 2. 4. 6. 8. 10.]
y = [[ 105. 150. 80.]
[ 105. 150. 80.]
[ 105. 150. 80.]
[ 105. 150. 80.]
[ 105. 150. 80.]
[ 105. 150. 80.]]
error: Error occurred while calling the Python function named sugar
更改Y0和T仍然導致同樣的錯誤的順序。我不明白這個錯誤,也沒有指導我解決這個問題。任何幫助將不勝感激。
更改'y0'和't'的順序仍導致相同的錯誤。 – user2855612