Changeset c0e12ddf6b61f571555d46c6f168e6bebabc80b1

Show
Ignore:
Timestamp:
04/11/12 20:06:38 (2 years ago)
Author:
Raymond Toy <toy.raymond@…>
Children:
d795ba718dc53f591c82994811f50250aceec1d7
Parents:
7c5a3186070096ee93e16f2ddf51b2c84e7c5895
git-committer:
Raymond Toy <toy.raymond@…> (04/11/12 20:06:38)
Message:

Correct the computation of alpha and beta.

Files:
1 modified

Legend:

Unmodified
Added
Removed
  • qd-bessel.lisp

    r7c5a31 rc0e12d  
    214214(defun alpha (n z) 
    215215  (let ((n (float n (realpart z)))) 
    216     (/ (cf-incomplete-gamma (1+ n) (/ z 2)) 
     216    (/ (incomplete-gamma (1+ n) (/ z 2)) 
    217217       (expt z (1+ n))))) 
    218218 
    219219(defun beta (n z) 
    220220  (let ((n (float n (realpart z)))) 
    221     (/ (- (cf-incomplete-gamma (1+ n) (/ z 2)) 
    222           (cf-incomplete-gamma (1+ n) (/ z -2))) 
     221    (/ (- (incomplete-gamma (1+ n) (/ z 2)) 
     222          (incomplete-gamma (1+ n) (/ z -2))) 
    223223       (expt z (1+ n))))) 
    224224