KF()BTC U/W-FOCAL:EIGVAL15.10.78?MB4BJC THIS PROGRAM FINDS THE EIGENVALUES OF THE MATRIX A(I,J)?M?CSBTC USING THE QR ALGORITHM.  THE MATRIX IS NOT REQUIRED TO?MC2B^C BE SYMMETRIC - COMPLEX EIGENVALUES WILL BE COMPUTED IF?MDQB(C NECESSARY.  THE ORIGINAL MA
RIX IS SAVED IN C(I,J) FOR?MD B2C USE LATER BY 	EIGVEC	.?MDB<?MD7CFF I=1,!; F J=1,!; S C(I,J)=A(I,J)?M?EIDJS K=A(!,!), M=1, N=!, %=1E-10?M?E^DTD .3; F L=M,N-1;D 3;N; F L=M,N-1;D 4?ME2D^S K=-K; F L=M,N; S A(L,L)=A(L,L)+K?MFJD(Y -L; IF (L-2).5; IF (%-FABS(A
L,L-1))).4?M?F(D2S K=A(N,N),M=L; I (L+1-N).2,.7; S X(N)=K;Z Y(N); G .9?M?GFD<S X(L)=X(N)=I,Y(N)=-Y(L)=FSQT(-J); T "COMPLEX"!; G .9?M?G$EFS I=(A(L,L)+K)/2; IF (J=I*I+A(L,N)*A(N,L)-A(L,L)*K).6?M?HEPS X(L)=I+J=FSQT(J)*FSGN(I+%), X(N)=I-J;Z Y(L) Y(N)?MH\EZI (
-N=L=L-1).4,.5;RETURN;EIGENVALUES IN X(I),Y(I)?MH,FJZ $;F I=L,N;S $=$+A(I,L)^2?MIBFTS B(L,L)=A(L,L)+$=FSGN(A(L,L))*FSQT($)?MIZF^S Q(L)=B(L,L)*$;F I=L+1,N;S B(I,L)=A(I,L)?M?I5F(F J=L,N;Z P(J);F I=L,N;S P(J)=P(J)+B(I,L)*A(I,J)?MJPF2F J=L,N;F I=L,N;S A(I,J)=
(I,J)-B(I,L)*P(J)/Q(L)?MJ+HJF I=M,N;Z P(I);F J=L,N;S P(I)=P(I)+A(I,J)*B(J,L)?MHTF I=M,N;F J=L,N;S A(I,J)=A(I,J)-P(I)*B(J,L)/Q(L)?M?ML D EIGVAL;L S EIGVAL;E?M? TURN NULLS INTO ZEROS?MLFJFOR I=4,-1,1;SET C(I)=-	?; CLEAR OUT PREVIOUS CHARACTERS?M?L
FTIF (255-C=FIN()),.3;IF (141-C),1.1;SET C(I=I+1)=-C;GO .4?ML?F^SET C(I)=-	?,I=I-FSGN(I);IF (C).2;GO .2+FOUT(92);RUBOUTS?M?F(IF (C(I)+	:),.*C(I)=-	?;GO .2; REPLACE A 	:	 WITH A NULL?M?L?L D HEADER;L S HEADER;E?M?()-	0)*2^(5-J)?M?ML C LTRHED?M?MO C?M?

Note: Lines longer than 256 characters were wrapped