All published worksheets from http://sagenb.org
f = x*(x-1)*(x-2)+1
p = plot(f, 0.5, 1.5, fill=True) + plot(f, 0, 2)
show(p)
