オレゴネータ(Oregonator)をドルマン・プリンス法で計算（Excel VBA)

化学反応のモデルとして有名なオレゴネータをルンゲクッタ8次であるDormand&Prince法のコード、DOP853のVBA移植版で計算しよう。式はこちら。

\[ \dot y_{1} = 77.27 (y_{2} + y_{1}(1-8.375\times10^{-6}y_{1}-y_{2})) \] \[ \dot y_{2} = \frac{1}{77.27} (y_{3} - (1+y_{1})y_{2}) \] \[ \dot y_{3} = 0.161 (y_{1} - y_{3}) \]

常微分方程式の数値解法I,IIに紹介されていたパラメータを使った。

本に載っていたようなグラフが再現できた。