3902.sagews

f = x^2 + 2*x
g = f.diff(x)

\int^\infty_0

print latex(g)

2 \, x + 2

I = integrate(f, x)
print I

1/3*x^3 + x^2