4.1 p88 (プログラムリスト1)

88ページのプログラムで、89-90 ページの計算実習1をやってみよう。

% ./p88
 INPUT : N=
4
 INPUT : EPS=
1E-5
 INPUT : M=
100
 INPUT : A(  1,  1)=
5
 INPUT : A(  1,  2)=
4
 INPUT : A(  1,  3)=
1
 INPUT : A(  1,  4)=
1
(中略)
 INPUT : A(  4,  1)=
1
 INPUT : A(  4,  2)=
1
 INPUT : A(  4,  3)=
2
 INPUT : A(  4,  4)=
4
 INPUT : X(  1)=
1
 INPUT : X(  2)=
0
 INPUT : X(  3)=
0
 INPUT : X(  4)=
0
 (NITER,EVAL,ERROR=           1   9.60465145       1.00000000
 (NITER,EVAL,ERROR=           2   9.92197895       3.19822803E-02
 (NITER,EVAL,ERROR=           3   9.98053360       5.86688565E-03
 (NITER,EVAL,ERROR=           4   9.99512196       1.45954755E-03
 (NITER,EVAL,ERROR=           5   9.99878025       3.65874090E-04
 (NITER,EVAL,ERROR=           6   9.99969292       9.12694377E-05
 (NITER,EVAL,ERROR=           7   9.99992371       2.30790938E-05
 (NITER,EVAL,ERROR=           8   9.99998188       5.81742370E-06
0NO. OF ITERATIONS=  8   EPS=  1.0000E-05
0EIGENVALUE=  1.0000E+01
0** EIGENVECTOR **
   I     X(I)      I     X(I)      I     X(I)      I     X(I)      I     X(I)
   1  6.3307E-01   2  6.3307E-01   3  3.1499E-01   4  3.1499E-01
[chronos:q1/progs/chapter4] mk%

この結果が 90 ページの記述と合っているだろうか?

桂田 祐史
2018-06-08