[Help-gsl] Fwd: Eigen vector problem

[belwafi kais]

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.





--