diff options
Diffstat (limited to 'CarpetExtra/IDHydroToy/src/erf.f')
| -rw-r--r-- | CarpetExtra/IDHydroToy/src/erf.f | 13 |
1 files changed, 13 insertions, 0 deletions
diff --git a/CarpetExtra/IDHydroToy/src/erf.f b/CarpetExtra/IDHydroToy/src/erf.f new file mode 100644 index 000000000..30bc84520 --- /dev/null +++ b/CarpetExtra/IDHydroToy/src/erf.f @@ -0,0 +1,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. |