function Y=transitiva(T)
% Y=transitiva(T) toma la matriz booleana T (una relacion binaria) y
%                    devuelve su clausura transitiva.

t=length(T);

Y=zeros(t);
for k=1:t
      P=potboole(T,k);
      Y=sumboole(Y,P);
end
