subroutine xdacosh (x, result)
      double precision x, result, dacosh
      result = dacosh (x)
      return
      end