cl/src/approx/maceps.lisp


1
2
3
4
5
6
7
8
9
10
11
12
13

(in-package :lizfcm.approx)

(defun compute-maceps (val f)
  (let* ((h val)
         (a val)
         (err val))
    (loop while (> err 0)
          do
          (progn
            (setf h (/ h 2)
                  err (abs (- (funcall f (+ a h))
                              (funcall f a)))))
          collect (list a h err))))