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