« 京都初詣シリーズ2013(その8 猿田彦神社) | トップページ | 奈良初詣シリーズ2013(その1 東大寺→奈良の大仏) »

2013年1月15日 (火)

KdV方程式をExcel VBAで計算してGIFアニメにしてみた。

最近、GIFアニメに凝っている。方程式の解が目に見える形にできるのがいいなあ、と。

今回はKdV方程式

∂u/∂t +u*∂u/∂x+δ^3 ∂^3u/∂x^3 = 0

を計算。

方法は、空間微分をO(h^6)という高次の差分で近似、かつ時間微分をルンゲクッタ8次のスキームであるDormand-Prince法で計算している。

(内容としてはオリジナルのザブスキーとクルスカルのものと同じ)。

Kdv


ソリトンが複数衝突を繰り返す様子がよくわかる。

« 京都初詣シリーズ2013(その8 猿田彦神社) | トップページ | 奈良初詣シリーズ2013(その1 東大寺→奈良の大仏) »

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

コメント

コメントを書く

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

トラックバック

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

この記事へのトラックバック一覧です: KdV方程式をExcel VBAで計算してGIFアニメにしてみた。:

« 京都初詣シリーズ2013(その8 猿田彦神社) | トップページ | 奈良初詣シリーズ2013(その1 東大寺→奈良の大仏) »

最近の記事

最近のコメント

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