zeroAngle = 15*pi/16;
z0 = exp(j*zeroAngle);
z1 = exp(-j*zeroAngle);
numer = [1 -(z0+z1) z0*z1];

r = 0.9;
poleAngle = pi/16;
p0 = r * exp(j*poleAngle);
p1 = r * exp(-j*poleAngle);
denom = [1 -(p0+p1) p0*p1];

%%% Normalize the DC response to 1 in linear units
%%% by setting H(z) evaluated at z = 1 to be 1
C = (denom * [1 1 1]') / (numer * [1 1 1]');

freqz(C*numer, denom);
