fig603.m

这是一个matlab应用的论文和源代码

% NaCl-NaH2PO4-H2O sys summary 2000 08 04 23:15

% ------------------------------------------------------------------------
% ORIGINAL soln, X10^6M
% H+
z1=[6.729 6.165 6.067 5.675 5.584 5.470 4.786 4.632 5.000 3.881];
z2=[3.837 3.784 3.772 4.045 2.421 2.328 2.244 2.517 2.437 1.496];
z3=[0.679 0.495 0.448 0.448 0.495 0.476 0.411];
zh1=[z1 z2 z3];
% ------------------------------------------------------------------------
% BALANCED soln, X10^6M

% H+
z1=[0.06683 0.05284 0.06025 0.04178 0.04786 0.05444 0.02857 0.03564 0.03614 0.02238];
z2=[0.02249 0.02660 0.02679 0.02294 0.02884 0.02393 0.02642 0.02741 0.02779 0.04027];
z3=[0.06591 0.04731 0.03846 0.03523 0.03184 0.03228 0.03435];
z=[z1 z2 z3];
zh2=z;
% ------------------------------------------------------------------------
%  Ion adsorption mol per g of geothite :  Vmol/g
zh=(zh1-zh2)*0.2;

fh=figure('Color','w','position',[200 100 500 350]);
h1=axes('position',[0 0 1 1],'visible','off');
h2=axes('position',[0.07 0.12 0.32 0.8]);

contourdraw27(zh1,'Ini [H^{+}], 10^{-6}M')

h3=axes('position',[0.55 0.12 0.32 0.8]);
contourdraw27(zh2,'Equ [H^{+}], 10^{-6}M')

str={'图6.3 吸附前后溶液中H^{ + }浓度变化（等高线图）'}; %the note of this picture! 
   set(gcf,'currentAxes',h1)
   text(0.5,0.03,str,'Fontsize',12,'horizontalAlignment','center')

[X,map]=capture(fh);
 imwrite(X,map,'fig603.jpeg') % succeed in 2000.10.01