3633.sagews

x = var('x')
p1 = plot((e^(-0.5*x)*cos(2*pi*x)),(x,0,6),rgbcolor=hue(0.4), legend_label='e^(-0.5*x)*cos(2*pi*x)')
p2 = plot((e^(-0.5*x)),(x,0,6),rgbcolor=hue(0.6), legend_label='e^(-0.5*x)')
x = text('x axis', (6,-0.1))
y = text('y axis', (0.5,1))
p = p1+p2+x+y
p.set_legend_options(shadow=true)
p.set_legend_options(loc=4)
show(p)

x = var('x')
f = sin(x)
find_maximum_on_interval(f,-0.5,1)

(0.841470971212, 0.9999999748357512)

f = animate([circle((1,i+1), 1-8/(i+6)) for i in srange(0.1,2)],
xmin=0,ymin=0,xmax=6,ymax=3,figsize=[5,5])
f.show()