Updates:
Status: Fixed
Comment #4 on issue 3343 by [email protected]: bug in eigenvect()
http://code.google.com/p/sympy/issues/detail?id=3343
And it seems to be fixed:
In [1]: from sympy import *
In [2]: M = Matrix([ [1, .6, .6], [.6, 1, .9], [.6, .9, 1] ])
In [3]: print 'M =', M
M = [ 1, 0.6, 0.6]
[0.6, 1, 0.9]
[0.6, 0.9, 1]
In [4]: print 'eig values =',M.eigenvals()
eig values = {-3*sqrt(41)/20 + 29/20: 1, 1/10: 1, 3*sqrt(41)/20 + 29/20: 1}
In [5]: print 'eig vectors =',M.eigenvects()
eig vectors = [(0.100000000000000, 1, [[ 0]
[-1.0]
[ 1.0]]), (0.489531364385073, 1, [[-2.35078105935821]
[ 1.0]
[ 1.0]]), (2.41046863561493, 1, [[0.850781059358212]
[ 1.0]
[ 1.0]])]
--
You received this message because you are subscribed to the Google Groups
"sympy-issues" group.
To post to this group, send email to [email protected].
To unsubscribe from this group, send email to
[email protected].
For more options, visit this group at
http://groups.google.com/group/sympy-issues?hl=en.