function Nayk_Diskr;
a=[1 1];
b=[0 1];
[bz,az]=impinvar(b,a,1);
z=tf(bz,az,1);
nyquist(z);
end
% [h,w]=freqs(b,a);
% plot(w,abs(h)); grid on