alpha=.9;
f1=2500;
f0=8000;

b=1-alpha;
a=[1 -2*alpha*cos(2*pi*f1/f0) alpha^2];

H=freqz(b,a);
a=a*max(abs(H))

plotZTP(b,a,200,8000)