[交流] 【求助】关于调试程序的

C234567
      SUBROUTINE S72(K,M,KK,NQ,I1,I2,C)
      COMMON /KP/KP(0:49,15),/KMN/KMN(31,0:11)
      DIMENSION KK(M),KI(25),K(9)
      K1=MAX(0,K(9)-K(5),K(9)-K(1))
      K2=MIN(K(2),K(4),K(9))
      I1=0
      DO 1 I=K1,K2
      I1=I1+(-1)**(K(1)-K(5)+I)*KMN(K(2)-I+1,I-K1)*
     1 KMN(K(4)-I+1,I-K1)*KMN(K(9)-I+1,I-K1)*
     1 KMN(I+1,K2-I)*KMN
     1 (K(5)-K(9)+I+1,K2-I)*KMN(K(1)-K(9)+I+1,K2-I)     
1     CONTINUE
      IF(I1) 2,3,2
3     C=0
      I1=0
      I2=1
      RETURN
2     CALL PRIMN(I1,M,KK,NQ)
      KI(1:M)=KK(1:M)-KP(K(2)-K1,1:M)-KP(K(4)-K1,
     1 1:M)-KP(K2,1:M)-KP(K(9)-K1,1:M)
     1 -KP(K(5)-K(9)+K2,1:M)-KP(K(1)-K(9)+K2,1:M)
      KK(1:M)=2*KI(1:M)-KP(K(1)+K(2)+K(3)+1,1:M)
      DO 4 I=1,9
      KK(1:M)=KK(1:M)+KP(K(I),1:M)
4      CONTINUE
      CALL I1I26J(M,KK,NQ,I1,I2)
      C=SQRT(REAL(I1)/REAL(I2))*ISIGN(1,NQ)
      RETURN
      END