Excelで双曲型偏微分方程式を解いてみる（その6） - 三次精度(Kawamura-Kuwahara)風上差分編

∂ｆ/∂ｔ + c∂f/∂x=0を解こうシリーズその6。3次精度(Kawamura-Kuwahara)は...

fn+1(i) = fn(i) - (c / 12) * (-fn(i + 2) + 8 * fn(i + 1) - 8 * fn(i - 1) + fn(i - 2))
            - (|c| / 4) * (fn(i + 2) - 4 * fn(i + 1) + 6 * fn(i) - 4 * fn(i - 1) + fn(i - 2))

とかなり長くなってきた。

やっぱり少しなまるけど、これが一番応用がききそうかな。