clc
//initialisation of variables
hf=-169297//btu/lbm
s1=536.7*(51.072-1.359-49.004)//Btu/lb mole
//CALCULATIONS
G=hf-s1//Btu/lb mole
//RESULTS
printf('The Gibbs function of=% fBtu/lb mole',G)