13 lines
498 B
Plaintext
13 lines
498 B
Plaintext
(defun gamma (x)
|
|
(let* ((cof #(76.18009172947146 -86.50532032941677
|
|
24.01409824083091 -1.231739572450155
|
|
0.001208650973866179 -0.000005395239384953))
|
|
(ser0 1.000000000190015)
|
|
(x55 (+ x 5.5))
|
|
(tmp (- x55 (* (+ x 0.5) (log x55))))
|
|
(ser (+ ser0 (sum [mapcar / cof (succ x)]))))
|
|
(exp (- (log (/ (* 2.5066282746310005 ser) x)) tmp))))
|
|
|
|
(each ((i (rlist 0.1..1.0..0.1 2..10)))
|
|
(put-line (pic "##.# ######.######" i (gamma i))))
|