% Testprogramm zu CharMap n = 5 for j=0:n-1 a(j+1) = (j==0)/4 + (3+2*cos(2*pi*j/n))/4/n; %a{j+1}(1,1) = (1+cos(2*pi*j/n))/n; %a{j+1}(1,1) = (j==0)*.98 + (.02 - 0.02*cos(2*pi*j/n))/n; end A = SubMat(a); [lambda, V] = CharMap(A); RingChar(V)