anchor.diff:  
   END DO
   
  -NATOM = N - 1
  -WRITE(6,*) NATOM
  +NSPHERE = N - 1
  +WRITE(6,"(/,A,I5,/)") " TOTAL SPHERES = ", NSPHERE
   
   EPS = 10.0**(-NEPS)
   S1 = 0.0
  
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:         WRITE(6,"(I5,E11.4)") IT, EIG(K)
  
          IF (ABS((EIG(K) - EA) / EIG(K)) .GE. EPS) THEN
              EA = EIG(K)
          ELSE
              EXIT
          END IF
  
          IF (IT .EQ. NMAX) THEN
              WRITE(9,"(A)") " MAXIMUM OF ITERATIONS REACHED"
              EXIT
          END IF
  
          IT = IT + 1
  
      END DO
  
  
  END DO
  
  