\int x^n e^{ax} dx = (-1)^n\frac{1}{a}\Gamma[1+n,-ax], \text{ where } \Gamma(a,x)=\int_x^{\infty} t^{a-1}e^{-t} dt