
r1 ,p1, k1 = residuez(-4 , [16 , 16 , 4]); % the first term
r2 ,p2, k2 =residuez([5] , [1 , 0 ,  0.64] ) ;% the third term

r = [r1' ,-2/4, r2'  ];
p = [p1' , 6/4 , p2'];
k = k1 + k2;

num , den = residuez(r , p , k)


