clc h3=75.1134 //in kJ/kg hf=12.5201 //in kJ/kg hg=176.1723 //in kJ/kg X=(h3-hf)/(hg-hf) mprintf("X=%f\n",X)//ans vary due to roundoff error