

def gamma(z, cof=_gammacoeff, g=7):
    z -= 1.0
    v1 = [cof[0]]
    v2 = [(cof[z] / (z + i)) for i in range(1, len(cof))]
    v1 = v1.__add__(v2)
    s = msum(v1)
    z += 0.5
    return (((((z + g) ** z) / exp((z + g))) * sqrt((2.0 * pi))) * s)

dataset/ETHPy150Open anandology/pyjamas/examples/libtest/RandomModuleTest.py gamma/VarMisuse@37/43 `i`->`z`