isotherm_testify.m

来自「有限差分法（crank_nicholson finite difference 」· M 代码 · 共 6 行

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

fplot('(10.9043*x/(1+0.0056*x))',[0,60],'b');hold on;

fplot('(54.3487*x/(1+1.5267*x)+8.7814*x/(1+0.0045*x))',[0,60],'g');hold on;
legend('1','2')