clc //initialisation of variables p= 4/3 //atm p1= 1 //atm R= 1.9872 //cal /mole K //CALCULATIONS S= 2*R*log(p/p1) //RESULTS printf (' increase in entropy= %.4f cal deg^-1',S)