>	restart:#"m09_p21"

Los datos termoquímicos tabulados contienen cierta redundancia que se desea comprobar. En particular, se pide:
a) Para cada compuesto tabulado (excluyendo los elementos de referencia), indicar la reacción de formación y calcular las entalpías, funciones de Gibbs y entropías de la reacción de formación en condiciones estandar.
b) Calcular la incertidumbre de los datos tabulados, comparando el residuo de la ecuación grÅ=hrÅ-TsrÅ con el valor de hrÅ, para cada compuesto.

Datos:

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

a) Para cada compuesto tabulado (excluyendo los elementos de referencia), indicar la reacción de formación y calcular las entalpías, funciones de Gibbs y entropías de la reacción de formación en condiciones estandar.
b) Calcular la incertidumbre de los datos tabulados, comparando el residuo de la ecuación grÅ=hrÅ-TsrÅ con el valor de hrÅ, para cada compuesto.

>

for ii from 1 to 29 do line:=readline(`../comb.asc`);if line<>0 then temp:=sscanf(%,`%s %f %f %f`):fi:comp:=convert(substring(temp[1],1..length(temp[1])-3),name);if member(comp,[C,H2,O2,N2,S]) then void; else if member(comp,Comp) then eq||ii:=eq_fit(c1*C+c2*H2+c3*O2+c4*N2+c5*S=comp);r25_:=hgs_r25(%); h25:=temp[2]*1e3*J_/mol_;g25:=temp[3]*1e3*J_/mol_;s25:=temp[4]*J_/(mol_*K_);print(comp,`       h25`=h25,'g25'=g25,'s25'=s25);Err:=100*subs(Const,r25_[1]-T25*r25_[3]-r25_[2])/r25_[1];print(`Eq:`,eq||ii,`   hr25-h25`=r25_[1]-h25,`gr25-g25`=r25_[2]-g25,`sr25`=evalf(r25_[3],3),`Err=`,evalf(Err,1),"%");fi:fi:od:

Se ve que el error es menor del 1%.

>