ADDED: superfactorial s(n) = prod(0..n,n\!)