subroutine xdatanh (x, result)
      double precision x, result, datanh
      result = datanh (x)
      return
      end