Example from Hao et al (2001)

Nz = 100  (z^3) + 2  (z^2) + 3  z + 11;
Dz = 
        100  (z^5) + 2  (z^4) + 5  (z^3) - 41  (z^2) + 52  z + 
          70;
Gw=Together[
        ExpandAll[
          (Nz/Dz) /. 
            z -> ((w + 1)/(w - 1))]]

Nw = Numerator[Gw];
Dw = Denominator[Gw];
pcw = Dw + k Nw;
tbl = 
        RouthTabulation[
          CoefficientList[pcw, w]];
cond = Apply[And, Map[(# > 0) &, tbl]]

Reduce[cond // N]