#exemple pour N=3
N=40
A=matrix.toeplitz( [2, -1] + [0]*(N-2), [-1] + [0]*(N-2) )
#A
#A.eigenvectors_right() #la liste des vecteurs propres
#A.eigenvalues() #la liste des valeurs propres
#A*vector([1,0,-1])-2*vector([1,0,-1]) #un test pour N=3 
#A*vector([1,-sqrt(2.),1])-(2.+sqrt(2.))*vector([1,-sqrt(2.),1]) #un autre test pour N=3
#Aeigv=list_plot(A.eigenvalues(),plotjoined=False,rgbcolor=(1.0,0.0,0.0)) #la liste des valeurs propres
#Aeigv #Afficher la liste
A.eigenvectors_right()
list_plot(A.eigenvectors_right()[20][1][0],plotjoined=True,rgbcolor=(1.0,0.0,0.0))#Afficher le 20ième vecteur propre: qui est celui de l'état findamental!