integrate (1/(1+x**2))

integrate(1/(x**2 + 1), x)

plot(atan(x))

solve(atan(x), x)

diff(atan(x), x)

series(atan(x), x, 0, 10)

Experiment with SymPy at SymPy Live.