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$)
Image HEAT1NI-v



桂田 祐史