移流方程式を高次ルンゲクッタで計算してみる（その2）

昨日は時間をルンゲクッタ5次、空間を対称な4次差分にしたけど矩形波の伝播はぐだぐだになった。そこで空間はやはり風上差分を使おう。3次精度風上差分を使ったときの結果がこれ。

まあまあ波形を保ったまま伝播してますな。

これつかっていろいろ非線形波動を計算してみますか。