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 で複雑な微分方程式を解く方法です。