octave - 在 Octave 中求解复微分方程组

Question

好吧，我想以以下形式解决一个复杂的 ODE 系统：

$ i\hbar \dfrac{\partial \rho}{\partial t} = \left[ H, \rho \right] $

http://mathbin.heroku.com/rq65idQ

其中 $\rho$ 和 $H$ 是 $nxn$ 矩阵。

我尝试使用以下方法定义通用函数：

function xdot = f(x,t)

i*xdot(1)=
i*xdot(2)=

endfunction
x0=[0;0];
t=linspace(0,20,200);
y=lsode("f",x0,t).

但我收到有关i的错误消息。那么我的问题是：如何在 Octave 中求解一个复杂的微分方程？

Answer 1

我认为问题是因为 i 位于等号的 LHS 而不是 RHS。

你有没有尝试过？

xdot(1) = <whatever your RHS expression is> / i;
xdot(2) = <whatever your RHS expression is> / i;