subroutine xdgami (a, x, result)
      double precision a, x, result, dgami
      result = dgami (a, x)
      return
      end