product (i,(i,1,k),(k,1,n))
(Product(factorial(k), (k, 1, n))).doit()
Want a full Python shell? Use SymPy Live.