1

私は次のようなODEを解いています

y''(t) + ay(t) + b = 0

Matlabのode45ソルバーで。エラー基準を満たす解が見つかるまで繰り返します。y比較のために各ステップの前の値を読み取りたい。

これを行う最良の方法は何ですか?

私のメインファイルは

[t,y] = ode45(@odefnc,tspan,[0 0]',options);

そして、odefnc

function dx=odefnc(t,x)

...

if history(end)<13
    dx=[x(2),-a*x(1)-b]'
else
    dx=[x(2),-c*x(1)-d]'
end


if flag==1
   history(end+1)=dx;
end
4

1 に答える 1