哪位大大有1stopt2.0以上版本，能不能帮我算算下面问题
parameter z1,z2,z3,z5,z6,z7;ConstStr xH2O1=11.04911698/(0.058639389*(1+1/17+(exp(z1+z2/272.15+z3/272.15^2)*1.0521e+005)/(1+exp(z1+z2/272.15+z3/272.15^2)*1.0521e+005)*2/17)+ 11.04911698*(1+1/17+(7.7299e-004*5.2719e+003)/(1+7.7299e-004*5.2719e+003)*2/17));
ConstStr xH2O21=0.058639389/(0.058639389*(1+1/17+(exp(z1+z2/272.15+z3/272.15^2)*1.0521e+005)/(1+exp(z1+z2/272.15+z3/272.15^2)*1.0521e+005)*2/17)+ 11.04911698*(1+1/17+(7.7299e-004*5.2719e+003)/(1+7.7299e-004*5.2719e+003)*2/17));
ConstStr yH2O21=xH2O21*(1/17+(exp(z1+z2/272.15+z3/272.15^2)*1.0521e+005)/(1+exp(z1+z2/272.15+z3/272.15^2)*1.0521e+005)*2/17)/(xH2O21*(1/17+(exp(z1+z2/272.15+z3/272.15^2)*1.0521e+005)/(1+exp(z1+z2/272.15+z3/272.15^2)*1.0521e+005)*2/17)+xH2O1*(1/17+(7.7299e-004*5.2719e+003)/(1+7.7299e-004*5.2719e+003)*2/17));
ConstStr xH2O2=3.972661765/(0.232587743*(1+1/17+(exp(z1+z2/272.15+z3/272.15^2)*2.0418e+005)/(1+exp(z1+z2/272.15+z3/272.15^2)*2.0418e+005)*2/17)+ 3.972661765*(1+1/17+(7.7299e-004*7.2066e+003)/(1+7.7299e-004*7.2066e+003)*2/17));
ConstStr xH2O22=0.232587743/(0.232587743*(1+1/17+(exp(z1+z2/272.15+z3/272.15^2)*2.0418e+005)/(1+exp(z1+z2/272.15+z3/272.15^2)*2.0418e+005)*2/17)+ 3.972661765*(1+1/17+(7.7299e-004*7.2066e+003)/(1+7.7299e-004*7.2066e+003)*2/17));
ConstStr yH2O22=xH2O22*(1/17+(exp(z1+z2/272.15+z3/272.15^2)*2.0418e+005)/(1+exp(z1+z2/272.15+z3/272.15^2)*2.0418e+005)*2/17)/(xH2O22*(1/17+(exp(z1+z2/272.15+z3/272.15^2)*2.0418e+005)/(1+exp(z1+z2/272.15+z3/272.15^2)*2.0418e+005)*2/17)+xH2O2*(1/17+(7.7299e-004*7.2066e+003)/(1+7.7299e-004*7.2066e+003)*2/17));
ConstStr xH2O3=8.573270588/(0.067891358*(1+1/17+(exp(z1+z2/272.15+z3/272.15^2)*1.3508e+005)/(1+exp(z1+z2/272.15+z3/272.15^2)*1.3508e+005)*2/17)+ 8.573270588*(1+1/17+(7.7299e-004*5.7826e+003)/(1+7.7299e-004*5.7826e+003)*2/17));
ConstStr xH2O23=0.067891358/(0.067891358*(1+1/17+(exp(z1+z2/272.15+z3/272.15^2)*1.3508e+005)/(1+exp(z1+z2/272.15+z3/272.15^2)*1.3508e+005)*2/17)+ 8.573270588*(1+1/17+(7.7299e-004*5.7826e+003)/(1+7.7299e-004*5.7826e+003)*2/17));
ConstStr yH2O23=xH2O23*(1/17+(exp(z1+z2/272.15+z3/272.15^2)*1.3508e+005)/(1+exp(z1+z2/272.15+z3/272.15^2)*1.3508e+005)*2/17)/(xH2O23*(1/17+(exp(z1+z2/272.15+z3/272.15^2)*1.3508e+005)/(1+exp(z1+z2/272.15+z3/272.15^2)*1.3508e+005)*2/17)+xH2O3*(1/17+(7.7299e-004*5.7826e+003)/(1+7.7299e-004*5.7826e+003)*2/17));
ConstStr xH2O4=11.08923335/(0.19612037*(1+1/17+(exp(z1+z2/273.15+z3/273.15^2)*1.6343e+005)/(1+exp(z1+z2/273.15+z3/273.15^2)*1.6343e+005)*2/17)+ 11.08923335*(1+1/17+(7.2746e-004*6.5278e+003)/(1+7.2746e-004*6.5278e+003)*2/17));
ConstStr xH2O24=0.19612037/(0.19612037*(1+1/17+(exp(z1+z2/273.15+z3/273.15^2)*1.6343e+005)/(1+exp(z1+z2/273.15+z3/273.15^2)*1.6343e+005)*2/17)+ 11.08923335*(1+1/17+(7.2746e-004*6.5278e+003)/(1+7.2746e-004*6.5278e+003)*2/17));
ConstStr yH2O24=xH2O24*(1/17+(exp(z1+z2/273.15+z3/273.15^2)*1.6343e+005)/(1+exp(z1+z2/273.15+z3/273.15^2)*1.6343e+005)*2/17)/(xH2O24*(1/17+(exp(z1+z2/273.15+z3/273.15^2)*1.6343e+005)/(1+exp(z1+z2/273.15+z3/273.15^2)*1.6343e+005)*2/17)+xH2O4*(1/17+(7.2746e-004*6.5278e+003)/(1+7.2746e-004*6.5278e+003)*2/17));
ConstStr xH2O5=5.119722253/(0.212892563*(1+1/17+(exp(z1+z2/273.15+z3/273.15^2)*1.9124e+005)/(1+exp(z1+z2/273.15+z3/273.15^2)*1.9124e+005)*2/17)+ 5.119722253*(1+1/17+(7.2746e-004*7.4119e+003)/(1+7.2746e-004*7.4119e+003)*2/17));
ConstStr xH2O25=0.212892563/(0.212892563*(1+1/17+(exp(z1+z2/273.15+z3/273.15^2)*1.9124e+005)/(1+exp(z1+z2/273.15+z3/273.15^2)*1.9124e+005)*2/17)+ 5.119722253*(1+1/17+(7.2746e-004*7.4119e+003)/(1+7.2746e-004*7.4119e+003)*2/17));
ConstStr yH2O25=xH2O25*(1/17+(exp(z1+z2/273.15+z3/273.15^2)*1.9124e+005)/(1+exp(z1+z2/273.15+z3/273.15^2)*1.9124e+005)*2/17)/(xH2O25*(1/17+(exp(z1+z2/273.15+z3/273.15^2)*1.9124e+005)/(1+exp(z1+z2/273.15+z3/273.15^2)*1.9124e+005)*2/17)+xH2O5*(1/17+(7.2746e-004*7.4119e+003)/(1+7.2746e-004*7.4119e+003)*2/17));
ConstStr xH2O6=8.939479578/(0.25310341*(1+1/17+(exp(z1+z2/274.15+z3/274.15^2)*1.7712e+005)/(1+exp(z1+z2/274.15+z3/274.15^2)*1.7712e+005)*2/17)+ 8.939479578*(1+1/17+(6.8492e-004*7.7655e+003)/(1+6.8492e-004*7.7655e+003)*2/17));
ConstStr xH2O26=0.25310341/(0.25310341*(1+1/17+(exp(z1+z2/274.15+z3/274.15^2)*1.7712e+005)/(1+exp(z1+z2/274.15+z3/274.15^2)*1.7712e+005)*2/17)+ 8.939479578*(1+1/17+(6.8492e-004*7.7655e+003)/(1+6.8492e-004*7.7655e+003)*2/17));
ConstStr yH2O26=xH2O26*(1/17+(exp(z1+z2/274.15+z3/274.15^2)*1.7712e+005)/(1+exp(z1+z2/274.15+z3/274.15^2)*1.7712e+005)*2/17)/(xH2O26*(1/17+(exp(z1+z2/274.15+z3/274.15^2)*1.7712e+005)/(1+exp(z1+z2/274.15+z3/274.15^2)*1.7712e+005)*2/17)+xH2O6*(1/17+(6.8492e-004*7.7655e+003)/(1+6.8492e-004*7.7655e+003)*2/17));
Function (1.0521e+005-(exp(1.13e+005*yH2O21*1.0224e-005/272.15)*z5*exp(z6/(272.15-z7))*0.0053^((-1)/(1/17)))*(1-(exp(z1+z2/272.15+z3/272.15^2)*1.0521e+005)/(1+exp(z1+z2/272.15+z3/272.15^2)*1.0521e+005))^2)=0;
(2.0418e+005-(exp(2.21e+005*yH2O22*1.0224e-005/272.15)*z5*exp(z6/(272.15-z7))*0.0453^((-1)/(1/17)))*(1-(exp(z1+z2/272.15+z3/272.15^2)*2.0418e+005)/(1+exp(z1+z2/272.15+z3/272.15^2)*2.0418e+005))^2)=0;
(1.3508e+005-(exp(1.45e+005*yH2O23*1.0224e-005/272.15)*z5*exp(z6/(272.15-z7))*0.0176^((-1)/(1/17)))*(1-(exp(z1+z2/272.15+z3/272.15^2)*1.3508e+005)/(1+exp(z1+z2/272.15+z3/272.15^2)*1.3508e+005))^2)=0;
(1.6343e+005-(exp(1.76e+005*yH2O24*1.0224e-005/273.15)*z5*exp(z6/(273.15-z7))*0.0176^((-1)/(1/17)))*(1-(exp(z1+z2/273.15+z3/273.15^2)*1.6343e+005)/(1+exp(z1+z2/273.15+z3/273.15^2)*1.6343e+005))^2)=0;
(1.9124e+005-(exp(2.07e+005*yH2O25*1.0224e-005/273.15)*z5*exp(z6/(273.15-z7))*0.0339^((-1)/(1/17)))*(1-(exp(z1+z2/273.15+z3/273.15^2)*1.9124e+005)/(1+exp(z1+z2/273.15+z3/273.15^2)*1.9124e+005))^2)=0;
(1.7712e+005-(exp(1.92e+005*yH2O26*1.0224e-005/274.15)*z5*exp(z6/(274.15-z7))*0.0247^((-1)/(1/17)))*(1-(exp(z1+z2/274.15+z3/274.15^2)*1.7712e+005)/(1+exp(z1+z2/274.15+z3/274.15^2)*1.7712e+005))^2)=0;
式子长了些，没办法，初学者不知道怎么样简单表示，其实就只有一个总表达式，那些数都可以用向量表示，我用的是1.5版，也不太会设置算法，所以目标函数值总是在60万以上，但我希望目标函数值在0附近。
页:
[1]