myfuncuo123.m

来自「在化学平衡计算中」· M 代码 · 共 13 行

function F=myfunCuO123(X)

global pH EnT ;
F(1)=X(1)-exp(22.8642-2.303*pH+log(X(3)));

F(2)=X(2)-exp(38.6351-4.606*pH+log(X(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);