Changeset 95aa580ce0dd62113535c30c4c0c82f8656c2d86 for qd-bessel.lisp

Timestamp:

04/08/12 09:37:22 (14 months ago)

Author:

Raymond Toy <toy.raymond@…>

Children:

c7fc989dad090ada2c046b28bf43406810474a77

Parents:

dbc1e376add1d492dc35c37b3098c2ffb796d6f1

git-committer:

Raymond Toy <toy.raymond@…> (04/08/12 09:37:22)

Message:

T is not a variable. Use a different name.

Files:

• qd-bessel.lisp (modified) (1 diff)

qd-bessel.lisp

@@ -321 +321 @@
 ;; So incomplete_gamma_tail(1-v-2*n, t*z) is
 ;;
-;;   (t*z)^(1-v-2*n)*exp(-t*z)/CF
+;;   (t*z)^(1-v-2*n)/CF
 ;;
 ;; which finally gives
 ;;
-;; integrate(exp(-t*z*s)*(1+s)^(-2*n-v), s, 0, inf)
-;;   = CF
-;;
-;; and I[n](t, z, v) = exp(-t*z)/t^(2*n+v-1)/CF
-(defun big-i (n t z v)
-  (/ (exp (- (* t z)))
-     (expt t (+ n n v -1))
+;;   integrate(exp(-t*z*s)*(1+s)^(-2*n-v), s, 0, inf)
+;;     = (t*z)^(1-v-2*n)/CF
+;;
+;; and I[n](t, z, v) = exp(-t*z)/CF
+(defun big-i (n theta z v)
+  (/ (exp (- (* theta z)))
      (let* ((a (- 1 v n n))
-            (z-a (- z a)))
+            (z-a (- (* theta z) a)))
        (lentz #'(lambda (n)
                   (+ n n 1 z-a))