« 「神様のパズル」を読んだ。 | トップページ | 「二の悲劇」を読んだ。 »

2011年4月28日 (木)

オレゴネータ(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に紹介されていたパラメータを使った。

Oreg01

Oreg02

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

« 「神様のパズル」を読んだ。 | トップページ | 「二の悲劇」を読んだ。 »

学問・資格」カテゴリの記事

コメント

コメントを書く

(ウェブ上には掲載しません)

トラックバック

この記事のトラックバックURL:
http://app.cocolog-nifty.com/t/trackback/512682/51418408

この記事へのトラックバック一覧です: オレゴネータ(Oregonator)をドルマン・プリンス法で計算(Excel VBA):

« 「神様のパズル」を読んだ。 | トップページ | 「二の悲劇」を読んだ。 »

最近の記事

2018年6月
          1 2
3 4 5 6 7 8 9
10 11 12 13 14 15 16
17 18 19 20 21 22 23
24 25 26 27 28 29 30
フォト
無料ブログはココログ