anchor.diff:  
         integer ma,mb,mp,np, ica,icb,icb_limit,ipa,ipb
         integer Ia,Ja,Ka, Ib,Jb,Kb, Ip,Jp,Kp
  -      integer npp1
   c
   c Compute the nuclear attraction integrals.
   c
  
positive.code:       Subroutine hfnai(E,R0,IJK,Vab,Nints,NPP,La,Lb,Li,Lp,Lp3,canAB)
  c $Id: hfnai.f,v 1.9 2004-07-22 16:31:41 edo Exp $
        
        Implicit none
        
  
negative.code:           Ib = Nxyz(1)
            Jb = Nxyz(2)
            Kb = Nxyz(3)
            call dcopy(NPP,0d0,0,VabP,1)
            do Ip = 0,Ia+Ib
              do Jp = 0,Ja+Jb
                do Kp = 0,Ka+Kb
  