function [M,C] = klein(G,D,n1,n2);

[S,T,Q,Z] = qz(G,D);
[S,T,Q,Z] = reorder(S,T,Q,Z);

Z11 = Z(1:n1,1:n1);
Z21 = Z(n1+1:end,1:n1);
Stt = S(1:n1,1:n1);
Ttt = T(1:n1,1:n1);

M = real(Z11*inv(Stt)*Ttt*inv(Z11));
C = real(Z21*inv(Z11));

alf = diag(S) + eps; 
bet = diag(T);
roots = bet./alf;

numb = 0;
for i = 1:length(roots);
     if abs(roots(i,1)) >= 1;
        numb = numb + 1;     
    end;
end;


numb; %numb tells me how many unstable eigenvalues there are%

Z = abs(eig(D,G));
numb1 = sum((Z>=1));