On Thu, Sep 4, 2008 at 2:16 PM, chris75de <[EMAIL PROTECTED]> wrote: > > Hi, > > the following code > > b=matrix(3,1,[random() for i in [1..3]]) > B=b*b.transpose() > lb,Sb=B.right_eigenvectors() > print lb > print > print Sb.transpose()*B*Sb > print > print Sb.transpose()*Sb > > give's in my sage version 3.0.6 the following result: > > > [1.844328223, 4.46105473802e-17, -3.42287439238e-16] > > [ 1.844328223 -1.11022302463e-16 2.91433543964e-16] > [-8.30521194918e-17 -3.99656191956e-18 2.92386038647e-19] > [ 2.73402117214e-16 2.893962297e-17 -5.58462128744e-17] > > [ 1.0 -1.66533453694e-16 -7.6327832943e-17] > [-1.66533453694e-16 1.0 0.709678542323] > [ -7.6327832943e-17 0.709678542323 1.0] > > I think the last matrix have to be the identity matrix, but there are > two times 0.709678542323 instead of zero.
The identity matrix would mean that Sb is orthogonal, which suggested that the eigenvectors were normalized, doesn't it? > > Christian > > > > > --~--~---------~--~----~------------~-------~--~----~ 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/sage-devel URLs: http://www.sagemath.org -~----------~----~----~----~------~----~------~--~---
