clear all; clc;

n=[-2 5.65 -2.88];
d=[1 -0.1 0.09 0.648];

xz = tf(n,d,0.1,'variable','z^-1')

[r,p,k] = residuez(n,d)
[bn,an] = residuez([r(1),r(2)],[p(1),p(2)],k);

xz1 = tf(bn,an,0.1,'variable','z^-1')
xz2=tf(r(3),[1 p(3)],0.1,'variable','z^-1')
