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ページの記述と合っているだろうか?
桂田 祐史