590
16 Computer Programs in Vibration Analysis
Table 16.6 (continued)
DO 166 IT=1,100
IF(MARK)168,108,168
108 DO 114 I=1,NM1
IP1=I+1
DO 114 J=IP1,N
ENO=ABS(A(I,J)+A(J,I))
IF(ENO-EPS)110,110,116
110 ENO=ABS(A(I,J)-A(J,I))
IF(ENO-EPS)114,114,112
112 ENO=ABS(A(I,I)-A(J,J))
IF(ENO-EPS*10.0)114,114,116
114 CONTINUE
GO TO 168
116 MARK=1
DO 164 K=1,NM1
KP1=K+1
DO 164 M=KP1,N
H=0.0
G=0.0
HJ=0.0
YH=0.0
DO 122 I=1,N
TE=A(I,K)*A(I,K)
TEE=A(I,M)*A(I,M)
YH=YH+TE-TEE
IF(I-K) 118,122,118
118 IF(I-M)120,122,120
120 AIK=A(I,K)
AKI=A(K,I)
AIM=A(I,M)
AMI=A(M,I)
H=H+AKI*AMI-AIK*AIM
TEP=AMI*AMI+TE
TEM=AKI*AKI+TEE
G=G+TEP+TEM
HJ=HJ-TEP+TEM
122 CONTINUE
H=2.0*H
D=A(K,K)-A(M,M)
TEP=A(K,M)
TEM=A(M,K)
C=TEP+TEM
E=TEP-TEM
IF(ABS(C)-EP)124,124,126
124 CC=1.0
SS=0.0
GO TO 134
126 BY=D/C
IF(BY) 128,130,130
128 SIG=-1.0
GO TO 132
130 SIG=1.0
132 COT=BY+(SIG*SQRT(BY*BY+1.0))
SS=SIG/SQRT(COT*COT+1.0)
(continued)