function chi (x)
  c december 1980 edition.  w. fullerton, bell labs.
  c
  c evaluate the hyperbolic cosine integral.  when x is negative, the
  c principal value is used.
  c
        external e1, ei
  c
        chi = 0.5 * (ei(x) - e1(x))
  c
        return
        end