%---- compute and sort eigenvalues/eigenvectors -----
hold on
A = scAGCT;
[vA, vpA] = eig(A);
[vpA, s] = sort(diag(vpA)); %sort eigenvalues
plot(vpA)
A = scATGC;
[vA, vpA] = eig(A);
[vpA, s] = sort(diag(vpA)); %sort eigenvalues
plot(vpA)
A = scAGCT;
[vA, vpA] = eig(A);
[vpA, s] = sort(diag(vpA)); %sort eigenvalues
plot(vpA)
A = scCGCG;
[vA, vpA] = eig(A);
[vpA, s] = sort(diag(vpA)); %sort eigenvalues
plot(vpA)
A = scAGCT;
[vA, vpA] = eig(A);
[vpA, s] = sort(diag(vpA)); %sort eigenvalues
plot(vpA)
A = scCGGC;
[vA, vpA] = eig(A);
[vpA, s] = sort(diag(vpA)); %sort eigenvalues
plot(vpA)
A = scCGTA;
[vA, vpA] = eig(A);
[vpA, s] = sort(diag(vpA)); %sort eigenvalues
plot(vpA)
A = scTAAT;
[vA, vpA] = eig(A);
[vpA, s] = sort(diag(vpA)); %sort eigenvalues
plot(vpA)
A = scTATA;
[vA, vpA] = eig(A);
[vpA, s] = sort(diag(vpA)); %sort eigenvalues
plot(vpA)
A = scTCGA;
[vA, vpA] = eig(A);
[vpA, s] = sort(diag(vpA)); %sort eigenvalues
plot(vpA)
hold off
line(158*ones(1,160),[1:160])
line(159*ones(1,5),[1:5])
line((158+9)*ones(1,160),[1:160])

%---- same, log scale -----
hold on
A = scAGCT;
[vA, vpA] = eig(A);
[vpA, s] = sort(diag(vpA)); %sort eigenvalues
plot(log(vpA))
A = scATGC;
[vA, vpA] = eig(A);
[vpA, s] = sort(diag(vpA)); %sort eigenvalues
plot(log(vpA))
A = scAGCT;
[vA, vpA] = eig(A);
[vpA, s] = sort(diag(vpA)); %sort eigenvalues
plot(log(vpA))
A = scCGCG;
[vA, vpA] = eig(A);
[vpA, s] = sort(diag(vpA)); %sort eigenvalues
plot(log(vpA))
A = scAGCT;
[vA, vpA] = eig(A);
[vpA, s] = sort(diag(vpA)); %sort eigenvalues
plot(log(vpA))
A = scCGGC;
[vA, vpA] = eig(A);
[vpA, s] = sort(diag(vpA)); %sort eigenvalues
plot(log(vpA))
A = scCGTA;
[vA, vpA] = eig(A);
[vpA, s] = sort(diag(vpA)); %sort eigenvalues
plot(log(vpA))
A = scTAAT;
[vA, vpA] = eig(A);
[vpA, s] = sort(diag(vpA)); %sort eigenvalues
plot(log(vpA))
A = scTATA;
[vA, vpA] = eig(A);
[vpA, s] = sort(diag(vpA)); %sort eigenvalues
plot(log(vpA))
A = scTCGA;
[vA, vpA] = eig(A);
[vpA, s] = sort(diag(vpA)); %sort eigenvalues
plot(log(vpA))
hold off
line(158*ones(1,5),[1:5])
line(159*ones(1,5),[1:5])
line((158+9)*ones(1,5),[1:5])