//exec('eqiir.2.ex')
om=[.25*%pi/10,4*%pi/10];
pdB=5;
sdB=120;
deltap=(1.0-10.0**(-0.05*pdB));
deltas=10.00**(-0.05*sdB);
[cells,fact,zers,pols]=eqiir('lp','butt',om,deltap,deltas);
cells
n=prod(cells(2));
d=prod(cells(3));
tr=n./d
