
function Df=jac_protolys1(x,nHA0,nA0,V,pKa)

% 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
      
Df=[Ka -H  0  -A      
    -V -V  0   0
     0  0 -H -OH
     0 -1 -1   1];
