« WiiのビリーズブートキャンプのCMがセイバンの「天使のはね」だというのは度肝を抜かれた! | トップページ | KdV方程式をDOP853(ルンゲクッタ8次)でExcel VBAで計算 »

2011年5月18日 (水)

拡散方程式をDOP853(ドルマンプリンス-ルンゲクッタ8次)で計算(Excel VBA)

ルンゲクッタ8次のルーチン、DOP853は偏微分方程式にも適用できる。まずは拡散方程式でもやってみよう。

∂u(x,t)/∂t = D ∆u(x,t)

du_i/dt = D (u_i+1 + u_i-1 - 2*u_i) / h^2

のように空間部分を差分化すると、これは普通の常微分方程式になった。でそのままDOP853で計算すると、

Diffdop853

と一瞬で計算できた。ではこれから偏微分方程式シリーズを少しやってみよう(続く)。

« WiiのビリーズブートキャンプのCMがセイバンの「天使のはね」だというのは度肝を抜かれた! | トップページ | KdV方程式をDOP853(ルンゲクッタ8次)でExcel VBAで計算 »

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

コメント

コメントを書く

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

トラックバック

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

この記事へのトラックバック一覧です: 拡散方程式をDOP853(ドルマンプリンス-ルンゲクッタ8次)で計算(Excel VBA):

« WiiのビリーズブートキャンプのCMがセイバンの「天使のはね」だというのは度肝を抜かれた! | トップページ | KdV方程式をDOP853(ルンゲクッタ8次)で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
フォト
無料ブログはココログ