subroutine xdasinh (x, result)
      double precision x, result, dasinh
      result = dasinh (x)
      return
      end