function [Entalpie,RosnyBod] = FceEntalpieSpalin(t,lambda)
%UNTITLED2 Summary of this function goes here
%   Detailed explanation goes here

%lam=1;
lam=lambda;
cs=1.1;         %měrná tepelná kapacita teoretických suchých spalin [kJ/(kg*K)]
ca=1;           %měrná tepelná kapacita suchého vzduchu [kJ/(kg*K)]
a= 1.077;       %součinitel teoretické potřeby vzduchu [-]
%t=60;           %teplota spalin [°C]
x=0.14;         %měrná vlhkost spalin [kg/kg t.s.s.]
p=100;           %atmosferický tlak [kPa]

%parciální tlak páry na mezi sytosti [kPa] - pro teploty 0 - 80 °C 
pp=10^-3*exp(23.58-4044.2./(235.6+t)); 
%[3] Chyský, J., Hemzal,K. a kol.: Technický průvodce "Větrání a klimatizace". Bolit, 1993

%parciální tlak páry v závislosti na měrné vlhkosti spalin
ppx=p/(0.622/(x*(1+a*(lam-1)))+1); %x=0.14

% měrná vlhkost spalin x pro teplotu pod rosným bodem nebude konstantní
xn=0.622*(1+a*(lam-1))*pp/(p-pp);

%teplota rosného bodu 
tr=(-4044.2/(-23.58+log(ppx)-log(0.001)))-235.6; 
%[3] Chyský, J., Hemzal,K. a kol.: Technický průvodce "Větrání a klimatizace". Bolit, 1993
%inverze pp



%Vztah pro entalpii spalin při teplotě nad rosným bodem
h1=(cs + a*(lam-1)*ca)*t+x*(2500+1.84*t);

%Vztah pro entalpii spalin při teplotě rosného bodu
h2=(cs + a*(lam-1)*ca)*tr+x*(2500+1.84*tr);

%Vztah pro entalpii spalin při teplotě pod rosným bodem 
h3=(cs + a*(lam-1)*ca)*t+xn*(2500+1.84*t);

RosnyBod=tr;

%% zobrazení do příkazového řádku
% 
% disp(['Teplota rosného bodu: ' num2str(tr) ' °C'] )
% 
if(t>=tr)
   Entalpie=h1;
   %disp(['Entalpie spalin při teplotě ' num2str(t) ' °C je ' num2str(Entalpie) ' kJ/kg teoretických suchých spalin'])
elseif (t<tr)
   Entalpie=h3;
   %disp(['Entalpie spalin při teplotě ' num2str(t) ' °C je ' num2str(Entalpie) ' kJ/kg teoretických suchých spalin'])
end

end

