蔵本・シバシンスキー方程式をExcel VBA&擬スペクトル法&高次ルンゲクッタ法で計算

昨日はKdV方程式を計算してみたけれど、今日は Kuramoto-Sivashinsky方程式をやってみよう。空間部分は擬スペクトル法、時間部分は高次ルンゲクッタ法（5次のCash-Karp法）で計算してみる。

これはさすがに時間がかかった。参考にしたのは「常微分方程式の数値解法II」に出ていた方法（コレ-エックマン-エプスタイン-シュトゥッベ）です。文献ではN=1024になっているんだが、そんなことしたらすさまじく時間がかかるのでN=200くらいで。

結果はこちら。なかなかきれいな絵がかけた。