求解Octave中复微分方程组的系统

Question

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

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

http://mathbin.heroku.com/rq65idQ

其中$ \ rho $和$ H $是$ n x n $ matrices。

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

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

我认为这个问题是因为我在等号的LHS而不是RHS。

你试过吗？

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