%exm0531_2.m

A=[3          -2         -0.9      2*eps
  -2           4         -1       -eps
  -eps/4       eps/2     -1        0
  -0.5        -0.5        0.1      1  ];
[V1,D1]=eig(A);
ER1=A*V1-V1*D1
[V2,D2]=eig(A,'nobalance');
ER2=A*V2-V2*D2