Listing of file='12' on disk='vmedia/glbr_22.wvd.zip'
# Sector 370, program filename = '12'
0010 DIM A(5,6),S(5,5)
0020 R=.00001
0030 PRINT "INPUT N"
: INPUT N
0040 PRINT "INPUT MATRIX"
0050 FOR I=1 TO N
0060 INPUT A(I,1),A(I,2),A(I,3),A(I,4),A(I,5),A(I,6),A(I,7),A(I,8),
A(I,9),A(I,10)
0075 NEXT I
0080 FOR I=1 TO N
: S(I,I)=1
: NEXT I
0090 FOR I=2 TO N
: FOR J=1 TO I-1
0100 I1=I1+2*A(I,J)^2
0110 NEXT J
: NEXT I
0130 N1=SQR(I1)
: N2=(R/N)*N1
: T=N1
0140 T=T/N
0150 FOR Q=2 TO N
: FOR P=1 TO Q-1
0160 IF ABS(A(P,Q))<=T THEN 320
0170 I2=1
0180 V1=A(P,P)
: V2=A(P,Q)
: V3=A(Q,Q)
0190 M1=(V1-V3)*.5
: IF M1<>0 THEN 210
0200 W=-1
: GOTO 220
0210 W=-SGN(M1)*V2/SQR(V2^2+M1^2)
0220 T1=W/SQR(2*(1+SQR(1-W/2)))
: T2=T1^2
0230 C1=SQR(1-T2)
: C2=C1^2
: T3=T1*C1
0240 FOR I=1 TO N
0250 I1=A(I,P)*C1-A(I,Q)*T1
: A(I,Q)=A(I,P)*T1+A(I,Q)*C1
0260 A(I,P)=I1
: I1=S(I,P)*C1-S(I,Q)*T1
0270 S(I,Q)=S(I,P)*T1+S(I,Q)*C1
: S(I,P)=I1
0280 NEXT I
0290 FOR I=1 TO N
: A(P,I)=A(I,P)
: A(Q,I)=A(I,Q)
: NEXT I
0300 A(P,P)=V1*C2+V3*T2-2*V2*T3
: A(Q,Q)=V1*T2+V3*C2+2*V2*T3
0310 A(P,Q)=(V1-V3)*T3+V2*(C2-T2)
: A(Q,P)=A(P,Q)
0320 NEXT P
0330 NEXT Q
0340 IF I2<>1 THEN 360
0350 I2=0
: GOTO 150
0360 IF T>N2 THEN 140
0370 PRINT " EIGENVALUE"," EIGENVECTOR MATRIX"
0380 FOR I=1 TO N
: PRINT A(I,I);
0390 FOR J=1 TO N
: PRINT TAB(16);S(J,I),
: NEXT J
0400 PRINT
: NEXT I
0999 END