%  w = kulmataajuus
%  x = kulmataajuus piin monikertoina eli w=0..pi -> x=0..1
%  A = vahvistus

w=0:0.01:pi;
x=w/pi;
A=abs((exp(j*w*2)-exp(j*w)+0.5)./(exp(j*2*w)-1.3376*exp(j*w)+0.8099));
plot(x,A), grid