function [Tt2,Pt2,Vh,swl]=wlot(Th,Ph,swl,Ma)
  
%[Tt1,pt1,Vh]=wlot(Th,ph,swl,Ma)
%
%Funkcja oblicza parametry spiêtrzenia strumienia korzystaj¹c z zale¿noœci na 
%sprê¿ izentropowy przy zadanych parametrach statycznych T,p i prêdkoœci Ma 
%Dane na wejœciu
%                    Th - temperatura statyczna
%                    ph – ciœnienie  statyczne
%                    swl – strata ciœnienia we wlocie
%                    Ma – prêdkoœæ lotu wyra¿ona liczb¹ Macha, nie jest wymagana gdy Ma=0,
%Dane na wyjœciu
%                    Tth, pth – temperatura ca³kowita i ciœnienie ca³kowite na wyjœciu z wlotu
%                     Vh prêdkoœæ lotu

if  exist('Ma','var')==0
   Ma=0;
end

if exist('swl','var')==0||isempty(swl)==1
    if Ma<1
        swl=.97;
    else
        swl=.97-0.02241*(Ma-1).^2-0.14561*(Ma-1).^3+0.086282*(Ma-1).^4-0.014342*(Ma-1).^5;
    end
end

Cp=1005 ;
k=1.4;
R=Cp*(1-1/k);  % indywidualna sta³a gazowa obliczana z równania R=cp-cv dla powietrza 287 J/kg/K
Tt2=Th*(1+(k-1)/2*Ma.^2);
Pt2=Ph*swl.*(1+(k-1)/2.*Ma.^2).^(k./(k-1));
Vh=sqrt(k*R*Th)*Ma;
