私は次のような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