function [H,P,U]=body_zvratu(x,alfa)
% funkce testuje nahodnost posloupnosti pomoci bodu zvratu
% synatxe: [H,P,D]=body_zvratu(x,alfa)
% vystupni parametry:
% H ... vysledek testu> 1 ... H0 zamitame, 0 ... H0 ... nezamitame
% P ... vypocitana p-hodnota
% U ... realizace testove statistiky
% vstupni parametry:
% x ... sloupcovy vektor hodnot testovane posloupnosti
% alfa ... hladina vyznamnosti
n=size(x,1);
Y=0;
for i=2:n-1
    if ((x(i)>x(i-1))&(x(i)>x(i+1)))|((x(i)<x(i-1))&(x(i)<x(i+1)))
        Y=Y+1;
    end
end;
U=(Y-(2*n-4)/3)/sqrt((16*n-29)/90);
P=2*min(normcdf(U,0,1),1-normcdf(U,0,1));
if P<=alfa 
    H=1; 
end
if P>alfa 
    H=0; 
end