3609.sagews

x = var('x')
c1 = plot(e^(-x/2), (x,0,6), rgbcolor=(0,0,1), legend_label='e^(-x/2)')
c2 = plot((e^(-x/2))*cos(2*pi*x), (x,0,6), rgbcolor=(3,2,5), legend_label='(e^(-x/2))*cos(2*pi*x)')
c = c1+c2
c.set_legend_options(shadow=True)
c.set_legend_options(loc=(4))
show(c, axes_labels=(['x', 'f(x)']))

f = lambda x: cos(x)
find_maximum_on_interval(f, -1,1)

(1.0, 9.8666666288175245e-09)

a = animate([circle((2,(abs(-1*i))+1), 1-1/(6), hue=i/10) for i in srange(-5,5,0.6)], xmin=0,ymin=0,xmax=10,ymax=10,figsize=[2,2])
show(a)