clc h1=2776.4; //kJ/kg h2=h1; h_f1=884.6; //kJ/kg h_fg1=1910.3; //kJ/kg x1=(h1-h_f1)/h_fg1; disp("Initial dryness fraction = ") disp(x1)