## Example 15.1     
# P  =   1493.19 MPa
# vG   =   0.000444 m^3/kg
# vD   =   0.000284 m^3/kg
# bTG   =   3.040E-07  / MPa
# bTD   =   1.600E-08  / MPa
# gD   =   2867.8 kJ/kmol
# n  =   12.011 kg/kmol
# error = (vG - vD) *P - (vG * bTG - vD * bTD)P^2 / 2 - gD / n ; kJ/mol 


# Example 15.2     
 T1   =   25 degC
 nH2   =   2.0 
 hfH2   =   0.000 kJ
 sfH2   =   130.678 kJ/degK
 nO2   =   1.0 
 hfO2   =   0.000 kJ
 sfO2   =   205.148 kJ/degK
 nH2O   =   -2.0 
 hfH2O   =   -241826.000 kJ
 sfH2O   =   188.834 kJ/degK
 DH1 = 483652.00 kJ
 DS1 = 88.836 kJ/degK
 DG1 = DH1 - T1 * DS1  ; kJ 
 T2   =   2000 degK
 nH2   =   2.0 
 hfH2   =   52942.000 kJ
 sfH2   =   188.419 kJ/degK
 nO2   =   1.0 
 hfO2   =   59176.000 kJ
 sfO2   =   268.748 kJ/degK
 nH2O   =   -2.0 
 hfH2O   =   -169038.000 kJ
 sfH2O   =   264.769 kJ/degK
 DH2 = 503136.00 kJ
 DS2 = 116.048 kJ/degK
 DG2 = DH2 - T2 * DS2  ; kJ 

# Example 15.3         
 R   =   8.3145 kJ/degK    
 T1   =   25 degC    
 DG1   =   457166 kJ    
 ln_K1 = - DG1 / R / T1
 T2   =   2000 degK    
 DG2   =   271040 kJ    
 ln_K2 = - DG2 / R/ T2  


## Example 15.5         
# K  =   0.1089     
# P_over_P0    =   1.0000     
# z  =   0.2189     
# error = (2* z / (1 - 2* z))^2 * (z / (1 + z)) - K / P_over_P0  
# nCO2   =   0.5622     
# yCO2 = nCO2 / (1 + z) 
# nCO   =   0.4378 
# yCO = nCO / (1 + z)  
# nO2   =   0.2189 
# yO2 = nO2 / (1 + z)  
# hfCO2   =   -393522 kJ
# hf_3000_CO2   =   152853 kJ
# hf_298_CO2   =   0 kJ
# hCO2 = n * (hf + hf_3000 - hf_298) ; kJ 
# hfCO   =   -110527 kJ
# hf_3000_CO   =   93504 kJ
# hf_298_CO   =   0 kJ
# hCO = n * (hf + hf_3000 - hf_298) ; kJ 
# hfO2   =   0 kJ
# hf_3000_O2   =   98013 kJ
# hf_298_O2   =   0 kJ
# hO2 = n * (hf + hf_3000 - hf_298) ; kJ 
# HP = -121301.7 kJ    
# HR   =   0 kJ    
# Q = HP - HR  ; kJ


# Example 15.6         
 K  =   0.1089     
 P_over_P0    =   1.0000     
 z  =   0.1553     
 error = (2* z / (1 - 2* z))^2 * (z / (1 + z)) - K / P_over_P0  
 nCO2   =   0.6894     
 yCO2 = nCO2 / (2 + z)    
 nCO   =   0.3106     
 yCO = nCO / (2 + z)    
 nO2   =   1.1553     
 yO2 = nO2 / (2 + z)  



# Example 15.7         
BEGIN         
 K1   =   0.0020620     
 K2   =   0.0028930     
 P_over_P0    =   1.0000     
 a  =   0.0534     
 b  =   0.0551     
 error1 = ((2*a+b)/(1-2*a-2*b))^2 * (a / (1 + a + b)) - K1 / P_over_P0  
 error2 = (2*a + b) / (1 + a + b) *(2*b / (1 - 2*a - 2*b))^2 - K2 / P_over_P0  
 
 yH2O = (1 - 2*a - 2*b) / (1 + a + b) 
 yH2 = (2*a + b) / (1 + a + b)  
 yO2 = a / (1 + a + b)    
 yOH = 2*b / (1 + a + b)  
END  


# Example 15.8         
BEGIN         
 K  =   0.00042     
 P_over_P0    =   0.0100     
 z  =   0.2008     
 yAr = (1 - z) / (1 + z)  
 yArplus = z / (1 + z)    
 yeminus = z / (1 + z)    
 error = yArplus * yeminus / yAr - K / P_over_P0
END         

# Homework Problems:

TBD