Create estimation.py