%         11z^2 - 15z + 6      A1         A2          B
% G(z) = ----------------- = ------- + --------- + -------
%        (z - 2) (z - 1)^2   (z - 1)   (z - 1)^2   (z - 2)

num = [11 -15 6];
den = conv([1 -2],conv([1 -1],[1 -1])); 
[res,pol] = respol(num,den) %User defined function

% % res =
% %    20.0000  <------
% %    -9.0000  <-----|--------
% %    -2.0000  <-----|-------|--------
% %                   |       |       |
% % pol =             |       |       |
% %     2.0000  <------       |       |
% %     1.0000  <--------------       |
% %     1.0000  <----------------------