20002
CONTINUE
PHIIJ(I1,I)=PHIIJ(I,I1)
20001
CONTINUE
20000 CONTINUE
DO 21000 J=1, NCAT
DO 21001 I=NC1,NN
SUMCA(J)=SUMCA(J)+XM(I)*(2.*(B(J,I))+ZZ*C(J,I))
21001
CONTINUE
21000 CONTINUE
DO 22000 J=1, NCAT
DO 10 J1=1,NCAT
IF (J.EQ.J1) GO TO 10
DO 22001 I=NC1, NN
SUM=SUM+XM(I)*PSI(J,J1,I)
22001
CONTINUE
SUMCAT(J)=SUMCAT(J)+XM(J1)*(SUM+2.*PHIIJ(J,J1))
SUM=0.
10
CONTINUE
22000 CONTINUE
DO 23000 J=1, NCAT
DO 23001 J1=NC1,NN-1
DO 23002 I=J1+1, NN
SUMAN(J)=SUMAN(J)+XM(J1)*XM(I)*PSI(J1,I,J)
23002
CONTINUE
23001
CONTINUE
23000 CONTINUE
SUM=0.
DO 24000 J=1,NCAT
DO 24001 I=NC1,NN
SUM=SUM+XM(J)*XM(I)*C(J,I)
24001
CONTINUE
24000 CONTINUE
DO 25000 J=1,NCAT
SUMZ(J)=SUM*DABS(Z(J))
25000 CONTINUE
DO 26000 J=1,NCAT
ACT(J+IP)=Z(J)**2*F+SUMCA(J)+SUMCAT(J)+SUMAN(J)+SUMZ(J)
25