Python:求解具有复杂初始条件的二阶微分方程

我想用类似于odeint的方法来解一个二阶变系数微分方程。这个问题是，如果初始条件很复杂（现在就是这样），它就不起作用

---

你知道用类似于odeint的东西来解上述方程的方法吗？

你总是可以使用真实的组件（

odeint

约定）

其中

f

表示二阶常微分方程的显式形式

如果我没记错的话，您可以在

scipy.integrate.ode

中使用的方法（“vzode”？）中有一种直接处理复杂状态变量。

不接受复杂变量。您可以使用：更新的解算器；较老的类与“零电压”积分器；或者，我编写的一个处理复数和矩阵值微分方程的

odeint

包装器。

“复数初始条件”-你能描述它们吗？什么意义上的复杂？它们有什么不同？自变量是什么？Initial的意思是“在时间零点”，所以我很难想象复杂性还能从哪里来。

def odesys(u,t):
    z = u[0]+1j*u[1]
    dz = u[2]+1j*u[3]
    d2z = f(t,z,dz)
    return [ dz.real, dz.imag, d2z.real, d2z.imag ]