(defun erf (x) (declare (type double-float x))
 (prog
  ((a (make-array '(4) :element-type 'float))
   (b (make-array '(4) :element-type 'float))
   (p (make-array '(8) :element-type 'float))
   (q (make-array '(8) :element-type 'float))
   (r (make-array '(5) :element-type 'float))
   (s (make-array '(5) :element-type 'float)) (erf 0.0d0) (c 0.0d0) (x2 0.0d0)
   (bot 0.0d0) (top 0.0d0) (t_ 0.0d0) (ax 0.0d0)
  )
  (declare (type (simple-array float (*)) a))
  (declare (type (simple-array float (*)) b))
  (declare (type (simple-array float (*)) p))
  (declare (type (simple-array float (*)) q))
  (declare (type (simple-array float (*)) r))
  (declare (type (simple-array float (*)) s)) (declare (type real erf))
  (declare (type double-float c)) (declare (type double-float x2))
  (declare (type double-float bot)) (declare (type double-float top))
  (declare (type double-float t_)) (declare (type double-float ax))
  (setq c 0.5641896) (replace a '(1.1283792) :end 0)
  (replace a '(0.102201134) :end 0) (replace a '(0.03253241) :end 0)
  (replace a '((- 1.6558183E-4)) :end 0) (replace b '(1.0) :end 0)
  (replace b '(0.42390674) :end 0) (replace b '(0.07013334) :end 0)
  (replace b '(0.0046498897) :end 0) (replace p '(300.45926) :end 0)
  (replace p '(451.91895) :end 0) (replace p '(339.32083) :end 0)
  (replace p '(152.98929) :end 0) (replace p '(43.162228) :end 0)
  (replace p '(7.211758) :end 0) (replace p '(0.5641955) :end 0)
  (replace p '((- 1.3686486E-7)) :end 0) (replace q '(300.45926) :end 0)
  (replace q '(790.9509) :end 0) (replace q '(931.35406) :end 0)
  (replace q '(638.9803) :end 0) (replace q '(277.58545) :end 0)
  (replace q '(77.00015) :end 0) (replace q '(12.782727) :end 0)
  (replace q '(1.0) :end 0) (replace r '(0.2820948) :end 0)
  (replace r '(4.658078) :end 0) (replace r '(21.36882) :end 0)
  (replace r '(26.237015) :end 0) (replace r '(2.1014414) :end 0)
  (replace s '(1.0) :end 0) (replace s '(18.012457) :end 0)
  (replace s '(99.01918) :end 0) (replace s '(187.1148) :end 0)
  (replace s '(94.15378) :end 0) (setf ax (abs x))
  (if (>= ax 0.5) (go label10)) (setf t_ (* x x))
  (setf top
   (+ (* (+ (* (+ (* (fref a 1) t_) (fref a 2)) t_) (fref a 3)) t_) (fref a 4))
  )
  (setf bot
   (+ (* (+ (* (+ (* (fref b 1) t_) (fref b 2)) t_) (fref b 3)) t_) (fref b 4))
  )
  (setf erf (/ (* x top) bot)) (go end_label) label10
  (if (> ax 4.0) (go label20))
  (setf top
   (+
    (*
     (+
      (*
       (+
        (*
         (+
          (*
           (+ (* (+ (* (+ (* (fref p 1) ax) (fref p 2)) ax) (fref p 3)) ax)
            (fref p 4)
           )
           ax
          )
          (fref p 5)
         )
         ax
        )
        (fref p 6)
       )
       ax
      )
      (fref p 7)
     )
     ax
    )
    (fref p 8)
  ))
  (setf bot
   (+
    (*
     (+
      (*
       (+
        (*
         (+
          (*
           (+ (* (+ (* (+ (* (fref q 1) ax) (fref q 2)) ax) (fref q 3)) ax)
            (fref q 4)
           )
           ax
          )
          (fref q 5)
         )
         ax
        )
        (fref q 6)
       )
       ax
      )
      (fref q 7)
     )
     ax
    )
    (fref q 8)
  ))
  (setf erf (+ 0.5 (+ 0.5 (/ (* (* -1 (exp (* (* -1 x) x))) top) bot))))
  (if (< x 0.0) (setf erf (- erf))) (go end_label) label20 (setf erf 1.0)
  (if (>= ax 5.6) (go label21)) (setf x2 (* x x)) (setf t_ (/ 1.0 x2))
  (setf top
   (+
    (*
     (+ (* (+ (* (+ (* (fref r 1) t_) (fref r 2)) t_) (fref r 3)) t_)
      (fref r 4)
     )
     t_
    )
    (fref r 5)
  ))
  (setf bot
   (+
    (*
     (+ (* (+ (* (+ (* (fref s 1) t_) (fref s 2)) t_) (fref s 3)) t_)
      (fref s 4)
     )
     t_
    )
    (fref s 5)
  ))
  (setf erf (/ (+ c (/ (* -1 top) (* x2 bot))) ax))
  (setf erf (+ 0.5 (+ 0.5 (* (* -1 (exp (- x2))) erf)))) label21
  (if (< x 0.0) (setf erf (- erf))) (go end_label) end_label (return erf)
))