      FUNCTION PNDF(X,IND)
C     ---------------
C     A = 1/SQRT(2)
C     C = SQRT(2/PI)
C     ---------------
      DATA A/.707106781186548/
      DATA C/.797884560802865/
C     ---------------
      T = A*X
      IF (IND .NE. 0) GO TO 20
      IF (X .LT. -8.0) GO TO 10
      PNDF = 0.5*ERFC1(0,-T)
      RETURN
   10 PNDF = C/ERFC1(1,-T)
      RETURN
   20 IF (X .GT. 8.0) GO TO 30
      PNDF = 0.5*ERFC1(0,T)
      RETURN
   30 PNDF = C/ERFC1(1,T)
      RETURN
      END
