function Rxy = jcorr(X, Y)
    N = min(length(X), length(Y));
%    X = X - mean(X);
%    Y = Y - mean(Y);
    for m=0:N-1
        Rxy(m+1,:) = sum(X(1+m:N) .* conj(Y(1:N-m)));
    end