« フィリピンの刑務所で受刑者が踊るThey don't care about usは感動的。 | トップページ | くじら雲 »

2010年7月 7日 (水)

陰的ルンゲクッタ(ラダウIIA)でOregonatorを計算(Excel VBA)

さて、スティフな(硬い)方程式の例としてまずはオレゴネータかな。8次のルンゲクッタで計算したときに、すさまじく計算しにくかったです。

http://www.scholarpedia.org/article/Oregonator

εdx/dt = q*y - x*y + x*(1-x)

ε'dx/dt = -q*y - x * y + f *z

   dz/dt = x - z

で計算結果はこちら。何も考えなくてもふつーに計算できる!!!!すごい!!!!

Oregonatorimplicit

これもソースをつけておきます。

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

« フィリピンの刑務所で受刑者が踊るThey don't care about usは感動的。 | トップページ | くじら雲 »

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

コメント

コメントを書く

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

トラックバック

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

この記事へのトラックバック一覧です: 陰的ルンゲクッタ(ラダウIIA)でOregonatorを計算(Excel VBA):

« フィリピンの刑務所で受刑者が踊るThey don't care about usは感動的。 | トップページ | くじら雲 »

最近のコメント

2018年11月
        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  
フォト
無料ブログはココログ