b=[1 1]
a=[1.9391 0.0609] 
[h,w]=freqz(b,a,15)
plot(w/2/pi,abs(h)),grid
disp(['w'])
(w/2/pi)'
abs(h)'