« シーサイド高浜でお土産を買ってきた。 | トップページ | グリーン車に初めて乗った。 »

2011年5月 7日 (土)

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

Bruss01

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

« シーサイド高浜でお土産を買ってきた。 | トップページ | グリーン車に初めて乗った。 »

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

コメント

コメントを書く

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

トラックバック

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

この記事へのトラックバック一覧です: ブラッセレータ(Brusselator)をドルマン・プリンスで計算(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
フォト
無料ブログはココログ