def gamma(z, cof=_gammacoeff, g=7): z -= 1.0 # Next line fails when not compiled with --operator-funcs #s = msum([cof[0]] + [cof[i] / (z+i) for i in range(1,len(cof))]) v1 = [cof[0]] v2 = [cof[i] / (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
def gamma(z, cof=_gammacoeff, g=7): z -= 1.0 # Next line fails when not compiled with --operator-funcs #s = msum([cof[0]] + [cof[i] / (z+i) for i in range(1,len(cof))]) v1 = [cof[0]] v2 = [cof[i] / (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
def gamma(z, cof=_gammacoeff, g=7): z -= 1.0 s = msum([cof[0]] + [cof[i] / (z+i) for i in range(1,len(cof))]) z += 0.5 return (z+g)**z / exp(z+g) * sqrt(2.0*pi) * s
def gamma(z, cof=_gammacoeff, g=7): z -= 1.0 s = msum([cof[0]] + [cof[i] / (z + i) for i in range(1, len(cof))]) z += 0.5 return (z + g)**z / exp(z + g) * sqrt(2.0 * pi) * s