A new computer programme for the computation of formation constants of complex species in solution is presented. The programme can be applied to all kinds of potentiometric titration data, including multi-reactant and multi-electrode systems. A statistical analysis is made of the results in order to assess their validity, and to assist in hypothesis testing.