Show
Ignore:
Timestamp:
04/08/12 09:37:22 (2 years 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:
1 modified

Legend:

Unmodified
Added
Removed
  • qd-bessel.lisp

    ra5a4c7 r95aa58  
    321321;; So incomplete_gamma_tail(1-v-2*n, t*z) is 
    322322;; 
    323 ;;   (t*z)^(1-v-2*n)*exp(-t*z)/CF 
     323;;   (t*z)^(1-v-2*n)/CF 
    324324;; 
    325325;; which finally gives 
    326326;; 
    327 ;; integrate(exp(-t*z*s)*(1+s)^(-2*n-v), s, 0, inf) 
    328 ;;   = CF 
    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)) 
     327;;   integrate(exp(-t*z*s)*(1+s)^(-2*n-v), s, 0, inf) 
     328;;     = (t*z)^(1-v-2*n)/CF 
     329;; 
     330;; and I[n](t, z, v) = exp(-t*z)/CF 
     331(defun big-i (n theta z v) 
     332  (/ (exp (- (* theta z))) 
    334333     (let* ((a (- 1 v n n)) 
    335             (z-a (- z a))) 
     334            (z-a (- (* theta z) a))) 
    336335       (lentz #'(lambda (n) 
    337336                  (+ n n 1 z-a))