A universal procedure for pH calculations of equilibrium solutions of one or more protolytes is given. The calculations are based on a mathematical treatment of the chemical model of the equilibrium solution. The method is suitable mainly for didactic purposes.