在Python中解决多维ODE

Question

我必须解决二阶数值的ODE。问题是我需要同时为两个对象（通常是N个对象）解决它，每个对象有4个变量（2个维度）：2个位置，2个速度。请参阅下面的代码以查看示例。

当尝试使用提供的方法在Python（Scipy）中解决它时，我总是遇到错误 - ＆＃34;引发ValueError（＆＃34; y0必须是1维的。＆＃34;） ValueError：y0必须是1维的。＆＃34;

y0需要的是：

array([[[1, 0],
    [0, 1]],

   [[2, 0],
    [0, 2]]])

我认为在Matlab中可能会这样做，但我想继续使用Python。另外，我想通过Runge-Kutta模型来解决它。

可能的解决方案是将所有内容重塑为1维，重新整形并再次重新整形，但这可能不那么有效和高效。

感谢您的帮助。请问我是否有不明确的事情