anchor.diff:      real(real64) :: a
   
       a = exp(x)
  -    nrm_func = exp(lfwa(a))  * a
  +    nrm_func = exp(log_fwa(a))  * a
       
     end function nrm_func
   
  
positive.code:     real(real64) :: a
  
      a = exp(x)
      m2_func = exp(lfwa(a)) * hdir(a)**2 * a
      
    end function m2_func
  
  
negative.code:     end do
      xs = exp(xs)
      
      fxs = log_fwa(xs)
      ! find amax such that the weight in the integral - fwa(alpha) - is maximal
      i = maxloc(fxs,1,fxs < largest)
      amax = xs(i)
  