function y=korel(f,g)
    z=repmat(g,1,length(f))
    d=zeros(length(f));
    for i = 0:(length(f)-1)
        d(i+1)=sum(f.*z((length(f)^2-length(f)-i+1):(length(f)^2-i)))
    end
    y=d'
endfunction

function y=korelFU(f,g)
    N=length(f)
    y=korel(f,g)/N
    //y=(ifft(fft(f).*conj(fft(g))))/N
endfunction

function y=korelFUN(f,g,t)
    N=length(f)
    k=t
    y=korel(f,g)./(N-abs(k))
    //y=(ifft(fft(f).*conj(fft(g))))./(N-abs(k))
endfunction