Распределительная диаграмма соли

[Snaut]

Здравствуйте!
Кто "в теме", проверьте мои рассуждения
мне нужно построить на компе распределительную диаграмму зависимости молярной доли от pH соли.
Как это делать для кислоты, скажем двухосновной - понятно
не сложно и везде есть формулы для мольных долей любого компонента
напр. H2CO3 = H+ + HCO3
HCO3- = CO3 2- + H+
уравнение материального баланса C(CO3 2-)=[HCO3-] + [CO3 2-] + [H2CO3]
далее расписываем K1=... и K2=... по двум уравнениям выше
совершая математ. действия получаем:
a0=[H+]^2/(K1*K2+K1*[H]+[H]^2)
a1=K1[H+]/(K1*K2+K1*[H]+[H]^2)
a2=[CO3 2-]/C(CO3 2-)=K1*K2/(K1*K2+K1*[H]+[H]^2)
получается график примерно как во вложении

Но как это же сделать для того же карбоната натрия?
Если сначала идет полностью диссоциация Na2CO3 = 2Na+ + CO3 2-
потом гидролиз CO3 2- + H2O = HCO3- + OH-
HCO3- + H2O = H2CO3 + OH-
так? но какие и откуда брать K1 и K2?
вместо иона водорода у нас будет фигурировать гидроксид ион? тогда для pH еще надо будет Kw=[H+]*[OH-]=10^-14