clc; rp=1.35; // Actual pressure ratio DelT_rise=30; // Actual temperature rise in K beta_1=47; // Inlet blade angle in degree beta_2=15; // Outlet blade angle in degree u=225; // Peripheral velocity in m/s ca=180; // Axial velocity in m/s T01=27+273; // Ambient temperature in kelvin Cp=1.005; // Specific heat in KJ/kg K r=1.4; // Specific heat ratio R=287; // Characteristic gas constant in J/kg K eff_s=(rp^((r-1)/r)-1)*T01/DelT_rise; wf=(DelT_rise*Cp*10^3)/(u*ca*(tand(beta_1)-tand(beta_2))); disp ("%",eff_s*100,"Stage Efficiency = "); disp (wf,"Work done factor = ");