home *** CD-ROM | disk | FTP | other *** search
/ Programmer 7500 / MAX_PROGRAMMERS.iso / INFO / FORTRAN / D2SOURCE.ZIP / ERFF.FOR < prev    next >
Encoding:
Text File  |  1985-12-03  |  495 b   |  19 lines
  1. C  ERROR FUNCTION (HASTINGS APPROXIMATION)
  2.       FUNCTION ERFF (A)
  3.       DATA A1/0.070523078/,A2/0.042282012/,A3/0.0092705272/
  4.       DATA A4/0.0001520143/,A5/0.0002765672/,A6/.0000430638/
  5.       DATA C/1./
  6.       C=SIGN(C,A)
  7.       Z=ABS(A)
  8.       IF (Z.LT..000001) GO TO 1
  9.       IF (Z.GT.4.) GO TO 2
  10.       ERR=((((((A6*Z+A5)*Z+A4)*Z+A3)*Z+A2)*Z+A1)*Z+1.)**16
  11.       ERR=ABS(1.-(1./ERR))
  12.       ERFF=ERR*C
  13.       RETURN
  14.  1    ERFF=0.
  15.       RETURN
  16.  2    ERFF=C
  17.       RETURN
  18.       END
  19.