follow script, how do i display the value of Vb and pH when "abs(M_excess)<eps" and following that how would i plot that point a figure.
Kw = 1e-14;Ka = 1.755388e-5;pKa = 4.756;%
Ca = 0.5;Cb = 0.1;Va = 100;Vb = 0.05:0.05:2000;%% This simulates initial molarity of CH3CO2H.
Ma = (Ca * Va) / 1000;% This simulates initial molarity of NaOH.
Mb = (Cb .* Vb) ./ 1000;for i = 1:length(Mb) M_excess = Ma - Mb(i);if abs(M_excess)<epsHplus = Ka * ((Ma_final * 0.999999) ./ Mb_final);elseif M_excess > 0 Ma_final = (M_excess * 1000) ./ (Va + Vb(i)); Mb_final = (Mb(i) * 1000) ./ (Va + Vb(i)); Hplus = Ka * (Ma_final ./ Mb_final);elseif M_excess < 0 OH = (M_excess * 1000 * (-1)) ./ (Va + Vb(i)); Hplus = Kw ./ OH;endpH(i) = -log10(Hplus);end
Best Answer