myfuncuo区域.asv

在化学平衡计算中

function F=myfunCuO123'(X)

global pH EnT ;
F(1)=X(1)-exp(22.8042-2.303pH+lnX(3));

F(2)=X(2)-exp(38.6351-4.606pH+lnX(3));

F(3)=X(1)+X(2)+X(3)-EnT;


%X0=[6.0;0.100;0.0500;1.00;0.1;0.1];
%options=optimset('Display','iter');
%[X,fval]=fsolve('mynewfun',X0,options);