Advanced Search

Question: call to Eigenvectors never returns

Question:call to Eigenvectors never returns

Lawrence Lifshitz 8 Maple
October 27 2006
0
Hi,
I'm trying to use Eigenvectors on a 7x7 matrix (in maple8 on linux). The process
never returns and its size continually grows (>200 MB before I kill the process).
I am new to maple, but have successfully used Eigenvectors on a different (sparser)
7x7 matrix. I'm assuming that either a general 7x7 matrix is too hard for
maple to solve analytically, or something about the degeneracy of my matrix is
causing problems? btw, Nullspace works fine on the matrix. Any help would be
appreciated. I added the assumptions hoping that would simplify the solution,
but it doesn't seem to make any difference. Below is the code.

restart;
with(LinearAlgebra);
Q := Matrix([[-k1on*ca, k1off, 0, 0, 0, 0, 0],
[k1on*ca,-k1off-alpha1-k2on*ca,k2off,0,bf,0,0],
[0,k2on*ca,-k2off-alpha2-k3on*ca,k3off,0,bf,0],
[0,0,k3on*ca,-k3off-alpha3,0,0,bf],
[0,alpha1,0,0,-bf,0,0],
[0,0,alpha2,0,0,-bf,0],
[0,0,0,alpha3,0,0,-bf]]);

assume(k1on>=0,k2on>=0,k3on>=0,k1off>=0,k2off>=0,k3off>=0,ca>=0,alpha1>=0,alpha2>=0,alpha3>=0,bf>=0);
(evals,evecs) := Eigenvectors(Q);
Please Wait...