47.sagews

x = var('x')
diff(x^2, x)

2*x

x = var('x')
p1 = plot(x^2, -5, 5, rgbcolor=hue(0.6))
show(p1)

x = var('x')
p1 = plot(x^2, -5, 5, rgbcolor=hue(0.6))
p2 = plot(2*x, -5, 5, rgbcolor=hue(0.3))
show(p1 + p2)