subroutine xdgamit (a, x, result)
      double precision a, x, result, dgamit
      result = dgamit (a, x)
      return
      end