      subroutine xgamma (x, result)
      external gamma
      real x, result, gamma
      result = gamma (x)
      return
      end
