anchor.diff:              B = 1.0 - FACT
   
               ! Calculation of the diagonal part of the Oseen tensor only
  -            DO 260  I=1,IP
  +            DO I=1,IP
                   V=A+(DIS(I)*DIS(I)*B/D2)
                   T(K,L,I)=V*6.0/(8.0*DISTAN)
                   T(L,K,I)=T(K,L,I)
               END DO
  +        END IF
   
       END DO
   END DO
  
positive.code:     E3 = E2 * E(K)
      S4 = S4 + E3
  
      DO L= 1, NATOM
          EL2 = E(L)**2
          IF (K .LT. L) THEN
              D2 = 0.0
  
negative.code: END DO
  
  NATOM = N - 1
  WRITE(6,*) NATOM
  
  EPS = 10.0**(-NEPS)
  S1 = 0.0
  