%Program 6_3
% Partial-Fraction Expansion of Rational z-Transform
% e.g G(z) = 18Z3/(18Z3 + 3Z2 -4Z -1)
% num = input('Type in numerator coefficients = ');
% den = input('Type in denominator coefficients = ');
num = [18 0 0 0];
den = [18 3 -4 -1];

[r,p,k] = residuez(num,den);
disp('Residues');disp(r')
disp('Poles');disp(p')
disp('Constants');disp(k)