如何用fortran求解复矩阵微分方程？

Question

我正在寻求求解 dy/dt=c.H.y 形式的一阶复矩阵微分方程，其中 H 是具有复元素的复矩阵，c 是复数。 Python 和 cpp 提供内置库，如“odeint”、“solve ivp”和“boost odeint”，可以轻松解决这个问题。但是，我面临着在 Fortran 中寻找合适的包以处理复数并有效求解此类 ODE 的挑战。我探索过诸如“odepack”之类的选项，但它无法处理复数。

非常感谢您的帮助。

Answer 1

存在一个名为slatec的库（https://www.netlib.org/slatec/），其中CDRIV可以求解复杂的ODE。