eigenvects function returning wrong answers for complex eigenvectors

[parisse]

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).