Change compute_equilibrium to accept bounds