
function f=fun_protolys1(x,nHA0,nA0,V,pKa)

% Följande ekvationer ska uppfyllas:
% HA <-> A + H    Ka = [A]*[H]/[HA]     Syrans protolysjämvikt
% nHA0 + nA0 = V*([HA] + [A])           Massbalans 
% H2O <-> OH+H    Kw = [OH]*[H]         Vattnets protolysjämvikt
% [H] = [A] - [A]0 + [OH]               Protonbalans (H kommer från vatten
%                                       som bildar OH och HA som bildar A,
%                                       eller omvänt)

% Parametrar för olika system:
% nHA0     Ämnesmängd tillsatt syra HA (i mol).
% nA0      Ämnesmängd tillsatt konjugerad bas A (i mol).
% V        Volym av vattenlösningen (i liter).
% pKa      pKa för aktuell fenol.

HA=x(1);   % Koncentrationen av syran HA i mol/liter vid jämvikt
A=x(2);    % Koncentrationen av den konjugerade basen A i mol/liter vid jämvikt
OH=x(3);   % Hydroxidjonkoncentrationen i mol/liter vid jämvikt
H=x(4);    % Vätejonkoncentrationen i mol/liter vid jämvikt

pKw=14;         
Ka=10^(-pKa);
Kw=10^(-pKw);     % Vattnets jonprodukt vid 25C

f=[Ka*HA-A*H           % Syrans protolysjämvikt
   nHA0+nA0-V*(HA+A)   % Balans av HA 
   Kw-OH*H             % Vattnets protolysjämvikt
   H+nA0/V-(A+OH)];    % Balans av H
