Wrote estimate_new_T function in perturbation.py