Changeset 95aa58 for qdbessel.lisp
 Timestamp:
 04/08/12 16:37:22 (3 years ago)
 Branches:
 master
 Children:
 c7fc98
 Parents:
 dbc1e3
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

qdbessel.lisp
ra5a4c7 r95aa58 321 321 ;; So incomplete_gamma_tail(1v2*n, t*z) is 322 322 ;; 323 ;; (t*z)^(1v2*n) *exp(t*z)/CF323 ;; (t*z)^(1v2*n)/CF 324 324 ;; 325 325 ;; which finally gives 326 326 ;; 327 ;; integrate(exp(t*z*s)*(1+s)^(2*nv), s, 0, inf) 328 ;; = CF 329 ;; 330 ;; and I[n](t, z, v) = exp(t*z)/t^(2*n+v1)/CF 331 (defun bigi (n t z v) 332 (/ (exp ( (* t z))) 333 (expt t (+ n n v 1)) 327 ;; integrate(exp(t*z*s)*(1+s)^(2*nv), s, 0, inf) 328 ;; = (t*z)^(1v2*n)/CF 329 ;; 330 ;; and I[n](t, z, v) = exp(t*z)/CF 331 (defun bigi (n theta z v) 332 (/ (exp ( (* theta z))) 334 333 (let* ((a ( 1 v n n)) 335 (za ( za)))334 (za ( (* theta z) a))) 336 335 (lentz #'(lambda (n) 337 336 (+ n n 1 za))
Note: See TracChangeset
for help on using the changeset viewer.