Octave 倍频程复微分方程组的求解

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

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

其中$\rho$和$H$是$n x n$矩阵

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

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的错误消息。那么，我的问题是：如何解一个倍频程的复微分方程？

我想问题是因为我在等号的左手边，而不是右手边

你试过了吗

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

xdot（1）=/i；
xdot（2）=/i；