function gershdisc(A) %% Spočte a vykreslí přesná vlastní čísla a Geršgorinovi kruhy pro matici A. R = ... % poloměry jsou sumy mimodiag. prvků % Kreslíme přes řádky A jednotlivé kruhy: B=256; % zvolíme body pro vykreslení kruhu t=(0:B)*2*pi/B; N = size(A,1); for i=1:N h= ... % souřadnice středu je diag. prvek k= ... plot(R(i)*cos(t)+h,R(i)*sin(t)+k,'-'); % vykreslení Gerš. kruhu hold on; plot(h, k,'+'); % střed Gerš. kruhu end axis equal; % čtvercové osy pro vizualizaci % Výpočet přesných vl.č.: ... end