
%% Fenol

% Redigera dessa värden för olika system:
nHA0=1e-4;       % Ämnesmängd tillsatt syra HA (i mol).
nA0=0;           % Ämnesmängd tillsatt konjugerad bas A (i mol).
V=1;             % Volym av vattenlösningen (i liter).
pKa=9.886;       % Detta är pKa för fenol.
x=1e-5*[1 1 1 1]'; % Startgissning

format short g, format compact
clc, disp('===== Fenol =====')

f=@(x)fun_protolys1(x,nHA0,nA0,V,pKa);
Df=@(x)jac_protolys1(x,nHA0,nA0,V,pKa);
x=newton(f,Df,x,1e-8)
pH=-log10(x(4))

%% 2-Klorofenol

nHA0=7.778e-5;   % Ämnesmängd tillsatt syra HA (i mol).
nA0=0;           % Ämnesmängd tillsatt konjugerad bas A (i mol).
V=1;             % Volym av vattenlösningen (i liter).
pKa=8.48;        % Detta är pKa för 2-Klorofenol.
x=1e-5*[1 1 1 1]'; % Startgissning

clc, disp('===== 2-Klorofenol =====')

f=@(x)fun_protolys1(x,nHA0,nA0,V,pKa);
Df=@(x)jac_protolys1(x,nHA0,nA0,V,pKa);
x=newton(f,Df,x,1e-8)
pH=-log10(x(4))
