N= 4 M= 4 N= 4 M= 4 1THE FULL MATRIX A OF ORDER 4 IS (PRINTED BY ROWS) 10. 0.I -35. 0.I 50. 0.I -24. 0.I 1. 0.I 0. 0.I 0. 0.I 0. 0.I 0. 0.I 1. 0.I 0. 0.I 0. 0.I 0. 0.I 0. 0.I 1. 0.I 0. 0.I 0THE FULL MATRIX B OF ORDER 4 IS (PRINTED BY ROWS) 1. 0.I 0. 0.I 0. 0.I 0. 0.I 0. 0.I 1. 0.I 0. 0.I 0. 0.I 0. 0.I 0. 0.I 1. 0.I 0. 0.I 0. 0.I 0. 0.I 0. 0.I 1. 0.I 100 CONTINUE L= 1 100 CONTINUE L= 2 100 CONTINUE L= 3 150 CONTINUE, LB= 2 160 CONTINUE, K= 1 150 CONTINUE, LB= 1 160 CONTINUE, K= 2 HESSENBERG RESULTS A 1 1 10.00000 -0.0000000 A 1 2 -35.00000 0.0000000 A 1 3 50.00000 0.0000000 A 1 4 -24.00000 0.0000000 A 2 1 1.000000 0.0000000 A 2 2 -0.0000000 0.0000000 A 2 3 0.0000000 0.0000000 A 2 4 0.0000000 0.0000000 A 3 1 0.0000000 0.0000000 A 3 2 1.000000 0.0000000 A 3 3 0.0000000 0.0000000 A 3 4 0.0000000 0.0000000 A 4 1 0.0000000 0.0000000 A 4 2 0.0000000 0.0000000 A 4 3 1.000000 0.0000000 A 4 4 0.0000000 0.0000000 B 1 1 1.000000 0.0000000 B 1 2 0.0000000 0.0000000 B 1 3 0.0000000 0.0000000 B 1 4 0.0000000 0.0000000 B 2 1 0.0000000 0.0000000 B 2 2 1.000000 0.0000000 B 2 3 0.0000000 0.0000000 B 2 4 0.0000000 0.0000000 B 3 1 0.0000000 0.0000000 B 3 2 0.0000000 0.0000000 B 3 3 1.000000 0.0000000 B 3 4 0.0000000 0.0000000 B 4 1 0.0000000 0.0000000 B 4 2 0.0000000 0.0000000 B 4 3 0.0000000 0.0000000 B 4 4 1.000000 0.0000000 Z 1 1 1.000000 0.0000000 Z 1 2 0.0000000 0.0000000 Z 1 3 0.0000000 0.0000000 Z 1 4 0.0000000 0.0000000 Z 2 1 0.0000000 0.0000000 Z 2 2 1.000000 0.0000000 Z 2 3 0.0000000 0.0000000 Z 2 4 0.0000000 0.0000000 Z 3 1 0.0000000 0.0000000 Z 3 2 0.0000000 0.0000000 Z 3 3 1.000000 0.0000000 Z 3 4 0.0000000 0.0000000 Z 4 1 0.0000000 0.0000000 Z 4 2 0.0000000 0.0000000 Z 4 3 0.0000000 0.0000000 Z 4 4 1.000000 0.0000000 L= 1 ITS 0 GO TO 70, ITS,L,EN = 1 1 4 L= 1 ITS 1 GO TO 70, ITS,L,EN = 2 1 4 L= 1 ITS 2 GO TO 70, ITS,L,EN = 3 1 4 L= 1 ITS 3 GO TO 70, ITS,L,EN = 4 1 4 L= 1 ITS 4 GO TO 70, ITS,L,EN = 5 1 4 L= 1 ITS 5 GO TO 70, ITS,L,EN = 6 1 4 L= 1 ITS 6 GO TO 70, ITS,L,EN = 7 1 4 L= 4 GO TO 60, EN = 3 L= 1 ITS 0 GO TO 70, ITS,L,EN = 1 1 3 L= 1 ITS 1 GO TO 70, ITS,L,EN = 2 1 3 L= 1 ITS 2 GO TO 70, ITS,L,EN = 3 1 3 L= 3 GO TO 60, EN = 2 L= 1 ITS 0 GO TO 70, ITS,L,EN = 1 1 2 L= 1 ITS 1 GO TO 70, ITS,L,EN = 2 1 2 L= 1 ITS 2 GO TO 70, ITS,L,EN = 3 1 2 L= 2 GO TO 60, EN = 1 L= 1 GO TO 60, EN = 0 ALFR(I) ALFI(I) BETA(I) 1 3.999998E+00 -1.110984E-07 1.000000E+00 2 3.000000E+00 -4.952930E-07 9.999999E-01 3 2.000000E+00 1.234871E-07 9.999997E-01 4 9.999998E-01 -1.660908E-09 9.999998E-01 COMPUTED EIGENVALUE AND EIGENVECTOR RESIDUAL 1 3.999997E+00 -1.110984E-07 4.86E-08 8.157446E-01 5.216080E-01 2.039362E-01 1.304020E-01 5.098408E-02 3.260052E-02 1.274603E-02 8.150132E-03 2 3.000000E+00 -4.952930E-07 3.56E-08 9.420009E-01 4.072810E-02 3.140001E-01 1.357604E-02 1.046667E-01 4.525354E-03 3.488888E-02 1.508455E-03 3 2.000001E+00 1.234871E-07 3.08E-08 -6.963045E-01 5.177853E-01 -3.481519E-01 2.588925E-01 -1.740758E-01 1.294462E-01 -8.703786E-02 6.472310E-02 4 9.999999E-01 -1.660908E-09 3.59E-08 -4.998143E-01 -1.364465E-02 -4.998138E-01 -1.364421E-02 -4.998136E-01 -1.364402E-02 -4.998136E-01 -1.364392E-02 AC( 1 , 1 )= (10.0000000000000,0.000000000000000) AC( 1 , 2 )= (-35.0000000000000,0.000000000000000) AC( 1 , 3 )= (50.0000000000000,0.000000000000000) AC( 1 , 4 )= (-24.0000000000000,0.000000000000000) AC( 2 , 1 )= (1.00000000000000,0.000000000000000) AC( 2 , 2 )= (0.000000000000000,0.000000000000000) AC( 2 , 3 )= (0.000000000000000,0.000000000000000) AC( 2 , 4 )= (0.000000000000000,0.000000000000000) AC( 3 , 1 )= (0.000000000000000,0.000000000000000) AC( 3 , 2 )= (1.00000000000000,0.000000000000000) AC( 3 , 3 )= (0.000000000000000,0.000000000000000) AC( 3 , 4 )= (0.000000000000000,0.000000000000000) AC( 4 , 1 )= (0.000000000000000,0.000000000000000) AC( 4 , 2 )= (0.000000000000000,0.000000000000000) AC( 4 , 3 )= (1.00000000000000,0.000000000000000) AC( 4 , 4 )= (0.000000000000000,0.000000000000000) E( 1 )= (3.99999666213989,-1.11098373167806E-7) V( 1 , 1 )= (0.815744638442993,0.521607995033264) V( 1 , 2 )= (0.203936234116554,0.130402013659477) V( 1 , 3 )= (0.0509840771555901,0.0326005183160305) V( 1 , 4 )= (0.0127460258081555,0.00815013237297535) LENGTH = 0.999999972066340 E( 2 )= (3.00000047683716,-4.95293022595433E-7) V( 2 , 1 )= (0.942000925540924,0.0407280959188938) V( 2 , 2 )= (0.314000129699707,0.0135760353878140) V( 2 , 3 )= (0.104666665196419,0.00452535413205624) V( 2 , 4 )= (0.0348888821899891,0.00150845525786281) LENGTH = 1.00000000543826 E( 3 )= (2.00000095367432,1.23487126302280E-7) V( 3 , 1 )= (-0.696304500102997,0.517785310745239) V( 3 , 2 )= (-0.348151862621307,0.258892506361008) V( 3 , 3 )= (-0.174075797200203,0.129446208477020) V( 3 , 4 )= (-0.0870378613471985,0.0647230967879295) LENGTH = 1.00000000340474 E( 4 )= (0.999999940395355,-1.66090818876086E-9) V( 4 , 1 )= (-0.499814331531525,-0.0136446505784988) V( 4 , 2 )= (-0.499813824892044,-0.0136442072689533) V( 4 , 3 )= (-0.499813556671143,-0.0136440200731158) V( 4 , 4 )= (-0.499813556671143,-0.0136439222842455) LENGTH = 1.00000003260125 ERR AV-EV at( 1 )= (0.00000201566596087490,0.00000253328620389915) ERR AV-EV at( 2 )= (0.00000202039842210208,8.26816510161495E-7) ERR AV-EV at( 3 )= (-0.00000174936697071343,5.01947825818809E-7) ERR AV-EV at( 4 )= (0.00000152005347823944,-0.00000220449488633051) V( 1 ) dot V( 2 )= 0.289297672516016 V( 1 ) dot V( 3 )= 0.148902583990480 V( 1 ) dot V( 4 )= 0.0979232852083672 V( 2 ) dot V( 3 )= 0.0979232026147917 V( 2 ) dot V( 4 )= 0.0779240917101936 V( 3 ) dot V( 4 )= 0.0695257915277268 ILAMBDA= 1 (3.99999666213989,-1.11098373167806E-7) SINGULARITY= 1 ILAMBDA= 2 (3.00000047683716,-4.95293022595433E-7) SINGULARITY= 1 ILAMBDA= 3 (2.00000095367432,1.23487126302280E-7) SINGULARITY= 1 ILAMBDA= 4 (0.999999940395355,-1.66090818876086E-9) SINGULARITY= 1 TOTAL EIGENVECTOR ERROR = 0.0000112849599758654 N= 0 M= 0 END OF DATA FOR SUBROUTINE RMATIN