3630.sagews

l_1 = plot(e^(-1/2*x),(x,0,6), rgbcolor='blue', legend_label='exp(-1/2*x)')
l_2 = plot(cos(2*pi*x)*e^(-1/2*x),(x,0,6), rgbcolor='yellow', legend_label='exp(-1/2*x)*cos(2*pi*x)')
l = l_1+l_2 
l.set_legend_options(shadow=True) 
l.set_legend_options(loc=(4)) 
l.axes_labels(['x - Achse','y - Achse']) 
l.axes_label_color('green')  
show (l)

k=plot((-x^2+4*x+5),(x,-1,6), rgbcolor=(1,0,0))
show(k)
f= lambda x:-x^2+4*x+5
find_maximum_on_interval(f, 0,5)

(9.0, 2.000000034598127)

a = animate([circle((i/2,abs((10-i/15)*sin(i/3))+i/5), 2/i^(2/3)) for i in srange(1,40)],xmin=-1,ymin=0,xmax=30,ymax=18,aspect_ratio=1,axes=false)
show(a)