Advanced Search

Question: bug in Eigenvectors

Question:bug in Eigenvectors

mbras 62 Maple
linearalgebra eigenvectors eigenvalues
March 16 2018
0

I believe I've found a bug, where Eigenvectors returns two eigenvectors which are the same (even though the eigenvalues are different). The expressions involves RootOf's, and it seems that the eigenvalues use "index" to distinguish themselves, yet the eigenvectors do not.

with(LinearAlgebra);

x := RootOf(z^2-t, z);
m := Matrix(2, (i,j) -> evala(add(a[j, k]*((-1)^(i-1))^k*x^k, k = 0..1)));
ev := Eigenvectors(m);
Please Wait...