function B=comprob_teor_cayley_hamilton_fun(A)
%primera linea de funcion MATLAB !!!

%nombre fichero IGUAL  nombre de funcion

%las variables  son locales


disp('Comprobación Teorema de Cayley-Hamilton')
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

%{
probar con 
A1=ceil(rand(n)*10))
%}


%no olvidar este end
end