integrate(x**a*exp(-x), (x, 0, oo))