A brief review of the numerical methods employed in the evaluation of acidity constants of two step overlapping equilibria is given. A new alternative to the method of three equations developed by Thamer is also proposed in this paper. Taking into account the Gauss-elimination method, several variants to the original method of Thamer are devised.