Index: qd-bessel.lisp
===================================================================
--- qd-bessel.lisp	(revision 6ab5226ca3e32b443d87934ec138ff0efc8aaecc)
+++ qd-bessel.lisp	(revision 9fd2ebcbeed3ecae899f732b15fd279f6fb0f14f)
@@ -157,5 +157,6 @@
 	  (bk (bk 0 p)
 	      (bk k p))
-	  (ratio v
+	  ;; Compute g[k](p)/(2*k)!, not r[2*k+1](p)/(2*k)!
+	  (ratio 1
 		 (* ratio (/ (+ v2 (expt (1- (* 2 k)) 2))
 			     (* 2 k (1- (* 2 k))))))
@@ -170,5 +171,5 @@
 	    (format t " term  = ~S~%" term)
 	    (format t " sum   - ~S~%" sum))
-	  (* sum #c(0 2) (/ (exp p) q)))
+	  (* sum 4 (exp p)))
       (when *debug-exparc*
 	(format t "k      = ~D~%" k)
@@ -391,4 +392,7 @@
 (defun bessel-j (v z)
   (let ((vv (ftruncate v)))
+    ;; Clear the caches for now.
+    (an-clrhash)
+    (%big-a-clrhash)
     (cond ((= vv v)
 	   ;; v is an integer