function F = myfun(x)
global P;
y=P;
F = [x(8)^2 - y*2*(4.58794*10^(-19))*(0.1057-x(5))/(x(4)*x(2)/2-x(5))-0.1276*10^(-9);
x(8)^2 - y*2*(1.49751*10^(-18))*(0.106-x(6))/(x(4)*(1-x(2))/2-x(6))-0.1276*10^(-9);
x(8)^2 - y*2*(1.01476*10^(-11))*(0.0582-x(7))/(x(4)/2-x(7))-0.1276*10^(-9);
x(5)*x(7)-x(2)*(2.71617*10^(-7))*exp(x(3)/x(8));
x(6)*x(7)-(1-x(2))*(5.71092*10^(-5))*exp(x(3)/x(8));
-x(1)+x(2)*1.24287*10^(-5)+(1-x(2))*1.13096*10^(-5);
-x(3)+2*(x(2)*0.7586+(1-x(2))*0.5718)*x(1)/873.15;
-x(4)+0.709165/(x(2)*1.24287+(1-x(2))*1.13096)];
clc
clear
x = [10^(-5) 0.5 10^(-8) 1 10^(-5) 10^(-5) 10^(-5) 10^(-9)]
x0=x
global P
for i=1:1:10
P=10^(i-1)
options = optimset('TolFun',1e-9);
[x,fval] = fsolve(@myfun,x0,options)
A(:,i)=x
x0=x
end
输出结果也不对,x5~x8输出应为幂的形式