方程如下:
T=1135;
Ptot=346;
P0=101;
Uf=0.704;
syms x y z;
Kr=exp(-2.63121*10^(-11)*T^4+1.24065*10^(-7)*T^3-2.25232*10^(-4)*T^2+1.95028*10^(-1)*T-66.1395);
Ks=exp(5.47301*10^(-12)*T^4-2.57479*10^(-8)*T^3+4.63742*10^(-5)*T^2-3.91500*10^(-2)*T+13.2097);
H2_in=16;CO_in=1.124;CO2_in=0.3;CH4_in=0.1;H2O_in=1.656;
n_tot=H2_in+CO_in+CO2_in+CH4_in+H2O_in;
eq1=(H2_in+3*x+y-z)^3*(CO_in-y+x)/((CH4_in-x)*(H2O_in-x-y+z)*n_tot^2)*(Ptot/P0)^2-Kr;
eq2=(H2_in+3*x+y-z)*(CO2_in-y)/((CO_in-y+x)*(H2O_in-x-y+z))-Ks;
eq3=z/(3*x+y)-Uf;
d=solve(eq1,eq2,eq3,'x','y','z');
求解后
警告: The solutions are parametrized by the symbols:
z = RootOf(z1^7 + (132884004*z1^6)/578125 +</