%---- Generalized eig ----
cdir = 'D:\Maths-EPFL-UNI\Travail de master\Matlab';
cd(cdir)
load stiffs
load opt
load filters
load optATGCNN2
A = scATGC;
a = sort(eig(A,A.*M11));
b = sort(eig(A,A.*M53));
c = sort(eig(A,optATGC));
d = sort(eig(A,A.*M97));
e = sort(eig(A,optATGC_97));
gvp = [a';b';c';d';e']';
plot(gvp);
xlabel('Eigenvalue No.')
ylabel('Eigenvalue size')
line([0,210],[1,1],'LineStyle','--','color','black')
legend(['M11  ';'M53  ';'opt53';'M97  ';'opt97'],'Location','NorthWest')
axis([0,210,-1,4.5])