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

Hirst-8-8.5 (last edited 2010-12-18 11:37:12 by SabareeshNikhil)