function ei (x)
c april 1977 version.  w. fullerton, c3, los alamos scientific lab.
c
      external e1
c
      ei = -e1(-x)
c
      return
      end