b=[1]
a=[1 1] 
[h,w]=freqs(b,a,15)
plot(w,abs(h)), grid
w'
abs(h)'