program done. Need to experiment rho_max and omega_r