Sunday, February 14, 2010

interesting Matlab code

N = 256; % Boettcher/Grudsky/Grudskaya matrix
r = [0 -6+13i 5-4i 3 zeros(1,N-9) 3 -4-3i 0];
A = toeplitz([0 fliplr(r(2:end))].',r);
figure, plot(eig(A),'k.'), axis equal

No comments:

Blog Archive