ブラッセレータ（Brusselator)をドルマン・プリンスで計算（Excel VBA)

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

\[ \dot y_{1} = 1+y_{1}^{2}y_{2}-4y_{1} \] \[ \dot y_{2} = 3y_{1} - y_{1}^{2}y_{2} \]

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

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