function pseudospectrum(A, epsilon)

% function pseudospectrum(A, epsilon)
%
% vstup: A - ctvercova matice
%        epsilon - parametr
%
% Vykresli aproximaci epsilon-pseudospectra vzorkovaci metodou.

max = 40;                             % pocet vzorku nahodnych matic

for j=1:max,
     E=randn(50) + 1i*randn(50);      % vygeneruj vzorek matice E
     E=epsilon*E/norm(E);
     
     ew=eig(A+E);                     % spocti vl. cisla A+E
     plot(real(ew),imag(ew),'*')     % vykresli
end;

end