Commit 95aa580ce0dd62113535c30c4c0c82f8656c2d86

Authored by Raymond Toy
1 parent dbc1e376

T is not a variable. Use a different name.

Showing 1 changed file with 7 additions and 8 deletions
qd-bessel.lisp
... ... @@ -320,19 +320,18 @@
320 320 ;;
321 321 ;; So incomplete_gamma_tail(1-v-2*n, t*z) is
322 322 ;;
323   -;; (t*z)^(1-v-2*n)*exp(-t*z)/CF
  323 +;; (t*z)^(1-v-2*n)/CF
324 324 ;;
325 325 ;; which finally gives
326 326 ;;
327   -;; integrate(exp(-t*z*s)*(1+s)^(-2*n-v), s, 0, inf)
328   -;; = CF
  327 +;; integrate(exp(-t*z*s)*(1+s)^(-2*n-v), s, 0, inf)
  328 +;; = (t*z)^(1-v-2*n)/CF
329 329 ;;
330   -;; and I[n](t, z, v) = exp(-t*z)/t^(2*n+v-1)/CF
331   -(defun big-i (n t z v)
332   - (/ (exp (- (* t z)))
333   - (expt t (+ n n v -1))
  330 +;; and I[n](t, z, v) = exp(-t*z)/CF
  331 +(defun big-i (n theta z v)
  332 + (/ (exp (- (* theta z)))
334 333 (let* ((a (- 1 v n n))
335   - (z-a (- z a)))
  334 + (z-a (- (* theta z) a)))
336 335 (lentz #'(lambda (n)
337 336 (+ n n 1 z-a))
338 337 #'(lambda (n)
... ...