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

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

今回はKdV方程式

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

を計算。

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

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

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