double precision function dei (x)
  c july 1977 edition.  w. fullerton, c3, los alamos scientific lab.
        double precision x, de1
        external de1
  c
        dei = -de1(-x)
  c
        return
        end