Drexel dragonThe Math ForumDonate to the Math Forum



Search All of the Math Forum:

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


Math Forum » Discussions » Software » comp.soft-sys.matlab

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

Advanced Search

Back to Topic List Back to Topic List Jump to Tree View Jump to Tree View   Messages: [ Previous | Next ]
Torsten Hennig

Posts: 2,419
Registered: 12/6/04
Re: matlab novice
Posted: Aug 1, 2013 4:59 AM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply

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.



Point your RSS reader here for a feed of the latest messages in this topic.

[Privacy Policy] [Terms of Use]

© Drexel University 1994-2014. All Rights Reserved.
The Math Forum is a research and educational enterprise of the Drexel University School of Education.