exp(-x)*sum(k=0 to n,n!/(n-k)!*x^(n-k)) + integral(exp(-x)*x^n,x)