陰的ルンゲクッタ法（ラダウIIA）でOregonatorを計算（Excel VBA）-その２

「常微分方程式の数値解法II」にもオレゴネータの例が出ていた。ちゃんと再現するか確認。

式は

y1 ' = 77.27 * (y2 + y1*(1-8.375E-6*y1-y2))

y2'=(y3-(1+y1)*y2)/77.27

y3' = 0.161*(y1-y3)

ソースはこちら。

「implicit-Oregonator2.txt」をダウンロード

実は！常微分方程式の数値解法IIでちょっと誤植があって、初期条件は

y1=1,y2=2,y3=3ではなく、y1=2,y2=1,y3=3ですよ。これなら本の内容を完全に再現しました。