The object of the present paper is to introduce a new approach, resting on a numerical analysis of the roots of the β function, which allows finding the optimal mathematical conditions to produce flat buffer index curves. It has been applied to solutions of 2, 3 and 4 weak acids, which are simple enough to be analyzed without numerical optimization techniques.