H15 = matkn06[610,15];

vap15 = Eigenvalues[H15];

vep15 = Eigenvectors[H15];

log15[x_,k_]:=If[x==vap15[[k]],0,Log[Abs[x-vap15[[k]]]]];

thouless15[x_]:=Sum[ N[log15[x,k]/610],{k,1,610}];

sth15[n_,m_]:=Sum[(fiboseq[k+1])^2*fiboseq[k]*Exp[-fiboseq[k+1]*(thouless15[vap15[[n]]]+thouless15[vap15[[m]]])],{k,15}];

s15[n_,m_]:= Sum[k^2*vep15[[n,k]]*vep15[[m,k]],{k,610}];
sth15bis[n_,m_]:=(fiboseq[15])^(-2)*Sum[(fiboseq[k+1])^2*fiboseq[k]*vep15[[n,fiboseq[k+1]]]*vep15[[m,fiboseq[k+1]]],{k,14}];

DiscretePlot[s15[1,n],{n,300}]

DiscretePlot[sth15[1,n],{n,300}]

rmsth15[t_]:=Sum[Sin[(vap15[[10*n]]-vap15[[10*m]])*t]*vep15[[10*n,21]]*vep15[[10*m,21]]*sth15[10*n,10*m],{n,1,61},{m,1,61}];

rmsth15[1]

DiscretePlot[{sth15bis[1,n],s15[1,n]},{n,300}]

ListPlot[Table[vep15[[50,n]],{n,610}]]

ListPlot[Table[{fiboseq[k],Abs[vep151[[k]]*vep15200[[k]]]},{k,15}]]

ListPlot[Table[{fiboseq[k],Exp[-fiboseq[k]*(thouless15[vap15[[1]]]+thouless15[vap15[[200]]])]},{k,15}]]

vep151:=Table[vep15[[1,fiboseq[k]]],{k,15}];

vep151

vep15200:=Table[vep15[[200,fiboseq[k]]],{k,15}];

fiboseq[15]