Computation of multi-component E-pH predominance diagrams