※バグあり メッシュ数を縦≠横にするとおかしくなる とりあえずソースと結果を貼っておく。 陰解法と書いたが、実際は陰解法と陽解法を平均した差分方程式を使っている。(クランクニコルソン法っぽくした) ↓の記事で導出した式と、陽解法を組み合わせている。…