Features and examples:
from scipy.special import gammaln [docs]def logfactorial(n): return gammaln(n + 1)