Create rho_max.py