2.1 熱方程式の差分法プログラム(GLSC)を1つ試す

ターミナルでやってみて下さい (GLSC版)

curl -O http://nalab.mind.meiji.ac.jp/~mk/program/fdm/heat1d-e-glsc.c
ls -l heat1d-e-glsc.c
cglsc heat1d-e-glsc.c
./heat1d-e-glsc

1. 区間の分割数 N = と出るので 40[enter] と入力する。
2. λ (=τ/h^2) =と出るので 0.5[enter] と入力する。
3. 最終時刻 Tmax = と出るので 0.5[enter] と入力する。
4. 数秒で表示が終わるはず。 終りました。X の場合はウィンドウをクリックして下さい。 と表示される。 グラフの表示されているウィンドウをクリックするとプログラムが終了する。

プログラムの解読はそのうちやってもらうが、 一言で説明すると、前回配ったプリント 「発展系の数値解析」 の3節の差分方程式をそのままプログラムにしたものである。