3.4 p75 (プログラムリスト4)

FINITE DIFFERENCE METHOD FOR -FNU*D2U/DX2+DU/DX=F

75ページのプログラムで、71ページの演習問題2をやってみよう。

% ./p75
 INPUT : FNU=
1
 INPUT : N=
10
 BACKWARD : METHOD=1,CENTRAL : =2,FORWARD : =3
 INPUT : METHOD=
2
0DIFFUSIVITY=  1.0000E+00   MESH SIZE=  1.0000E-01   METHOD= 2
0** NODAL VALUES OF SOLUTION **
   I     U(I)      I     U(I)      I     U(I)      I     U(I)      I     U(I)
   1  6.1405E-02   2  1.1875E-01   3  1.7160E-01   4  2.1949E-01   5  2.6189E-01
   6  2.9823E-01   7  3.2787E-01   8  3.5011E-01   9  3.6415E-01  10  3.6915E-01
%

これは $N=10$ , $\nu=1$ として中心差分近似を用いて計算したものであるが、 77ページの表にある結果と一致しているだろうか？