function s=irredble(A)
% IRREDBLE   s=irredble(A)
%            Contesta s=1 ´o s=0 si la matriz no negativa nxn A es irreducible o
%            reducible.

% El algoritmo empleado aqui esta  expuesto en el libro de F. Robert:
% "Matrices nonegatives et normes vectorielles", Universite Scientifique et
% Medicale de Grenoble, 1973-1974, pagina 15, Proposition 3.

B=sign(A);
[m,n]=size(B); B=B | eye(n);
for k=2:n
   if potboole(B,k)==potboole(B,k-1)
   P=potboole(B,k-1);
     if P==ones(n)
         s=1;
     else
         s=0;
     end
   break
   end
end

