n=20; figure; hold on; for i = 1:n % subplot(n,1,i); v = r.L^(3*i) * [1; zeros(r.n-1,1)]; v = r.L^(3*i) * [zeros(r.n-1,1); 1]; v = v / norm(v); p=plot(v); h = i / n set(p, 'color', [h rand 0.5]); if i == 15 set(p, 'color', 'k'); end pause end