#stability

陰解法とトーマス法 — 1次元拡散を無条件安定で解くひと筋のトリック

rが0.6のとき陽解法は60ステップで1e12まで跳ねる。陰解法は三重対角を1回掃くだけで、同じ刻みでもr=2でも有限のまま残る。