Hello,
I try to implement the eigenvector of a nonsymetric matrix. To do this I
use gsl_eigen_genv function. And I try to compare the output of this
function with the output of Matbal function ( [c d]=eig(A,B), A and B are
identical matrix, all of these element is one; A=ones(4,4,)). I get two
different results:
Eigenvector
Eigen value
Matlab result
c =
-1 1 0 0
0 -1 0 0
0 0 1 0
0 0 0 1
d =
1 0 0 0
0 0 0 0
0 0 NaN 0
0 0 0 NaN
GSL result
c =
-1 -0 -0 -0
0.7 -1.1 -0.7 0
0 0 1 0
0 0 0 1
d =
1 0 0 0
0 inf 0 0
0 0 -NaN 0
0 0 0 -NaN
Please can you give me any tips to correct this problem?
Thank you.
--
