subroutine xdgamma (x, result)
      double precision x, result, dgamma
      result = dgamma (x)
      return
      end