Changeset 95aa58 for qd-bessel.lisp


Ignore:
Timestamp:
04/08/12 16:37:22 (3 years ago)
Author:
Raymond Toy <toy.raymond@…>
Branches:
master
Children:
c7fc98
Parents:
dbc1e3
Message:

T is not a variable. Use a different name.

File:
1 edited

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))
Note: See TracChangeset for help on using the changeset viewer.