4318.sagews

var('x')
f =  cos(x)*x^2 
fp = f.plot(xmin=0, xmax=20, rgbcolor=(0,0,.5))

(fp).show()