clc //initialisation of variables HfAlO2= -399.1 //kcal HfFe2O2= -196.5 //kcal //CALCULATIONS dH= HfAlO2-HfFe2O2 //RESULTS printf ('Enthalpy change= %.1f kcal mole^-1',dH)