function [] = my_gerschgorin_24_03(a) n=length(a); for k=1:n r=0; for j=1:k-1 r=r+abs(a(k,j)); end for j=k+1:n r=r+abs(a(k,j)); end my_circle_24_03(a(k,k),r,1000); hold on; end ee=eig(a); plot(real(ee),imag(ee),'d'); hold off; end