i have the script and i want to add three condition, not sure how to do it.
Kw = 1e-14;Ka = 1.755388e-5;pKa = 4.756;%
Ca = 0.5;Cb = 0.1;Va = 100;Vb = 0.1:0.1:1500;%Ma = (Ca * Va) / 1000;%Mb = (Cb .* Vb) ./ 1000;%for i = 1:length(Mb) M_excess = Ma - Mb(i);
this is where im stuck!!!!!!!!!
if M_excess > 0Ma_final = (M_excess * 1000) ./ (Va + Vb(i)); Mb_final = (Mb(i) * 1000) ./ (Va + Vb(i));Hplus = Ka * (Ma_final ./ Mb_final);if M_excess < 0OH = (M_excess * 1000 * (-1)) ./ (Va + Vb(i));Hplus = Kw ./ OH;if M_excess = 0Hplus = Ka * ((Ma_final * 0.999999) ./ Mb_final); ./ (Va + Vb(i));endpH(i) = -log10(Hplus);end
Best Answer