diff --git a/qd-bessel.lisp b/qd-bessel.lisp
index 37a4c8e1548e4e954c3ebd804c58478832bdb0c5..bf3ca04067d28dbfb69491a0531b9a45e4f67c25 100644
--- 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