function gersh(A)
% plot dei cerchi di Gershgorin di una matrice quadrata A
% necessita la funzione
% http://www.dm.unibo.it/~achilles/calc/programmi/cerchio.m
for k = 1 : length(A)
   c = abs(A(k,k));
   r = sum(abs(A(k,:))) - c;
   cerchio(A(k,k),r);
   hold on;
end
lambda = eig(A);
for k = 1 : length(A)
   plot(real(lambda(k)), imag(lambda(k)),'*');
end;
end
