3722.sagews

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

43 ! = 60415263063373835637355132068513997507264512000000000 TEST: 60415263063373835637355132068513997507264512000000000

l = [] 
for i in srange(1,2,0.99):
     l.append(circle((2,abs(-1*i)), 1) ) 
a= animate(l, xmin=0,ymin=0,xmax=4,ymax=4,
figsize=[4,4])
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)