isotherm_testify.asv

有限差分法（crank_nicholson finite difference ) 求解色谱模型

fplot('13.7257*x/(1+0.0196*x)+0.1058*x',[0,60],'r');hold on;

fplot('54.3487*x/(1+1.5267*x)+8.7814*x/(1+0.0045*x)',[0,60],'r');hold on;
a=[0,0;0.114,1.268430333;0.306,3.784844923;0.512,5.989867194;1.08,7.097584304;2.8,8.010158349;8.14,12.391429;10.41,13.15064521;19.62,16.81114992;29.6,21.08834174;49.7,28.09674423;52.97,28.98944838;];
plot(a(:,1),a(:,2),'ko');
legend('1','2','3','experimental data')