exec("thermochim.sci")
data=readData("thermodata.txt");
            
data('SO_2(g)').H
[DH,DS]=reaction(data,"2*SO_2(g)+O_2(g)=2*SO_3(g)")
[DH,DS]
constanteK(DH,DS,700)
plotA=scf();
plotK(DH,DS,300,1500,100)
