This function is responsible to compute the entire equilibrium problem. It works in the following way.
First, the control factor that should limit the speed of
convergence to avoid negative mol number for gazes.
if it is a non-fixed temperature problem,
If there is condensed species which concentration is negative, they should be discard from the list of considered species and convergence to a new equilibrium is obtain.