It's a shame factorial isn't consider "basic" enough to be included in NumPy.