function Nepr_Disk;
a=[1 -3 4 -2 1];
b=[1 3 3 1];
% [bz,az]=bilinear(b,a,1);
% [bz,az]=impinvar(b,a);