CarpetExtra/IDHydroToy/src/erf.f


1
2
3
4
5
6
7
8
9
10
11
12
13

      FUNCTION erf(x)
      implicit none
      REAL*8 erf,x
CU    USES gammp
      REAL*8 gammp
      if(x.lt.0)then
        erf=-gammp(.5d0,x**2)
      else
        erf=gammp(.5d0,x**2)
      endif
      return
      END
C  (C) Copr. 1986-92 Numerical Recipes Software t4.