add_equilibrium function

Compute equilibrium values