Create factorize.py