%
clc;
clear?all

close?all

cu=xlsread（‘asxishu‘‘sheet1‘‘A2:A66‘）;??X1=cu;
xishu=xlsread（‘asxishu‘‘sheet1‘‘B2:B66‘）;??Y=xishu;
%diye=xlsread（‘cusolid‘‘sheet1‘‘C2:C66‘）;??Y=diye;


x=[X1?];?
y=Y;

P0=[0.02?0.03?];?%?P0的第一個的ORP的系數，第二個是入口鈷（取對數后）的系數???1.0.*x（:3）
%??f=@（px）???x（:2）.*p（1）.*exp（-p（2）.*（x（:1）*（430/150）））?原，430是反應器體積，150是流量
f=@（px）??p（1）*log（x（:1）.*1）+p（2）%（x（:1）.^2）.*p（1）+x（:1）.*p（2）+p（3）