function [ndt, R] = corb(a, b, f)
%   Correlation functing using convolution
N = length(a);
R = conv(a, fliplr(b))/(N + 1);
ndt = (-(N - 1): (N - 1))/f;
end

