The factorial function
f = factorial ( n )
a matrix of doubles. Must contain positive integers.
a matrix of doubles.
Returns the factorial of n, that is, the product of all integers 1 * 2 * ... * n.
This function overflows as soon as n>170.
Introduction to discrete probabilities, Michael Baudin