clc
//initialisation of variables
g=10000//K
p=0.01//atm
p1=0.00042 //atm
z=0.201//Btu
Ya=0.666//Btu
Yb=0.167//btu
Yc=0.167 //btu
//CALCULATIONS
Y=Yc//Btu
//RESULTS
printf('the ionization equilibrium constant=% f btu',Y)