fig605.m

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

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

% ------------------------------------------------------------------------
% ORIGINAL soln, X10^6M
tic
% Na+
z1=[0.0015 0.00135 0.01135 0.0012 0.0112 0.0212 0.0009 0.0209 0.0409 0.0006];
z2=[0.0106 0.02060 0.04060 0.0606 0.0003 0.0203 0.0403 0.0603 0.0803 0.00015];
z3=[0 0.01 0.02 0.04 0.06 0.08 0.1];
z=[z1 z2 z3];
zna1=z*10^6;

% ------------------------------------------------------------------------
% BALANCED soln, X10^6M

% Na+
z1=[0.0015 0.00135 0.01135 0.0012 0.0112 0.0212 0.0009 0.0209 0.0409 0.0006];
z2=[0.0106 0.02060 0.04060 0.0606 0.0003 0.0203 0.0403 0.0603 0.0803 0.000232];
z3=[0.0000935 0.01 0.02 0.04 0.06 0.08 0.1];
z=[z1 z2 z3];
z=z*10^6;
zna2=z;
% ------------------------------------------------------------------------

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]);

draw27(zna1,'Ini [Na^{+}], 10^{-6}M')

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


str={'图6.5 吸附前后溶液中Na^{ + }浓度变化（三维网格图）'}; %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,'fig605.jpeg') % succeed in 2000.10.01
toc
% n=-150 is OK