function freqhist_bfw(nf)
%
load arraydata dd cc nnr nni
%
sr=750000/76;
zr=size(nnr);
freq=(nf-1)*sr/2/(zr(1)-1);
%
xx=[1:zr(2)];
yy=cc(nf)*(2*xx-(zr(2)+1))/(zr(2)-1);
%
% real part of fft
%
subplot(2,1,1)
bar(yy,nnr(nf,:),1)
xlabel('real part fft')
ylabel('number of samples')
%title(['Histograms at ',sprintf('%8.2f',freq),' Hz  from ',num2str(dd),'  blocks'])
%
% imaginary part of fft
%
subplot(2,1,2)
bar(yy,nni(nf,:),1)
xlabel('imaginary part fft')
ylabel('number of samples')