3708.sagews

factorial(38)

523022617466601111760007224100074291200000000

a=1
for i in range (38):
    b=a*(i+1)
    a=b
print a

523022617466601111760007224100074291200000000

a=[]
for i in srange(-5,5,0.6):
    a.append(circle((i,abs(-i^2+5)), 0.5))
show(animate(a,xmin=-8,ymin=-5,xmax=2.5,ymax=8,axes=false))

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