subroutine xderfc (x, result)
      double precision x, result, derfc
      result = derfc (x)
      return
      end