Calculeu el pH d' una dissolució ca = 0.1 M d' àcid H2A 
Dades: k1 = (10^(-4)); k2 = (10^(-5)); ca = 0.1; kw = (10^(-14));

(* Amb Eliminate *)

ClearAll["Global`*"];
eq=Eliminate[{kw==oh*h,
			k1==ha*h/h2a,
			k2==a*h/ha,
			h==oh+ha+2*a,
			ca==h2a+ha+a},{a,ha,h2a,oh}]
k1=(10^(-4));k2=(10^(-5));ca=0.1;kw=(10^(-14));
sol1=Solve[eq,h]
Print["pH = ", -Log[10,sol1[[4,1,2]]]]

(*Amb Eliminate més Solve amb condicions de solució positiva*)

sol2=Solve[{eq,h>0},h]
Print["pH = ", -Log[10,sol2[[1,1,2]]]]

(*Directament Solve i condicions solucions positives*)

ClearAll["Global`*"];
k1=(10^(-4));k2=(10^(-5));ca=0.1;kw=(10^(-14));
sol3=Solve[{kw==oh*h,
			k1==ha*h/h2a,
			k2==a*h/ha,
			h==oh+ha+2*a,
			ca==h2a+ha+a,
			a>0,ha>0,h2a>0,oh>0,h>0},{h,a,ha,h2a,oh}]
Print["pH = ", -Log[10,sol3[[1,1,2]]]]