バグを直した＆メッシュ数に応じて出力画像の比率を変えるようにした。 ■出力画像 ■Rustソース use std::time::{Duration, Instant}; use plotters::prelude::*; const PLT_RATIO: usize = 20; const T_MAX: usize = 1000; const X_MAX: usize = 10; const Y…

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