Excelでスペクトル法でKdV方程式を。

次はKdV方程式をやってみよう。

∂u/∂t + u ∂u/∂x = μ ∂^3u/∂x^3

で、初期条件をcos(πx)として計算。時間発展は蛙跳び（leap-frog)。

μ=0.000474の時の計算結果がこれ。

 

ザブスキーとクルスカルのオリジナルの文献に近い図になっているな。