aiman-problem1

pkg load io;
airproperties = xlsread ('Air properties.xlsx');
temperature = airproperties(:,1);
enthalpy  = airproperties(:,2);
pressure = airproperties(:,3);
u = airproperties(:,4);
vr = airproperties(:,5);
s = airproperties(:,6);
temperature1 = 300; temperature3 = 1400; flowrate = 5; p1 = 100; ratio1 = 10; ratio2 = 0.1;
x = 0;
for x = 1:numel(temperature)
    if temperature(x) == temperature1
        break
    end
end
pr1 = pressure(x);
enthalpy1 = enthalpy(x);
pr2 = pr1 * (ratio1);
enthalpy2 = interp1(pressure, enthalpy, pr2);
pr3 = interp1(temperature, pressure, temperature3);
enthalpy3 = interp1(pressure, enthalpy, pr3);
pr4 = pr3 * (ratio2);
enthalpy4 = interp1 (pressure, enthalpy, pr4);
efficiency = ((enthalpy3-enthalpy4)-(enthalpy2-enthalpy1))/(enthalpy3-enthalpy2)
backworkratio = (enthalpy2 - enthalpy1)/(enthalpy3 - enthalpy4)
mass_dot = ((flowrate*p1*(10^3 ))/((8.314/28.97)*temperature1))*(1/10^3);
w_dot_cycle = mass_dot*((enthalpy3 - enthalpy4) - (enthalpy2 - enthalpy1))