myfuncuo乙二胺'.asv

在化学平衡计算中

global ClT EnT ;
int Row Col;
Row=26;
Col=26;
double Cu[Row][Col];double pH[Row][Col];double En[Row][Col];double Cl[Row][Col];

%double Cu1[Row][Col];double pH1[Row][Col];double En1[Row][Col];double Cl1[Row][Col];

double Fval1[Row][Col];
double Fval2[Row][Col];
double Fval3[Row][Col];
double Fval4[Row][Col];

for i=1:Row
  for j=1:Col
      ClT=j/5-0.2;
      EnT=i/2.5-0.4;
    X0=[4;0.1;0.1;0.1;ClT;EnT];
      %X0=[abs(pH(i,j));abs(Cl(i,j));abs(En(i,j));abs(Cu(i,j));ClT;EnT];
    
   %X0=[pH(i,j);Cl(i,j);NH(i,j);Ni(i,j);ClT;NHT];
    % X0=[pH1(i,j);Cl1(i,j);En1(i,j);Cu1(i,j);ClT;NHT];
   options=optimset('Display','off');
    [X,fval]=fsolve('myfunCuO',X0,options);
    pH(i,j)=X(1);   Cl(i,j)=X(2); En(i,j)=X(3);   Cu(i,j)=X(4);
    %pH1(i,j)=X(1);    Cl1(i,j)=X(2);   En1(i,j)=X(3);     Cu1(i,j)=X(4);

    Fval1(i,j)=fval(1);
    Fval2(i,j)=fval(2);
    Fval3(i,j)=fval(3);
    Fval4(i,j)=fval(4);

end
end