c **********************************************************************
      function f0(x)
c **********************************************************************
c
c     function f0(x)=1.0d0/sqrt(x)
c
c     0 < x <1.0d-32 it returns 1.0d0/sqrt(1.0d-32)
c     x < 0          it returns 0
c
c **********************************************************************
      implicit real*8 (a-h,o-z)
c
      eps=1.0d-32
      if(x.lt.0) then
         f0=0.0d0
      else if(x.lt.eps) then
         f0=1.0d0/sqrt(eps)
      else
         f0=1.0d0/sqrt(x)
      endif
c
      end