The log-factorial function
flog = specfun_factoriallog ( n )
a matrix of floating point integers, must be positive
a matrix of floating point integers
Returns the logarithm of the factorial of n.
This function overflows as soon as x>1.e306.
"Introduction to discrete probabilities", Michael Baudin, 2010