x*(16*x**7 + 192*x**6 + 1008*x**5 + 3024*x**4 + 5670*x**3 + 6804*x**2 + 5103*x + 2187)
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)
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)
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)
Want a full Python shell? Use SymPy Live.