function HPOutput = HPTurbine(m3,T3,P3,T4,P)
Cpg = P(1); % Parameters

P4mid = P3*(T4/T3)^(1/0.23); % Computing intermediate pressure after HP turbine
T4x1 = T3*(P4mid/P3)^0.23; % Computing intermediate temperature after HP turbine
Ws1 = m3*Cpg*(T4x1-T3); % Computing power output from HP turbine

HPOutput = [P4mid; T4x1; Ws1];
end

