[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Help-gsl] wrong eigenvalues?
From: |
Michael Gross |
Subject: |
[Help-gsl] wrong eigenvalues? |
Date: |
Fri, 07 Nov 2003 00:58:03 +0100 |
User-agent: |
Mozilla/5.0 (X11; U; Linux i686; en-US; rv:1.4) Gecko/20030624 |
Hello
(gsl 1.4, suse 8.2, make check is ok)
When I replace the Matrix in the example on
http://sources.redhat.com/gsl/ref/gsl-ref_14.html#SEC235
with:
double data[] = { 1, 2, 3,
4, 5, 6,
7, 8, 9};
(and replace the 4 in the gsl-calls with 3) then I get the eigenvalues:
18.8302, -3.15747 and -0.672768.
The polynom of the matrix is 18*l+15*l^2-l^3 so one eigenvalue has to be 0.
Matlab (and Maple) calculates (eig([[1 2 3];[4 5 6];[7 8 9]])) the
eigenvalues 16.1168, -1.1168 and -0.0000.
What is wrong?
thx
Michael
- [Help-gsl] wrong eigenvalues?,
Michael Gross <=