disp('Comprobación Teorema de Cayley-Hamilton')
A=input('introduce matriz  cuadrada ')
PA=charpoly(A) %polinomio caracteristico
n=size(A,1)
B=PA(1); %comienzo
for i=1:n
    B=B*A+PA(i+1)*eye(n) %Alg. Horner o Regla Ruffini
end