function Rmaxf = demod_Rmaxf (ifreq_N)
N = length(ifreq_N);
%mean_ifreq = (sum( ifreq_N ))/N;
%ifreq_N = ifreq_N - mean_ifreq;  
%ifreq_N = ifreq_N/Fd; 

Y2 = fft(ifreq_N);
YY2 = Y2.* conj(Y2);

Rmaxf = max(YY2) /N;