grzegorzmazur / yacas

Computer calculations made easy
http://www.yacas.org
GNU Lesser General Public License v2.1
120 stars 23 forks source link

EigenVectors return equations not vectors #320

Open Isomorph70 opened 3 years ago

Isomorph70 commented 3 years ago

Out> {{6,2,-2},{2,6,-2},{-2,-2,6}} In> ata24:=% Out> {{6,2,-2},{2,6,-2},{-2,-2,6}} In> EigenValues(%) Out> {4,4,10} In> e24:=% Out> {4,4,10} In> EigenVectors(ata24,e24) Out> {{k1==k1,k2==(4*k1+2*k3-6*k1)/2,k3==k3},{k1==k1,k2==(4*k1+2*k3-6*k1)/2,k3==k3},{k1== -k3,k2==(10*k1+2*k3-6*k1)/2,k3==k3}}

Verify([Local(a);a:={{6,2,-2},{2,6,-2},{-2,-2,6}};Length(Difference(EigenVectors(a,EigenValues(a)),{{1,0,1},{-1,1,0},{-1,-1,1}}))=0;],True)