function lineseek(f)

load results_before

dn=51.2;fn=f*256/125+1;
fi=round(max(fn-dn,1));fo=round(min(fn+dn,4097));
fr=[fi-1:fo-1]*125/256;

subplot(2,1,1)
tmp=sqrt(varr(fi:fo)+vari(fi:fo));
% Philip, I suppose the means should be taken 
% into account here, it I really want sqrt(power)!
semilogy(fr,tmp);

t=sprintf('(corrected) Amplitude spectra at %6.2f Hz', f);
title(t);
xlabel('Frequency (Hz)');
ylabel('Amplitude before LR');


load results_after

subplot(2,1,2)
tmp=sqrt(varr(fi:fo)+vari(fi:fo));
semilogy(fr,tmp);

xlabel('Frequency (Hz)');
ylabel('Amplitude after LR');