Excelで双曲型偏微分方程式を解いてみる（その2） - 蛙跳び編

∂ｆ/∂ｔ + c∂f/∂x=0を解こうシリーズその２。FTCSはだめだということを見たが、次は蛙跳び法（leap-frog)を見てみよう。

χ=c Δｔ/Δxと置くと、

fn+1(i) = fn-1(i) + χ（fn(i+1) - fn(i-1))

となる。

これは拡散項がつくとそこそこいいはずなんだけど、純粋な双曲型ではまだだめだな。次いってみよう。