eig returns nonorthogonal eigenvectors

[chris2]

Hello,
  The following commands show a case where octave returns non-orthogonal
eigenvectors.  Mathematically, the eigenvalue is degenerate, so nothing
prohibits this, even for a symmetric matrix, but it would be very useful if
it did.  Is there any way to force octave to orthogonalize the basis?
--thanks, Chris

octave> [q,r]=qr([1, -1, 1, 0; 2, 3, 0, 1; 4, 3, 0, 0; -1, 2, 0, 0])
octave> a=q*diag([5,5,-7,9])*q'
octave> [va,da]=eig(a)
octave> va'*va




--
View this message in context: 
http://octave.1599824.n4.nabble.com/eig-returns-nonorthogonal-eigenvectors-tp4680711.html
Sent from the Octave - General mailing list archive at Nabble.com.