Changeset c0e12ddf6b61f571555d46c6f168e6bebabc80b1
- Timestamp:
- 04/11/12 20:06:38 (14 months 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:
-
Legend:
- Unmodified
- Added
- Removed
-
|
r7c5a31
|
rc0e12d
|
|
| 214 | 214 | (defun alpha (n z) |
| 215 | 215 | (let ((n (float n (realpart z)))) |
| 216 | | (/ (cf-incomplete-gamma (1+ n) (/ z 2)) |
| | 216 | (/ (incomplete-gamma (1+ n) (/ z 2)) |
| 217 | 217 | (expt z (1+ n))))) |
| 218 | 218 | |
| 219 | 219 | (defun beta (n z) |
| 220 | 220 | (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))) |
| 223 | 223 | (expt z (1+ n))))) |
| 224 | 224 | |
