N= 4 M= 4 N= 4 M= 4 1THE FULL MATRIX A OF ORDER 4 IS (PRINTED BY ROWS) 1. 4.I 2. 4.I 3. 4.I 4. 4.I 2. 4.I 2. 3.I 3. 3.I 4. 3.I 3. 4.I 3. 3.I 3. 2.I 4. 2.I 4. 4.I 4. 3.I 4. 2.I 4. 1.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= 1 150 CONTINUE, LB= 2 160 CONTINUE, K= 1 150 CONTINUE, LB= 1 160 CONTINUE, K= 2 HESSENBERG RESULTS A 1 1 1.000000 4.000000 A 1 2 4.270862 7.332989 A 1 3 1.798025 -1.143815 A 1 4 0.2880772 -0.6023436 A 2 1 8.774965 0.0000000 A 2 2 12.13101 -2.396817 A 2 3 -1.860792 -3.028406 A 2 4 -0.6744961 -0.6177906 A 3 1 0.0000000 0.0000000 A 3 2 1.100300 0.0000000 A 3 3 1.007121 -0.05441273 A 3 4 0.1483012 -0.1753773 A 4 1 0.0000000 0.0000000 A 4 2 0.0000000 0.0000000 A 4 3 0.2057403 0.0000000 A 4 4 0.5089182 -0.2034380 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 0.7071067 -0.7071067 B 2 3 -1.827395E-8 4.319577E-8 B 2 4 -1.029549E-8 1.792370E-8 B 3 1 0.0000000 0.0000000 B 3 2 0.0000000 0.0000000 B 3 3 -0.8320503 0.5547001 B 3 4 -1.353349E-8 1.606769E-8 B 4 1 0.0000000 0.0000000 B 4 2 0.0000000 0.0000000 B 4 3 0.0000000 0.0000000 B 4 4 -0.4350202 0.9004208 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 0.4834937 0.1611646 Z 2 3 -0.8267329 -0.1667183 Z 2 4 0.1702277 0.0000000 Z 3 1 0.0000000 0.0000000 Z 3 2 0.5640761 0.08058229 Z 3 3 0.1825459 0.006331075 Z 3 4 -0.7856658 0.1571331 Z 4 1 0.0000000 0.0000000 Z 4 2 0.6446584 0.0000000 Z 4 3 0.5012102 -0.06436594 Z 4 4 0.5401452 -0.1931428 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= 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= 2 GO TO 60, EN = 1 L= 1 GO TO 60, EN = 0 ALFR(I) ALFI(I) BETA(I) 1 1.215895E+01 1.215895E+01 1.000000E+00 2 -1.340416E+00 -1.340416E+00 1.000000E+00 3 -5.280409E-01 -5.280409E-01 1.000000E+00 4 -2.904939E-01 -2.904939E-01 1.000000E+00 COMPUTED EIGENVALUE AND EIGENVECTOR RESIDUAL 1 1.215895E+01 1.215895E+01 7.76E-08 4.269603E-01 3.033994E-01 4.202161E-01 2.215222E-01 4.480319E-01 1.578640E-01 5.126958E-01 1.071891E-01 2 -1.340416E+00 -1.340416E+00 4.08E-08 5.542877E-02 -6.332651E-01 -8.853878E-02 -2.967555E-01 -1.664530E-01 2.611445E-01 -1.201872E-01 6.242211E-01 3 -5.280408E-01 -5.280408E-01 2.16E-08 3.141842E-02 4.708746E-01 8.321620E-02 -5.199665E-01 -2.258013E-02 -5.260986E-01 -8.561458E-02 4.640910E-01 4 -2.904938E-01 -2.904938E-01 1.80E-08 2.198985E-01 -1.844457E-01 -5.241373E-01 3.780266E-01 5.361247E-01 -3.608252E-01 -2.491761E-01 1.424325E-01 AC( 1 , 1 )= (1.00000000000000,4.00000000000000) AC( 1 , 2 )= (2.00000000000000,4.00000000000000) AC( 1 , 3 )= (3.00000000000000,4.00000000000000) AC( 1 , 4 )= (4.00000000000000,4.00000000000000) AC( 2 , 1 )= (2.00000000000000,4.00000000000000) AC( 2 , 2 )= (2.00000000000000,3.00000000000000) AC( 2 , 3 )= (3.00000000000000,3.00000000000000) AC( 2 , 4 )= (4.00000000000000,3.00000000000000) AC( 3 , 1 )= (3.00000000000000,4.00000000000000) AC( 3 , 2 )= (3.00000000000000,3.00000000000000) AC( 3 , 3 )= (3.00000000000000,2.00000000000000) AC( 3 , 4 )= (4.00000000000000,2.00000000000000) AC( 4 , 1 )= (4.00000000000000,4.00000000000000) AC( 4 , 2 )= (4.00000000000000,3.00000000000000) AC( 4 , 3 )= (4.00000000000000,2.00000000000000) AC( 4 , 4 )= (4.00000000000000,1.00000000000000) E( 1 )= (12.1589498519897,12.1589469909668) V( 1 , 1 )= (0.426960319280624,0.303399384021759) V( 1 , 2 )= (0.420216053724289,0.221522241830826) V( 1 , 3 )= (0.448031872510910,0.157863974571228) V( 1 , 4 )= (0.512695789337158,0.107189096510410) LENGTH = 1.00000000198310 E( 2 )= (-1.34041595458984,-1.34041607379913) V( 2 , 1 )= (0.0554287731647491,-0.633265137672424) V( 2 , 2 )= (-0.0885387808084488,-0.296755522489548) V( 2 , 3 )= (-0.166453033685684,0.261144459247589) V( 2 , 4 )= (-0.120187208056450,0.624221086502075) LENGTH = 1.00000000507702 E( 3 )= (-0.528040766716003,-0.528040766716003) V( 3 , 1 )= (0.0314184240996838,0.470874577760696) V( 3 , 2 )= (0.0832161977887154,-0.519966542720795) V( 3 , 3 )= (-0.0225801281630993,-0.526098608970642) V( 3 , 4 )= (-0.0856145843863487,0.464090973138809) LENGTH = 1.00000001171746 E( 4 )= (-0.290493845939636,-0.290493816137314) V( 4 , 1 )= (0.219898506999016,-0.184445679187775) V( 4 , 2 )= (-0.524137318134308,0.378026604652405) V( 4 , 3 )= (0.536124706268311,-0.360825181007385) V( 4 , 4 )= (-0.249176144599915,0.142432540655136) LENGTH = 0.999999947816466 ERR AV-EV at( 1 )= (-0.00000125048575227993,0.00000596295013366444) ERR AV-EV at( 2 )= (-0.00000144830911708027,-8.02528624888055E-7) ERR AV-EV at( 3 )= (-1.68092070040871E-7,9.80752391299688E-7) ERR AV-EV at( 4 )= (-0.00000103186709399239,9.48437297587645E-8) V( 1 ) dot V( 2 )= 0.000608932741155220 V( 1 ) dot V( 3 )= 0.000149367561244979 V( 1 ) dot V( 4 )= -0.00394079529726255 V( 2 ) dot V( 3 )= 0.00368805345568815 V( 2 ) dot V( 4 )= 0.000149362855503335 V( 3 ) dot V( 4 )= 0.000608898273934629 ILAMBDA= 1 (12.1589498519897,12.1589469909668) SINGULARITY= 1 ILAMBDA= 2 (-1.34041595458984,-1.34041607379913) SINGULARITY= 1 ILAMBDA= 3 (-0.528040766716003,-0.528040766716003) SINGULARITY= 1 ILAMBDA= 4 (-0.290493845939636,-0.290493816137314) SINGULARITY= 1 TOTAL EIGENVECTOR ERROR = 0.0000117882489800697 N= 0 M= 0 END OF DATA FOR SUBROUTINE RMATIN