Hello,
I am trying to determine the limiting behavior of eigenvectors of a matrix and having problem seeing the behavior I expect. The matrix QMQ depends upon k1 and k2 and seems to be a continuous function of k1 and k2. For k1=0 and k2=-1 say it becomes a diagonal matrix and eigen vectors are three coordinate axis vectors. But I can not see the expected limiting behavior ( I know that eigen vectors are scalars invariant). I should be able to show that EA2[2] is equal to (upto a scalar) [1,0,0] in the limiting case (k1=0,k2=-1). For reference other constants are positive and are
a:8.0665;
b:192.6988;
c:3.1721;
d:2.7281;
g:2.4489;
f:4.9619;
Thank you.
Ram
-----------------------------------Code -----------------------------------------------------
load(eigen);
QMQ:matrix([f*k2^2+b*k1^2,f*k1*k2+d*k1*k2,0],
[f*k1*k2+d*k1*k2,a*k2^2+f*k1^2,0],[0,0,g*k2^2+f*k1^2]);
EA2:radcan(eivects(QMQ));
EA2[2];