b=fir1(500,[0.49 0.51],'stop')';
[h f] = freqz( b, [1], 512, 8000);
plot( f, abs(h));
xlabel 'Frequency (Hz)'
ylabel 'Magnitude'
title 'Filter Frequency Response'
print -djpeg95 filter.jpeg
