3629.sagews

p1 = plot(exp(-1/2*x)*cos(2*pi*x),(x,0,6), rgbcolor=(1,0.75,0), legend_label='exp(-1/2*x)*cos(2*pi*x)')
p2 = plot(exp(-1/2*x),(x,0,6),rgbcolor=(1,0,0), legend_label='exp(-1/2*x)')
p = p1 + p2
p.set_legend_options(loc = 4, shadow=True)
p.axes_label_color((0,0,1))
p.axes_labels(['$x$','$y$'])
show (p)

f = (x^3)
find_maximum_on_interval(f,-1,1.1)

(1.33099985888, 1.0999999611244449)

a = animate([circle((1,i+0.5), 0.5, hue=i/10) for i in srange(0,1,0.1)], 
               xmin=0,ymin=0,xmax=2,ymax=2,figsize=[2,2])
a.show()

Animation with 21 frames