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