a=[1 2 1];
b=[0 0 1];
[bz,az]=bilinear(b,a,1);
[bz,az]=impinvar(b,a,1);