double precision function dli (x)
c sept 1983 edition.  w. fullerton, c3, los alamos scientific lab.
      double precision x, sqeps, dei, d1mach, dlog, dsqrt
      external d1mach, dei, dlog, dsqrt
      data sqeps / 0.0d0 /
c
      if (sqeps.eq.0.0d0) sqeps = dsqrt(d1mach(3))
c
      if (x.le.0.d0) call seteru (
     1  41hdli     log integral undefined for x le 0, 41, 2, 2)
      if (x.eq.1.d0) call seteru (
     1  40hdli     log integral undefined for x = 1, 40, 3, 2)
      if (dabs(1.d0-x).lt.sqeps) call seteru (57hdli     answer lt half
     1precision because x too close to 1, 57, 1, 1)
c
      dli = dei (dlog(x))
c
      return
      end