sum-an was summing one too many terms.

diff --git a/qd-bessel.lisp b/qd-bessel.lisp
index 37a4c8e..bf3ca04 100644 (file)
--- a/qd-bessel.lisp
+++ b/qd-bessel.lisp
@@ -350,7 +350,7 @@
                   (an n k v)))
          (sum term (+ sum term)))
         ((or (<= (abs term) (* eps (abs sum)))
-             (> k big-n))
+             (>= k big-n))
          sum))))
 
 ;; SUM-AB computes the series