3.3 p73 (プログラムリスト3)

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

% ./p73
 INPUT : N=
10
 INPUT : M=
1000
 INPUT : EPS=
1E-5
 JACOBI : METHOD=0, SOR : =1
 INPUT : METHOD=
0
 (NITER,ERROR)=           1   1.00000000    
 (NITER,ERROR)=           2  0.500000000    
 (NITER,ERROR)=           3  0.333333373    
 (NITER,ERROR)=           4  0.250000030    
 (NITER,ERROR)=           5  0.200000092    
(中略)
 (NITER,ERROR)=         580   1.06710586E-05
 (NITER,ERROR)=         581   1.04721094E-05
 (NITER,ERROR)=         582   1.04720166E-05
 (NITER,ERROR)=         583   1.02067979E-05
 (NITER,ERROR)=         584   1.02067088E-05
 (NITER,ERROR)=         585   9.94150014E-06
0NO. OF ITERATIONS= 585   EPS=  1.0000E-05
0MESH SIZE=  1.0000E-01
0** NODAL VALUES OF SOLUTION **
   I     U(I)      I     U(I)      I     U(I)      I     U(I)      I     U(I)
   1  8.9947E-02   2  1.6990E-01   3  2.3985E-01   4  2.9980E-01   5  3.4976E-01
   6  3.8973E-01   7  4.1970E-01   8  4.3968E-01   9  4.4967E-01  10  4.4966E-01
%

テキストには589回で停止したと書いてあったが、585回で停止し、 得られた解もわずかに異なっている。 ちなみにCに書き換えたプログラムは 587 回で停止した。