题目:
d(XA)/d(Z)=RA/(3.90469442*0.3884/106)
d(T)/d(Z)=(RA*(-30850)-4*20.02*(T-923.15)/0.123)/(3.90469442*2.22)
RA=C1*0.247*0.95*KE*((PE-PS*PH/KP)+PE*(KT+KB))/(1+LEME*PE+LEMS*PS)^2*1100
PE=34.1952*(1-XA)/(356.452+37.477*XA)*1.529731063
PS=(31.801*XA+0.0208)/(356.452+37.477*XA)*1.529731063
PH=38.982*XA/(356.452+37.477*XA)*1.529731063
KE=5.7e8*exp(-34420/(1.987*T))
KT=1.3e3*exp(-18880/(1.987*T))
KB=9.1e7*exp(-38190/(1.987*T))
KP=exp(16.12-15350/T)
LEME=3.6e-3*exp(-40780/(1.987*T))
LEMS=1.5e-3*exp(-16430/(1.987*T))
C1=((3624.641-61.549*XA)/2265.401*(1-XA)/2.45)^0.25
G=9430.217/(250*3.14159/4*0.123^2)