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