Search All of the Math Forum:

Views expressed in these public forums are not endorsed by NCTM or The Math Forum.

Notice: We are no longer accepting new posts, but the forums will continue to be readable.

Topic: matlab novice
Replies: 12   Last Post: Aug 1, 2013 7:42 AM

 Messages: [ Previous | Next ]
 Torsten Hennig Posts: 2,419 Registered: 12/6/04
Re: matlab novice
Posted: Aug 1, 2013 4:59 AM

I don't know which values you want to keep for further use; I assume it's only the pH-array.

Kw = 1e-14;
Ka = 1.755388e-5;
pKa = 4.756;
%
Ca = 0.5;
Cb = 0.1;
Va = 100;
Vb = 0:0.1:1500;
%
% 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 M_excess > 0
% This simulates final molarity of CH3CO2H.
Ma_final = (M_excess * 1000) / (Va + Vb);
% This simulates final molarity of CH3CO2-.
Mb_final = (Mb(i) * 1000) / (Va + Vb);
% This simulates HPlus ions.
Hplus = Ka * (Ma_final / Mb_final);
else
OH = (M_excess * 1000 * (-1)) / (Va + Vb);
% This simulates HPlus ions.
Hplus = Kw / OH;
end
% This simulates the pH content.
pH(i) = -log10(Hplus);
end

Best wishes
Torsten.

Date Subject Author
8/1/13 harley
8/1/13 Torsten Hennig
8/1/13 harley
8/1/13 harley
8/1/13 Torsten Hennig
8/1/13 harley
8/1/13 harley
8/1/13 harley
8/1/13 Torsten Hennig
8/1/13 harley
8/1/13 Torsten Hennig
8/1/13 harley
8/1/13 harley