clear all

A = [ 7 -4  1  0  0 ;
     -4  6 -4  1  0 ;
      1 -4  6 -4  1 ;
      0  1 -4  6 -4 ;
      0  0  1 -4  5 ] ;
  
[LambdaMax    xmax] = Mises (A) ;
[LambdaMaxInv xmin] = Mises (inv(A)) ;
LambdaMinMax = [1/LambdaMaxInv LambdaMax]
xminmax      = [xmin xmax]
[X Lambda]   = eig (A)