3710.sagews

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

1 2 6 24 120 720 5040 40320 362880 3628800 39916800 479001600 6227020800 87178291200 1307674368000 20922789888000 355687428096000 6402373705728000 121645100408832000 2432902008176640000 51090942171709440000 1124000727777607680000 25852016738884976640000 620448401733239439360000 15511210043330985984000000 403291461126605635584000000 10888869450418352160768000000 304888344611713860501504000000 8841761993739701954543616000000 265252859812191058636308480000000 8222838654177922817725562880000000 263130836933693530167218012160000000 8683317618811886495518194401280000000 295232799039604140847618609643520000000 10333147966386144929666651337523200000000 371993326789901217467999448150835200000000

factorial(36)

371993326789901217467999448150835200000000

a=[] 
for i in srange(0,2,0.8): 
    a.append(circle((1,i), 1-1/(i+30))) 
    show(animate(a,xmin=0,ymin=0,xmax=2,ymax=3,axes=false))

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

b = [0.99*(n+1) for i in range (150),sin(2*pi*0.99*(n+1))for i in range (150), plot(2*pi*x)] 
list_plot (b)