eigenvects function returning wrong answers for complex eigenvectors

[medwatt]

(01-02-2020 09:05 AM)parisse Wrote:  The reason is that eigVl and eigenvects do not necessarily sort the eigenvectors the same way. You should use jordan to compute both simultaneously:
a:=[[0,0.1],[-1,-0.2]];
p,d:=jordan(a);
p*d*p^-1;
In addition, this is more efficient (for large matrices).

I'm actually not interested in using the diagonalization formula. What I'm interested in is the result from the eigenvects function. Why doesn't the function return a simple answer?

[[-0.1-0.3*i, -0.1+0.3*i],[1,1]] - this is the simple answer I expect. Instead, what I got is something that doesn't even tell me that the eigenvectors are complex conjugates. I've tried it on Wolfram alpha and lots of other websites and all of them return [[-0.1-0.3*i, -0.1+0.3*i],[1,1]].