■計算の実行■
上のプログラムを起動させ、k=1, Δt=0.1, NSTEP=10 を入力し、α＝1, 0, 0.5 で実行すると、x=L/2 での温度T(x=π/2,t)は、下図の様な結果になります。NSTEPは、時間ステップ数のことで、t=NSTEP*ΔtでのT(x,t)の計算を終了すると、プログラムの実行は終了されます。

Crank-Nicolson metod は、厳密解とほぼ同じ結果を出していますね。