f1=500;     % notch frequency in Hertz
fs=8000;
h=[1 -2*cos(2*pi*f1/fs) 1]
