1 files changed, 55 insertions, 0 deletions
diff --git a/Utilib/src/AFERF.f b/Utilib/src/AFERF.f
new file mode 100644
index 0000000..9e92961
--- /dev/null
+++ b/Utilib/src/AFERF.f
@@ -0,0 +1,55 @@
+*DECK AFERF
+      FUNCTION AFERF (XX)
+C
+C-----------------------------------------------------------------------
+C
+C ERROR FUNCTION.
+C
+C-----------------------------------------------------------------------
+C
+      IMPLICIT NONE
+C----
+C FUNCTION
+C----
+      REAL  AFERF
+C----
+C VARIABLES
+C----
+      REAL     A0,A1,A2,A3,A4,A5,A6,B0,B1,B2,B3,
+     &         B4,B5,B6,B7,D1,D2,X,XX,Z,Y
+      INTEGER  IS
+C----
+      DATA A0 / 0.12837912 /, A1 / -0.37612384 /, A2 / 0.11281417 /
+     1 ,   A3 / -0.26782287E-01 /, A4 / 0.50805667E-02 /
+     2 ,   A5 / -0.72514300E-03 /, A6 / 0.58167427E-04 /
+     3 ,   B0 /  0.39141751E-02 /, B1 / -0.17582889E-01 /
+     4 ,   B2 /  0.35873539E-01 /, B3 / -0.42869095E-01 /
+     5 ,   B4 /  0.32161925E-01 /, B5 / -0.11846550E-01 /
+     6 ,   B6 /  0.30705572E-02 /, B7 /  0.59813837E-02 /
+     7 ,   D1,D2 / 1.317, 2.040001 /
+C
+      X=XX
+      IF (X.GE.0.) THEN
+         IS=2
+      ELSE
+         IS=1
+         X=-X
+      ENDIF
+      IF (X.GT.D1) GO TO 8
+      Z=X*X
+      Y=(((((A6*Z+A5)*Z+A4)*Z+A3)*Z+A2)*Z+A1)*Z+A0
+      Y=X*Y+X
+      GO TO 6
+   8  IF (X.GT.D2) GO TO 10
+      X=X-D2
+      Y=((((((B7*X+B6)*X+B5)*X+B4)*X+B3)*X+B2)*X+B1)*X+B0
+      GO TO 7
+  10  Y=1.
+      GO TO 6
+   7  Y=1.-Y
+   6  GO TO (17,18),IS
+  17  AFERF=-Y
+      GO TO 30
+  18  AFERF=Y
+  30  RETURN
+      END