Commit c0e12ddf authored by Raymond Toy's avatar Raymond Toy
Browse files

Correct the computation of alpha and beta.

parent 7c5a3186
...@@ -213,13 +213,13 @@ ...@@ -213,13 +213,13 @@
(defun alpha (n z) (defun alpha (n z)
(let ((n (float n (realpart z)))) (let ((n (float n (realpart z))))
(/ (cf-incomplete-gamma (1+ n) (/ z 2)) (/ (incomplete-gamma (1+ n) (/ z 2))
(expt z (1+ n))))) (expt z (1+ n)))))
(defun beta (n z) (defun beta (n z)
(let ((n (float n (realpart z)))) (let ((n (float n (realpart z))))
(/ (- (cf-incomplete-gamma (1+ n) (/ z 2)) (/ (- (incomplete-gamma (1+ n) (/ z 2))
(cf-incomplete-gamma (1+ n) (/ z -2))) (incomplete-gamma (1+ n) (/ z -2)))
(expt z (1+ n))))) (expt z (1+ n)))))
;; a[0](k,v) := (k+sqrt(k^2+1))^(-v); ;; a[0](k,v) := (k+sqrt(k^2+1))^(-v);
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment