Return the factorial of n where n is a positive integer. If
n is a scalar, this is equivalent to prod (1:n). For
vector or matrix arguments, return the factorial of each element in the
array. For non-integers see the generalized factorial function
gamma.