https://trspos.com/factorial-en-python-numpy/