DO 3I-K.N P=ARK, RKJ AR(S.J) R(S.J》 AI(K.I)AI(S.] AI(IS,D=P 30 CONTINUE D0 35J-1.M PR(K BRCK.D-BRUIS.D RQS.J)-F P-BIK,J S》一上 35 CONTINUE D 50I-1.N P=ARQ,K) AR1,K)=A.SK》 AROIS(K=P P=AICK) AI(I.K)-AI(I.JS(K)) K-P 50 DO60J-K-IN P=AR(K,J】ARK,K Q-A(KJ)AI(K,K) S-(AR(K,K)-AI(K,K))+(AR(K.J)+AICK.J)) AR(K.J)-(PQ)/ A0K.J)=5-P-Q)石 CONTINUE D065J=1M Q=-BK,J·0K,K) S=(ARK,K》一A1(K,K)e(BR(X,J》风K,I) BR(K=(P-0)/D B团(K,J》=sP-Q/D CONTINUE D090I=1,N 压L.NE.)THE 080J=K+1. P-ARKL.K)+AR(K.D)