function resplot(f,b,m,ifr,ofr,k)
fi=round(ifr/f(2));fo=min(round(ofr/f(2)),length(f));

%---------- plot phase serie --------------------
%
 subplot(2,1,1)
if k==1
	plot(f(fi:fo),angle(b(fi:fo))) 
else  
	plot(f(fi:fo),unwrap(angle(b(fi:fo)),k*pi)) 
end
	ylabel('phase')
 xlabel('frequency [Hz]')
%title(t0s)
%
%------ plot amplitude series--------------------
% 
% b = fft(a);
% m = abs(b(1:length(b)/2));
 subplot(2,1,2)
 plot(f(fi:fo),log10(m(fi:fo)))
 ylabel('amplitude')
 xlabel('frequency [Hz]')
%title(['FFT for ',channel])
