Changeset c0e12d


Ignore:
Timestamp:
04/12/12 03:06:38 (3 years ago)
Author:
Raymond Toy <toy.raymond@…>
Branches:
master
Children:
d795ba
Parents:
7c5a31
Message:

Correct the computation of alpha and beta.

File:
1 edited

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