function [ psi,opr,mpsi,os_red ] = zapisnik( sm,psm,psk )
% vypocet zapisniku
% function [ prumery,opravy,std_prumeru ] = zapisnik( smery,poc_smeru,poc_skupin )

% ======================================================================= %
% vypocet -----------------------------------------------------------------
for isk = 1:psk
    bR = 2*(isk-1)*psm+1;   % zacatek I rady v sm
    eR = 2*(isk-1)*psm+psm;   % konec I rady v sm
% prumer z poloh
    for ism = eR+psm:-1:bR+psm
        if( sm(ism,3)<200 )
            sm2pol(ism,3) = sm(ism,3)+200;
        else
            sm2pol(ism,3) = sm(ism,3)-200;
        end
    end
    osnova(:,isk) = (  sm( bR:eR, 3 ) + sm2pol( ism, 3 )  )./2;
% redukce
    os_red(:,isk) = osnova(:,isk)-osnova(1,isk);
end

% prumer ze skupin
psi = mean (os_red, 2);

% smerodatna odchylka radkoveho prumeru
mpsi = std(osnova')/sqrt(psk);
mpsi = mpsi';

% opravy prumeru
opr = repmat(psi,1,psk)-os_red;
