b = [0,1]; a = [3, -4, 1];                         % polynomials coefficients

[R,p,c] = residuez(b,a)

[b,a] = residuez(R,p,c)