load SymmSnowflake
[V,D]=eig(B);
jj=592:600;
jj=jj(end:-1:1)
for k=1:length(jj)
subplot(3,3,k);
imagesc(D(jj(k),jj(k))*V(:,jj(k))*V(:,jj(k))');
axis equal; axis off;
end
A1=V(:,1:100)*D(1:100,1:100)*V(:,1:100)';
A2=V(:,500:600)*D(500:600,500:600)*V(:,500:600)';
imagesc(A1+A2)
axis equal
axis off