3718.sagews

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

print 'TEST:' 

factorial(42)

42 ! = 1405006117752879898543142606244511569936384000000000 TEST: 1405006117752879898543142606244511569936384000000000

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)