C121.m

num=[1 0]; den=[1 -0.8];
W=-pi:pi/100:pi;
H=freqz(num,den,W);
mag=abs(H);
phase=180/pi*unwrap(angle(H));
subplot(2,1,1);
plot(W,mag);grid;
subplot(2,1,2);
plot(W,phase);grid