//Parameters ;
Parameter x(1:7)[0,14.7];
Constant P=101325,k=1.3806488*10^(-23),e=1.60217733*10^(-19),T=5000,Z1=1836,Z2=4774,Z3=4.041,Z4=8.834,Z5=5.804,Z6=1.000,Ze=7.9354*10^(26),Zn=5.7232*10^(32);
Constant E1=9.756,E2=15.581,E3=14.53413,E4=29.60125,E5=47.4453;
//Variable ;
//Function ;
Function (x1+x2+x3+x4+x5+x6+x7)*10^(23)-P/k/T=0;
x2+x4+2*x5+3*x6-x7=0;
x3*x3*10^(23)-x1*Z3*Z3*Zn/Z1*exp(-E1*e/k/T)=0;
x2*x7*10^(23)-x1*Z2*Ze/Z1*exp(-E2*e/k/T)=0;
x4*x7*10^(23)-x3*Z4*Ze/Z3*exp(-E3*e/k/T)=0;
x5*x7*10^(23)-x4*Z5*Ze/Z4*exp(-E4*e/k/T)=0;
x6*x7*10^(23)-x5*Z6*Ze/Z5*exp(-E5*e/k/T)=0;
ÕâÊÇ´úÂ룬ÓõÄ1.5°æ±¾£¬²»ÖªµÀÊDz»ÊÇÒòΪ°æ±¾µÍµÄÎÊÌ⣬»»Á˺öàÖÖËã·¨£¬×îºóµÄÎó²î¶¼Ìرð´ó
Çó´óÉñ°ïæÅÜһϣ¬¿´Äܲ»Äܵóö±È½ÏÀíÏëµÄ½â£¬¸Ð¼¤²»¾¡~~~ |