3647.sagews

p1=plot(e^(-1/2*x)*cos(2*pi*x),(0,6),rgbcolor=(0,1,0),legend_label='$e^{(-1/2*x)}*cos(2*\pi*x)$')
p2=plot(e^(-1/2*x),(0,6),rgbcolor=(1,0,0),legend_label='$e^{(-1/2*x)}$') 
p=p1+p2 
p.axes_labels(['$x$','$f(x)$']) 
p.axes_label_color((0,0,1)) 
p.set_legend_options(shadow=True) 
p.set_legend_options(loc=(4)) 
show(p)

f = (x+1)
find_maximum_on_interval(f, -1,3)

(3.99999992962, 2.9999999296172466)

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