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;