[Numpy-discussion] how to compute a Gamma function with Numpy?