x**7-x**5+4*x**3+1
solve(x**7 - x**5 + 4*x**3 + 1, x)
diff(x**7 - x**5 + 4*x**3 + 1, x)
series(x**7 - x**5 + 4*x**3 + 1, x, 0, 10)
Want a full Python shell? Use SymPy Live.