SymPy Gamma

SymPy:

integrate ((2*x+3)**7)

$\int \left(2 x + 3\right)^{7}\, dx$

x*(16*x**7 + 192*x**6 + 1008*x**5 + 3024*x**4 + 5670*x**3 + 6804*x**2 + 5103*x + 2187)

$x \left(16 x^{7} + 192 x^{6} + 1008 x^{5} + 3024 x^{4} + 5670 x^{3} + 6804 x^{2} + 5103 x + 2187\right)$

integrate((2*x + 3)**7, x)

Plot:

plot(16*x**8 + 192*x**7 + 1008*x**6 + 3024*x**5 + 5670*x**4 + 6804*x**3 + 5103*x**2 + 2187*x)

Roots:

solve(16*x**8 + 192*x**7 + 1008*x**6 + 3024*x**5 + 5670*x**4 + 6804*x**3 + 5103*x**2 + 2187*x, x)

$x =$

diff(16*x**8 + 192*x**7 + 1008*x**6 + 3024*x**5 + 5670*x**4 + 6804*x**3 + 5103*x**2 + 2187*x, x)

$\frac{d}{d x} \left(16 x^{8} + 192 x^{7} + 1008 x^{6} + 3024 x^{5} + 5670 x^{4} + 6804 x^{3} + 5103 x^{2} + 2187 x\right) =$

series(16*x**8 + 192*x**7 + 1008*x**6 + 3024*x**5 + 5670*x**4 + 6804*x**3 + 5103*x**2 + 2187*x, x, 0, 10)

See what Wolfram|Alpha has to say.

Experiment with SymPy at SymPy Live.