SciPyのODEINT関数の配列入力

Question

配列をODEに渡して、ODEINT関数を使用してそのODEを解決しようとしています。しかし、私は次のエラーを取得する：ここで

RuntimeError: The size of the array returned by func (50) does not match the size of y0 (1)

は私のサンプルコードは次のとおりです。

from scipy.integrate import odeint 
import numpy as np 
import matplotlib.pyplot as plt 


#First ODE 
def eq1 (x,t): 
    return (0.5*(x)+2*x**2) 

#Second ODE 
def eq2 (y,t,m): 
    p = np.sqrt(y)+m*y**2 
    return p 

t = np.linspace(0,1,50) 
IC = [2] # Initial Condition 
m = [0.1]*len(t) # A random variable 

#Solver for the first ODE 
one = odeint(eq1,IC,t) 
plt.plot (t,one) 

#Solver for the Second ODE 
two = odeint (eq2,IC,t,(m,)) 
plt.plot (t,two) 

最初のソルバーは正常に動作しますが、もう一つはエラーをスローします。このエラーメッセージに基づいて、私はICの次元が変数mの次元と同じではないことを理解しています。したがって、最初のODEのソルバーに関連する行をコメントアウトし、ICをIC = [2]*len(m)と変更すると、エラーメッセージが表示されなくなります。方法3、two.shapeは（50,50）です。私は解のどの次元（2つ）が実際の解答なのか分かりません。

私は本当にいくつかのガイダンスに感謝します。 ありがとうございました！

Answer 1

私の知る限り、配列をPythonのODE関数のパラメータとして渡す方法はありません。 Funnily十分に私はそれがMatlabで可能だと思います。あなたの問題を解決する唯一の方法は（しかし非効率的に）FORループを使うことです。