2 試してみよう1

• nfunc(1..5)= に対して、1から5までの整数値を入力できるが、 これにより初期条件 $f$ を選ぶようになっている。 それぞれどういう関数であるか読み取ること。
• $\theta=0$ とすると陽解法 (前進 Euler 法) になる。 安定性の条件 $0<\lambda\le 1/2$ について調べよ。
• 安定性の条件が成り立っているとき、 ある程度以上 $N$ を大きくすると、 ほとんど結果が変わらないことを確かめよ。
• $f(x)=\sin \pi x$ とするとき、 $u(x,t)=e^{-\pi^2 t}\sin \pi x$ であることを利用して、 $N$ が大きくなるとともに誤差がどのように減少するか調べよ。 (前回配布した資料にある実験を追試せよ。)
• 今日中に一回は図を印刷してみること。