fig610.m

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

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

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

% H2PO4-
z1=[1500 1350 1350 1200 1200 1200 900 900 900 600];
z2=[600 600 600 600 300 300 300 300 300 150];
z3=[0 0 0 0 0 0 0];
z=[z1 z2 z3];
zp1=z;

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

% H2PO4-
z1=[1085.8 855.2 831.7 689 672 659.3 634.4 470.3 113.9 440.5];
z2=[320.2 192.2 53.3 20 243.5 43.6 48.8 8.9 10.2 119.8];
z3=[0 0 0 0 0 0 0];
z=[z1 z2 z3];
zp2=z;

% ------------------------------------------------------------------------
%  Ion adsorption mol per g of geothite :  Vmol/g

zp=(zp1-zp2)*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]);
draw27(zp,'P adsptn 10^{-6}mol/g')
set(h2,'Zlim',[0 200])

h3=axes('position',[0.55 0.12 0.32 0.8]);
draw27(zp2,'P Equ. 10^{-6}M')
set(h3,'Zlim',[0 1500])


str={'图6.10 磷吸附量与平衡浓度的对比'}; %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,'fig610.jpeg') % succeed in 2000.10.01
% n=-150 OK