>	restart:#"m16_p04"

>	read`../therm_eq.m`:read`../therm_chem.m`:with(therm_chem);with(therm_proc):

Calcular la concentración de equilibrio de NO en el aire a presión ambiente y a 300 K, 1000 K y 3000 K suponiendo que sólo hay N2, O2 y NO.

>	su1:="Aire":su2:="H2O":dat:=[T1=300*K_,T2=1000*K_,T3=3000*K_];

Eqs. const.:

>	Adat:=get_gas_data(su1):dat:=op(dat),op(subs(g=g0,[Const])),Adat,SI2,SI1:

a) Calcular la concentración de equilibrio de NO en el aire.

Incóg.: a,xNO,xN2,xO2

Ec.: eqNX,eqBO,eqBN,eqEQ

>

eq:=eqMIX(a*c79*N2+a*c21*O2=[2,3,11]);sol:=solve({eqNX,eqBN,eqBO},{a,x[Comp[2]],x[Comp[3]]});eq1:=eqEQ(N2+O2=2*NO);eq1_:=subs(sol,dat,eq1);xNO:=solve(%,x[NO])[2];xN2:=subs(sol,x[NO]=xNO,dat,x[N2]);xO2:=subs(sol,x[NO]=xNO,dat,x[O2]);evalf(seq(subs(T=i*1000*K_,[T_=T,x[NO]=xNO,x[N2]=xN2,x[O2]=xO2]),i=1..3));plot(subs(SI0,[xNO,xN2,xO2]),T=1000..4000,x=0..1);

Pero a alta T pueden aparecer muchos otros.

>