3746.sagews

n=44
f=1
for i in range (1, n+1):
    f= f * i
print n, '! =', f
print 'TEST:'
factorial(44)

44 ! = 2658271574788448768043625811014615890319638528000000000 TEST: 2658271574788448768043625811014615890319638528000000000

l = []
for i in srange (0,10,0.1):
    l.append (circle((i,1+abs(3*sin(i))), 1) )
a = animate (l, xmin=-1, ymin=0, xmax=11, ymax=5, figsize=[12,5])
show(a)

x_coords = [i for i in srange(1, 150, 1)] 
y_coords = [sin(2*0.99*pi*i) for i in srange(1,150, 1)] 
list_plot(zip(x_coords, y_coords))

x_coords = [0.99*i for i in srange(1, 150, 1)] 
y_coords = [sin(2*0.99*pi*i) for i in srange(1,150, 1)] 
l=list_plot(zip(x_coords, y_coords)) 
p=plot(sin(2*pi*x), (x,0,150)) 

show (l+p)