(defun gamma (x)
  (* (sqrt (/ (* 2 %pi%)
              x))
     (expt (/ x %e%) x)))

(each ((i 1..11))
  (put-line (pic "##.######" (gamma (/ i 3.0)))))