a variacional del estado de transición a la - Páxinas persoais - USC ...

C. Programas 143CDO 30 J=IP1,N30 SUM=SUM+A(I+(J*(J-1))/2)*W(2,J)40 W(1,I)=SUM/H50 SUMM=W(1,I)*W(2,I)+SUMMU=SUMM*0.5D0/HDO 60 J=KP1,NW(1,J)=W(2,J)*U-W(1,J)DO 60 I=KP1,J60 A(I+(J*(J-1))/2)=W(1,I)*W(2,J)+W(1,J)*W(2,I)+A(I+(J*(J-1))/2)70 A((K*(K+1))/2)=H80 W(2,NM1)=A((NM1*(NM1+1))/2)W(2,N)=A((N*(N+1))/2)W(1,NM1)=A(NM1+(N*(N-1))/2)W(1,N)=0.D0GERSCH=ABS(W(2,1))+ABS(W(1,1))DO 90 I=1,NM190 GERSCH=MAX(ABS(W(2,I+1))+ABS(W(1,I))+ABS(W(1,I+1)),GERSCH)DEL=EPS*GERSCHDO 100 I=1,NW(3,I)=W(1,I)E(I)=W(2,I)100 V(I,M)=E(I)IF(DEL.EQ.ZERO) GOTO 210QR-METHOD WITH ORIGIN SHIFTK=N110 L=K120 IF(ABS(W(3,L-1)).LT.DEL) GOTO 130L=L-1IF(L.GT.1) GOTO 120130 IF(L.EQ.K) GOTO 160WW=(E(K-1)+E(K))*0.5D0R=E(K)-WWZ=SIGN(SQRT(W(3,K-1)**2+R*R),R)+WWEE=E(L)-ZE(L)=EEFF=W(3,L)R=SQRT(EE*EE+FF*FF)J=LGOTO 150140 R=SQRT(E(J)**2+W(3,J)**2)W(3,J-1)=S*REE=E(J)*CiFF=W(3,J)*Ci150 Ci=E(J)/RS=W(3,J)/RWW=E(J+1)-ZE(J)=(FF*Ci+WW*S)*S+EE+ZE(J+1)=Ci*WW-S*FFJ=J+1IF(J.LT.K) GOTO 140W(3,K-1)=E(K)*SE(K)=E(K)*Ci+ZGOTO 110160 K=K-1IF(K.GT.1) GOTO 110