% MATLAB FUNCTION polroots TO COMPUTE THE
% ROOTS OF A SQUARE POLYNOMIAL MATRIX M

% H. Kwakernaak, July, 1990


function rootsM = polroots(M,degM)

[rM,cM] = polsize(M,degM); 
if rM ~= cM
   disp('polroots: Input matrix not square');
end

[d,degd] = poldet(M,degM);
dt = zeros(1,degd+1);
for j = 1:degd+1
    dt(1,j) = d(1,degd+2-j);
end
rootsM = roots(dt)';