計算海水的pH值

Hp = (-K1*pCO2 - √(K1^2*pCO2^2 - 4*K1*K2*pCO2*(pCO2 - DIC/K0)))/(2*(pCO2 - DIC/K0));
pH = -log(Hp);
H₂CO₃ = pCO₂*K₀;
HCO₃⁻ = H₂CO₃*K₁/Hp;
CO₃²⁻ = HCO₃⁻*K₂/Hp;

print(Hp);
print(pH);
print(H₂CO₃);
print(HCO₃⁻);
print(CO₃²⁻);
print(H₂CO₃*10⁶);
print(HCO₃⁻*10⁶);
print(CO₃²⁻*10⁶);
print(DIC*10⁶);
print(H₂CO₃/DIC*100);
print(HCO₃⁻/DIC*100);
print(CO₃²⁻/DIC*100);