添付プログラムのおかしい点を指摘してほしい。

Question

Cir 2t 2021-11-17

评论： Cir 2t 2021-11-18

以前に自分で動かしていたんですが、ひょんなことからデータを消してしまい、添付プログラムの動かし方を忘れてしまいました。

どこか変なのですが、教えていただきたいです。

HH.mでHHscheme.mを動かして微分方程式を解くのですが、下記のエラーが出ます。

エラー odearguments (line 91)

HHSCHEME は列ベクトルを返さなければなりません。

エラー ode15s (line 149)

[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...

どうかよろしくお願いいたします。

Answer 1

Atsushi Ueno 2021-11-17

HHSCHEME は列ベクトルを返さなければなりません。

を真に受けてHHscheme.mの最期を下記の様に列ベクトルにしたら動きました。

dVdt=[dV(1)';dV(2)';dV(3)';dV(4)'];

Cir 2t 2021-11-18

ありがとうございます。

助かりました。