1次元熱方程式を差分法で解いた。以下のことが分かった。
heatT.edp で の場合に、 時間刻みを大きくしても安定性が保たれるか調べる。
[ |
l]キーボード入力が出来るように改造する
real Tmax=10, tau=0.01, t, theta=1; // thetaを加える。=1の時は実は後退Euler // コメント・アウトすると、m, tau, theta が実行時に変更できるようになる // cout << "m dt theta: "; cin >> m >> tau >> theta; // cout << "m=" << m << ", tau=" << tau << ", theta=" << theta << endl; mesh Th=square(m,m);Tmax も大きめにしておくと分かりやすい。 |