3 heat1n-i-glsc.c

heat1n-i-glsc.c は、同次Neumann境界条件 ( $u_x(0,t)=u_x(1,t)=0$) のもとでの1次元熱方程式を $\theta$法という陰解法 (implicit shceme) で解くためのプログラムである。

桂田 [4] の 1.4 節第3項を見ながら解読すると良い。

使い方はほぼ heat1d-i-glsc.c と同じである。

図 3: heat1n-i-glsc の実行結果 (nfunc$=1$, $N=100$, $\lambda =1/2$, $T_{\text {max}}=1$)