3716.sagews

s = 45
m = 1 
for i in srange(1,s+1,1): 
    m = m*i 
print(m)

119622220865480194561963161495657715064383733760000000000

factorial(45)

119622220865480194561963161495657715064383733760000000000

s=[] 
for i in srange(-5,5,0.6):
     s.append(circle((2,(abs(-1*i))+1), 1-1/(6)))
show(animate(s,xmin=0,ymin=0,xmax=10,ymax=10))

a = [sin(2*pi*0.99*(n+1)) for n in range (150)]
list_plot(a)

t = var('t') 
a = plot(sin(2*pi*t), (t,0,150)) 
b = [(0.99*n) for b in range(150)] 
c = [(sin(0.99*n*2*pi)) for c in range(150)] 
list_plot((b), rgbcolor=(4,5,9)),show(a)
list_plot((c), rgbcolor=(5,3,3))