anchor.diff:  
          E_now(iE) = 0.0d0
          do idist = 1, ndist(iE)
  -    
  +
             call eval_fitfunct(dist(idist,iE), ftype, np, param, E, &
                                dEdr, d2Edr2, d3Edr3)
   
  
positive.code:        d3y = d3fx - d3tx
  
      end if
      
    end subroutine eval_polyexp
  
    !--------------------------------------------------------------------!
  
negative.code:                 end do
                  if (newdist) then
                     if (ndist(nE) >= ndist_max) call more_memory()
                     ndist(nE) = ndist(nE) + 1
                     p_dist(ndist(nE), nE)  =  r
                     p_count(ndist(nE), nE) =  1
                 end if
  