function [Tt2,pt2,m2] = turbina(Tt1,pt1,m1,Pt,et,mch,Tch) %et[et,etR] Cp1=1200; R=287; k1=1/(1-R/Cp1); Cp2=1100; k2=1/(1-R/Cp2); if exist('et','var')==0 et=[1,0]; mch=0; Tch=0; end if exist('mch','var')==0 mch=0; Tch=0; end if exist('Tch','var')==0 Tch=mch; mch=et; et=[1,0]; end Tt2=(m1*Cp1*Tt1-Pt+mch*Cp2*Tch)./(m1*Cp1+mch*Cp2); m2=m1+mch; if et(2)==0 Pit=(et(1)./(et(1)+Tt2./Tt1-1)).^(k1./(k1-1)); elseif et(2)==1 Pit=(Tt1./Tt2)^(k1/(k1-1)/et(1)); end pt2=pt1./Pit; end