def integrate_monomial_over_enr2(k): if numpy.any(k % 2 == 1): return 0 return numpy.prod([math.gamma((kk+1) / 2.0) for kk in k])