5.2 p104 (プログラムリスト2)

104ページのプログラムで、103-106ページの計算実習2をやってみよう。

% ./p104
 INPUT : N=
4
 INPUT : SHIFT=
0
 INPUT : EPS=
1E-6
 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   1.10236216       1.00000000
 (NITER,EVAL,ERROR)=           2   1.00199914      0.100162774
 (NITER,EVAL,ERROR)=           3   1.00005829       1.94073329E-03
 (NITER,EVAL,ERROR)=           4   1.00000215       5.61474553E-05
 (NITER,EVAL,ERROR)=           5   1.00000012       2.02655769E-06
 (NITER,EVAL,ERROR)=           6  0.999999881       2.38418608E-07
0NO. OF ITERATIONS=  6   EPS=  1.0000E-06   SHIFT=  0.0000E+00   NPIVOT=  0
0EIGENVALUE=  1.0000E+00
0** EIGENVECTOR **
   I     X(I)      I     X(I)      I     X(I)      I     X(I)      I     X(I)
   1  7.0712E-01   2 -7.0710E-01   3 -1.7818E-05   4 -1.7818E-05

ここでは 106 ページの1の場合を計算している。 106ページの記述と合っているだろうか?

桂田 祐史
2018-06-08