var('x n')
print "gamma(n) =",integrate(x**n*e**x,x,0,1)