SymPy Gamma logo SymPy Gamma

integrate (x**a*exp (-x),(x,0,oo ))
integrate(x**a*exp(-x), (x, 0, oo))
plot(Piecewise((gamma(a + 1), re(a) > -1), (Integral(x**a*exp(-x), (x, 0, oo)), True)))
solve(Piecewise((gamma(a + 1), re(a) > -1), (Integral(x**a*exp(-x), (x, 0, oo)), True)), x)
diff(Piecewise((gamma(a + 1), re(a) > -1), (Integral(x**a*exp(-x), (x, 0, oo)), True)), x)
series(Piecewise((gamma(a + 1), re(a) > -1), (Integral(x**a*exp(-x), (x, 0, oo)), True)), x, 0, 10)
See what Wolfram|Alpha has to say.

Want to compute something more complicated? Try a full Python/SymPy console at SymPy Live.