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)
Want to compute something more complicated? Try a full Python/SymPy console at SymPy Live.