subroutine xderf (x, result)
      double precision x, result, derf
      result = derf (x)
      return
      end