1 files changed, 13 insertions, 0 deletions
diff --git a/cl/src/approx/maceps.lisp b/cl/src/approx/maceps.lisp
new file mode 100644
index 0000000..cd24482
--- /dev/null
+++ b/cl/src/approx/maceps.lisp
@@ -0,0 +1,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))))