Eigenvectors and Eigenvalues
Worksheet by Russell Blyth and Mike May, S.J.
restart:with(plots):with(LinearAlgebra):
Warning, the name changecoords has been redefined
Using Maple commands to compute eigenvectors
We first look at the Maple commands for computing eigenvalues and eigenvectors. These commands let us work with eigenvectors and eigenvalues before we work through the detailed computations.
We start with a randomly generated 2 by 2 matrix.
A := RandomMatrix(2,2,generator=rand(-5..5));
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5USJBRigvJSdmYW1pbHlHUTBUaW1lc35OZXd+Um9tYW5GKC8lJXNpemVHUSMxMkYoLyUlYm9sZEdRJmZhbHNlRigvJSdpdGFsaWNHUSV0cnVlRigvJSp1bmRlcmxpbmVHRjgvJSpzdWJzY3JpcHRHRjgvJSxzdXBlcnNjcmlwdEdGOC8lK2ZvcmVncm91bmRHUSpbMCwwLDI1NV1GKC8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdRigvJSdvcGFxdWVHRjgvJStleGVjdXRhYmxlR0Y4LyUpcmVhZG9ubHlHRjsvJSljb21wb3NlZEdGOC8lKmNvbnZlcnRlZEdGOC8lK2ltc2VsZWN0ZWRHRjgvJSxwbGFjZWhvbGRlckdGOC8lMGZvbnRfc3R5bGVfbmFtZUdRKjJEfk91dHB1dEYoLyUqbWF0aGNvbG9yR0ZELyUvbWF0aGJhY2tncm91bmRHRkcvJStmb250ZmFtaWx5R0YyLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGKC8lKW1hdGhzaXplR0Y1LUkjbW9HRiU2M1EjOj1GKC8lJWZvcm1HUSZpbmZpeEYoLyUmZmVuY2VHRjgvJSpzZXBhcmF0b3JHRjgvJSdsc3BhY2VHUS90aGlja21hdGhzcGFjZUYoLyUncnNwYWNlR0ZbcC8lKXN0cmV0Y2h5R0Y4LyUqc3ltbWV0cmljR0Y4LyUobWF4c2l6ZUdRKWluZmluaXR5RigvJShtaW5zaXplR1EiMUYoLyUobGFyZ2VvcEdGOC8lLm1vdmFibGVsaW1pdHNHRjgvJSdhY2NlbnRHRjgvJTBmb250X3N0eWxlX25hbWVHRlgvJSVzaXplR0Y1LyUrZm9yZWdyb3VuZEdGRC8lK2JhY2tncm91bmRHRkctRiQ2JS1GX282M1EiW0YoL0Zjb1EncHJlZml4RigvRmZvRjtGZ28vRmpvUS50aGlubWF0aHNwYWNlRigvRl1wRl9yL0ZfcEY7RmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiMtSSdtdGFibGVHRiU2JC1JJG10ckdGJTYkLUkkbXRkR0YlNiMtSSNtbkdGJTY5USkmbWludXM7MkYoRjBGM0Y2L0Y6RjhGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduL0ZqblEnbm9ybWFsRihGXG8tRltzNiMtRl5zNjlRIjJGKEYwRjNGNkZhc0Y8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GYnNGXG8tRmhyNiQtRltzNiMtRl5zNjlRIjVGKEYwRjNGNkZhc0Y8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GYnNGXG8tRltzNiMtRl5zNjlRKSZtaW51czszRihGMEYzRjZGYXNGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmJzRlxvLUZfbzYzUSJdRigvRmNvUShwb3N0Zml4RihGXXJGZ29GXnIvRl1wUTJ2ZXJ5dGhpbm1hdGhzcGFjZUYoRmFyRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxNyMtX0YpSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSJBR0YoLUknUlRBQkxFR0YoNiUiKCFHIz4mLUknTUFUUklYR0YoNiM3JDckISIjIiIjNyQiIiYhIiRJJ01hdHJpeEc2JCUqcHJvdGVjdGVkR0YqNyMtRmJ2NiMvSSQlaWRHRihGZ3U=
Eigenvectors(A);
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtRiQ2JS1JI21vR0YlNjNRIltGKC8lJWZvcm1HUSdwcmVmaXhGKC8lJmZlbmNlR1EldHJ1ZUYoLyUqc2VwYXJhdG9yR1EmZmFsc2VGKC8lJ2xzcGFjZUdRLnRoaW5tYXRoc3BhY2VGKC8lJ3JzcGFjZUdGPS8lKXN0cmV0Y2h5R0Y3LyUqc3ltbWV0cmljR0Y6LyUobWF4c2l6ZUdRKWluZmluaXR5RigvJShtaW5zaXplR1EiMUYoLyUobGFyZ2VvcEdGOi8lLm1vdmFibGVsaW1pdHNHRjovJSdhY2NlbnRHRjovJTBmb250X3N0eWxlX25hbWVHUSoyRH5PdXRwdXRGKC8lJXNpemVHUSMxMkYoLyUrZm9yZWdyb3VuZEdRKlswLDAsMjU1XUYoLyUrYmFja2dyb3VuZEdRLlsyNTUsMjU1LDI1NV1GKC1GJDYjLUknbXRhYmxlR0YlNiQtSSRtdHJHRiU2Iy1JJG10ZEdGJTYjLUYkNiYtRi82M1EqJnVtaW51czA7RihGMi9GNkY6RjgvRjxRJDBlbUYoL0Y/RmhvL0ZBRjpGQkZERkdGSkZMRk5GUEZTRlZGWS1JJm1mcmFjR0YlNiotSSNtbkdGJTY5USI1RigvJSdmYW1pbHlHUTBUaW1lc35OZXd+Um9tYW5GKC8lJXNpemVHRlUvJSVib2xkR0Y6LyUnaXRhbGljR0Y6LyUqdW5kZXJsaW5lR0Y6LyUqc3Vic2NyaXB0R0Y6LyUsc3VwZXJzY3JpcHRHRjovJStmb3JlZ3JvdW5kR0ZYLyUrYmFja2dyb3VuZEdGZW4vJSdvcGFxdWVHRjovJStleGVjdXRhYmxlR0Y6LyUpcmVhZG9ubHlHRjcvJSljb21wb3NlZEdGOi8lKmNvbnZlcnRlZEdGOi8lK2ltc2VsZWN0ZWRHRjovJSxwbGFjZWhvbGRlckdGOi8lMGZvbnRfc3R5bGVfbmFtZUdGUi8lKm1hdGhjb2xvckdGWC8lL21hdGhiYWNrZ3JvdW5kR0Zlbi8lK2ZvbnRmYW1pbHlHRmRwLyUsbWF0aHZhcmlhbnRHUSdub3JtYWxGKC8lKW1hdGhzaXplR0ZVLUZfcDY5USIyRihGYnBGZXBGZ3BGaXBGW3FGXXFGX3FGYXFGY3FGZXFGZ3FGaXFGW3JGXXJGX3JGYXJGY3JGZXJGZ3JGaXJGW3NGXnMvJS5saW5ldGhpY2tuZXNzR1EiMUYoLyUrZGVub21hbGlnbkdRJ2NlbnRlckYoLyUpbnVtYWxpZ25HRmhzLyUpYmV2ZWxsZWRHRjpGVkZZLUYvNjNRIitGKC9GM1EmaW5maXhGKEZmb0Y4L0Y8UTBtZWRpdW1tYXRoc3BhY2VGKC9GP0ZjdEZqb0ZCRkRGR0ZKRkxGTkZQRlNGVkZZLUYkNiUtRlxwNiotRl9wNjlGSUZicEZlcEZncEZpcEZbcUZdcUZfcUZhcUZjcUZlcUZncUZpcUZbckZdckZfckZhckZjckZlckZnckZpckZbc0Zec0Zgc0Zjc0Zmc0Zpc0ZbdEZWRlktRi82M1ExJkludmlzaWJsZVRpbWVzO0YoRmB0RmZvRjhGZ29GaW9Gam9GQkZERkdGSkZMRk5GUEZTRlZGWS1GJDYjLUkmbXNxcnRHRiU2JS1GX3A2OVEjNDFGKEZicEZlcEZncEZpcEZbcUZdcUZfcUZhcUZjcUZlcUZncUZpcUZbckZdckZfckZhckZjckZlckZnckZpckZbc0Zec0ZWRlktRlxvNiMtRl9vNiMtRiQ2JkZjb0ZbcC1GLzYzUSgmbWludXM7RihGYHRGZm9GOEZidEZkdEZqb0ZCRkRGR0ZKRkxGTkZQRlNGVkZZRmV0LUYvNjNRIl1GKC9GM1EocG9zdGZpeEYoRjVGOEY7L0Y/UTJ2ZXJ5dGhpbm1hdGhzcGFjZUYoRkBGQkZERkdGSkZMRk5GUEZTRlZGWS1GLzYzUSIsRihGYHRGZm8vRjlGN0Znby9GP1EzdmVyeXRoaWNrbWF0aHNwYWNlRihGam9GQkZERkdGSkZMRk5GUEZTRlZGWS1GJDYlRi4tRiQ2Iy1GaW42JC1GXG82JC1GX282Iy1GXHA2Ki1GJDYjRmBzLUYkNiMtRiQ2JkZjb0ZndEZddEZldEZjc0Zmc0Zpc0ZbdEZWRlktRl9vNiMtRlxwNipGaHctRiQ2Iy1GJDYmRmNvRmd0Rlx2RmV0RmNzRmZzRmlzRlt0RlZGWS1GXG82JC1GX282I0ZpdEZoeEZfdjcjLV9GKUksbXByaW50c2xhc2hHRig2JDckLUknUlRBQkxFR0YoNiUiKDdwUCMtSSdNQVRSSVhHRig2IzckNyMsJiMiIiYiIiMhIiIqJiMiIiJGXHpGYHopIiNURl96RmB6RmB6NyMsJkZqeUZdekZeekZdeiZJJ1ZlY3Rvckc2JCUqcHJvdGVjdGVkR0YqNiNJJ2NvbHVtbkdGKC1GYXk2JSIoU283KC1GZXk2IzckNyQsJComRlx6RmB6LCZGX3pGXXpGXnpGYHpGXXpGYHosJComRlx6RmB6LCZGX3pGXXpGXnpGXXpGXXpGYHo3JEZgekZgekknTWF0cml4R0ZnejckLUZlejYjL0kkJWlkR0YoRmN5LUZpW2w2Iy9GXlxsRl1bbA==
The output is an ordered pair, where the first entry is a vector of the Eigenvalues, and the second entry is a matrix whose columns are the corresponding eigenvectors.
If we specify that the output should be a list, we get a set of triples, where the first entry of each triple is an eigenvalue, the second entry is its multiplicity as a root of the characteristic polynomial, and the third entry is a basis for the eigenspace corresponding to that eigenvector. We can extract this data if we first place these two triples into a sequence.
evA := Eigenvectors(A, output='list');
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5USRldkFGKC8lJ2ZhbWlseUdRMFRpbWVzfk5ld35Sb21hbkYoLyUlc2l6ZUdRIzEyRigvJSVib2xkR1EmZmFsc2VGKC8lJ2l0YWxpY0dRJXRydWVGKC8lKnVuZGVybGluZUdGOC8lKnN1YnNjcmlwdEdGOC8lLHN1cGVyc2NyaXB0R0Y4LyUrZm9yZWdyb3VuZEdRKlswLDAsMjU1XUYoLyUrYmFja2dyb3VuZEdRLlsyNTUsMjU1LDI1NV1GKC8lJ29wYXF1ZUdGOC8lK2V4ZWN1dGFibGVHRjgvJSlyZWFkb25seUdGOy8lKWNvbXBvc2VkR0Y4LyUqY29udmVydGVkR0Y4LyUraW1zZWxlY3RlZEdGOC8lLHBsYWNlaG9sZGVyR0Y4LyUwZm9udF9zdHlsZV9uYW1lR1EqMkR+T3V0cHV0RigvJSptYXRoY29sb3JHRkQvJS9tYXRoYmFja2dyb3VuZEdGRy8lK2ZvbnRmYW1pbHlHRjIvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YoLyUpbWF0aHNpemVHRjUtSSNtb0dGJTYzUSM6PUYoLyUlZm9ybUdRJmluZml4RigvJSZmZW5jZUdGOC8lKnNlcGFyYXRvckdGOC8lJ2xzcGFjZUdRL3RoaWNrbWF0aHNwYWNlRigvJSdyc3BhY2VHRltwLyUpc3RyZXRjaHlHRjgvJSpzeW1tZXRyaWNHRjgvJShtYXhzaXplR1EpaW5maW5pdHlGKC8lKG1pbnNpemVHUSIxRigvJShsYXJnZW9wR0Y4LyUubW92YWJsZWxpbWl0c0dGOC8lJ2FjY2VudEdGOC8lMGZvbnRfc3R5bGVfbmFtZUdGWC8lJXNpemVHRjUvJStmb3JlZ3JvdW5kR0ZELyUrYmFja2dyb3VuZEdGRy1GJDYlLUZfbzYzUSJbRigvRmNvUSdwcmVmaXhGKC9GZm9GO0Znby9Gam9RLnRoaW5tYXRoc3BhY2VGKC9GXXBGX3IvRl9wRjtGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtRiQ2JS1GJDYlRmhxLUYkNictRiQ2Ji1GX282M1EqJnVtaW51czA7RihGW3JGZW9GZ28vRmpvUSQwZW1GKC9GXXBGXnNGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtSSZtZnJhY0dGJTYqLUkjbW5HRiU2OVEiNUYoRjBGM0Y2L0Y6RjhGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduL0ZqblEnbm9ybWFsRihGXG8tRmRzNjlRIjJGKEYwRjNGNkZnc0Y8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaHNGXG8vJS5saW5ldGhpY2tuZXNzR1EiMUYoLyUrZGVub21hbGlnbkdRJ2NlbnRlckYoLyUpbnVtYWxpZ25HRmJ0LyUpYmV2ZWxsZWRHRjhGYnFGZHEtRl9vNjNRIitGKEZib0Zlb0Znby9Gam9RMG1lZGl1bW1hdGhzcGFjZUYoL0ZdcEZbdUZecEZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1GJDYlLUZhczYqLUZkczY5RmdwRjBGM0Y2RmdzRjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZoc0Zcb0Zqc0ZddEZgdEZjdEZldEZicUZkcS1GX282M1ExJkludmlzaWJsZVRpbWVzO0YoRmJvRmVvRmdvRl1zRl9zRl5wRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiMtSSZtc3FydEdGJTYlLUZkczY5USM0MUYoRjBGM0Y2RmdzRjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZoc0Zcb0ZicUZkcS1GX282M1EiLEYoRmJvRmVvL0Zob0Y7Rl1zL0ZdcFEzdmVyeXRoaWNrbWF0aHNwYWNlRihGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHFGYXVGXnYtRiQ2JS1GX282M1EifGZyRihGW3JGXXJGZ29GXnJGYHJGYXJGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtRiQ2Iy1GJDYlRmhxLUYkNiMtSSdtdGFibGVHRiU2JC1JJG10ckdGJTYjLUkkbXRkR0YlNiMtRmFzNiotRiQ2I0Zqcy1GJDYjLUYkNiZGanJGX3VGZ3RGXXVGXXRGYHRGY3RGZXRGYnFGZHEtRmN3NiMtRmZ3NiNGYXUtRl9vNjNRIl1GKC9GY29RKHBvc3RmaXhGKEZdckZnb0Zeci9GXXBRMnZlcnl0aGlubWF0aHNwYWNlRihGYXJGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtRl9vNjNRInxockYoRmd4Rl1yRmdvRl5yRml4RmFyRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxRmR4Rl52LUYkNiVGaHEtRiQ2Jy1GJDYmRmpyRmBzLUZfbzYzUSgmbWludXM7RihGYm9GZW9GZ29GanRGXHVGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHFGXXVGXnZGYXVGXnYtRiQ2JUZmdi1GJDYjLUYkNiVGaHEtRiQ2Iy1GYHc2JC1GY3c2Iy1GZnc2Iy1GYXM2KkZqdy1GJDYjLUYkNiZGanJGX3VGZHlGXXVGXXRGYHRGY3RGZXRGYnFGZHFGYHhGZHhGW3lGZHhGZHg3Iy1fRilJLG1wcmludHNsYXNoR0YoNiQ3Iz5JJGV2QUdGKDckNyUsJiMiIiYiIiMhIiIqJiMiIiJGaFtsRlxcbCkiI1RGW1xsRlxcbEZcXGxGXFxsPCMtSSdSVEFCTEVHRig2JSIoMyVwQy1JJ01BVFJJWEdGKDYjNyQ3IywkKiZGaFtsRlxcbCwmRltcbEZpW2xGaltsRlxcbEZpW2xGXFxsNyNGXFxsJkknVmVjdG9yRzYkJSpwcm90ZWN0ZWRHRio2I0knY29sdW1uR0YoNyUsJkZmW2xGaVtsRmpbbEZpW2xGXFxsPCMtRmFcbDYlIignSDIjKS1GZVxsNiM3JDcjLCQqJkZoW2xGXFxsLCZGW1xsRmlbbEZqW2xGaVtsRmlbbEZcXGxGXF1sRl1dbDcjNyQ3JUZlW2xGXFxsPCMtRl1dbDYjL0kkJWlkR0YoRmNcbDclRmRdbEZcXGw8Iy1GXV1sNiMvRmdebEZoXWw=
The eigenvalues:
evA[1][1]; evA[2][1];
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiYtSSNtb0dGJTYzUSomdW1pbnVzMDtGKC8lJWZvcm1HUSdwcmVmaXhGKC8lJmZlbmNlR1EmZmFsc2VGKC8lKnNlcGFyYXRvckdGNS8lJ2xzcGFjZUdRJDBlbUYoLyUncnNwYWNlR0Y6LyUpc3RyZXRjaHlHRjUvJSpzeW1tZXRyaWNHRjUvJShtYXhzaXplR1EpaW5maW5pdHlGKC8lKG1pbnNpemVHUSIxRigvJShsYXJnZW9wR0Y1LyUubW92YWJsZWxpbWl0c0dGNS8lJ2FjY2VudEdGNS8lMGZvbnRfc3R5bGVfbmFtZUdRKjJEfk91dHB1dEYoLyUlc2l6ZUdRIzEyRigvJStmb3JlZ3JvdW5kR1EqWzAsMCwyNTVdRigvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XUYoLUkmbWZyYWNHRiU2Ki1JI21uR0YlNjlRIjVGKC8lJ2ZhbWlseUdRMFRpbWVzfk5ld35Sb21hbkYoLyUlc2l6ZUdGUi8lJWJvbGRHRjUvJSdpdGFsaWNHRjUvJSp1bmRlcmxpbmVHRjUvJSpzdWJzY3JpcHRHRjUvJSxzdXBlcnNjcmlwdEdGNS8lK2ZvcmVncm91bmRHRlUvJStiYWNrZ3JvdW5kR0ZYLyUnb3BhcXVlR0Y1LyUrZXhlY3V0YWJsZUdGNS8lKXJlYWRvbmx5R1EldHJ1ZUYoLyUpY29tcG9zZWRHRjUvJSpjb252ZXJ0ZWRHRjUvJStpbXNlbGVjdGVkR0Y1LyUscGxhY2Vob2xkZXJHRjUvJTBmb250X3N0eWxlX25hbWVHRk8vJSptYXRoY29sb3JHRlUvJS9tYXRoYmFja2dyb3VuZEdGWC8lK2ZvbnRmYW1pbHlHRlxvLyUsbWF0aHZhcmlhbnRHUSdub3JtYWxGKC8lKW1hdGhzaXplR0ZSLUZnbjY5USIyRihGam5GXW9GX29GYW9GY29GZW9GZ29GaW9GW3BGXXBGX3BGYXBGZHBGZnBGaHBGanBGXHFGXnFGYHFGYnFGZHFGZ3EvJS5saW5ldGhpY2tuZXNzR1EiMUYoLyUrZGVub21hbGlnbkdRJ2NlbnRlckYoLyUpbnVtYWxpZ25HRmFyLyUpYmV2ZWxsZWRHRjVGU0ZWLUYtNjNRIitGKC9GMVEmaW5maXhGKEYzRjYvRjlRMG1lZGl1bW1hdGhzcGFjZUYoL0Y8RlxzRj1GP0ZBRkRGR0ZJRktGTUZQRlNGVi1GJDYlLUZaNiotRmduNjlGRkZqbkZdb0Zfb0Zhb0Zjb0Zlb0Znb0Zpb0ZbcEZdcEZfcEZhcEZkcEZmcEZocEZqcEZccUZecUZgcUZicUZkcUZncUZpcUZcckZfckZickZkckZTRlYtRi02M1ExJkludmlzaWJsZVRpbWVzO0YoRmlyRjNGNkY4RjtGPUY/RkFGREZHRklGS0ZNRlBGU0ZWLUYkNiMtSSZtc3FydEdGJTYlLUZnbjY5USM0MUYoRmpuRl1vRl9vRmFvRmNvRmVvRmdvRmlvRltwRl1wRl9wRmFwRmRwRmZwRmhwRmpwRlxxRl5xRmBxRmJxRmRxRmdxRlNGVjcjLCYjIiImIiIjISIiKiYjIiIiRmN0Rmd0KSIjVEZmdEZndEZndA==
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiYtSSNtb0dGJTYzUSomdW1pbnVzMDtGKC8lJWZvcm1HUSdwcmVmaXhGKC8lJmZlbmNlR1EmZmFsc2VGKC8lKnNlcGFyYXRvckdGNS8lJ2xzcGFjZUdRJDBlbUYoLyUncnNwYWNlR0Y6LyUpc3RyZXRjaHlHRjUvJSpzeW1tZXRyaWNHRjUvJShtYXhzaXplR1EpaW5maW5pdHlGKC8lKG1pbnNpemVHUSIxRigvJShsYXJnZW9wR0Y1LyUubW92YWJsZWxpbWl0c0dGNS8lJ2FjY2VudEdGNS8lMGZvbnRfc3R5bGVfbmFtZUdRKjJEfk91dHB1dEYoLyUlc2l6ZUdRIzEyRigvJStmb3JlZ3JvdW5kR1EqWzAsMCwyNTVdRigvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XUYoLUkmbWZyYWNHRiU2Ki1JI21uR0YlNjlRIjVGKC8lJ2ZhbWlseUdRMFRpbWVzfk5ld35Sb21hbkYoLyUlc2l6ZUdGUi8lJWJvbGRHRjUvJSdpdGFsaWNHRjUvJSp1bmRlcmxpbmVHRjUvJSpzdWJzY3JpcHRHRjUvJSxzdXBlcnNjcmlwdEdGNS8lK2ZvcmVncm91bmRHRlUvJStiYWNrZ3JvdW5kR0ZYLyUnb3BhcXVlR0Y1LyUrZXhlY3V0YWJsZUdGNS8lKXJlYWRvbmx5R1EldHJ1ZUYoLyUpY29tcG9zZWRHRjUvJSpjb252ZXJ0ZWRHRjUvJStpbXNlbGVjdGVkR0Y1LyUscGxhY2Vob2xkZXJHRjUvJTBmb250X3N0eWxlX25hbWVHRk8vJSptYXRoY29sb3JHRlUvJS9tYXRoYmFja2dyb3VuZEdGWC8lK2ZvbnRmYW1pbHlHRlxvLyUsbWF0aHZhcmlhbnRHUSdub3JtYWxGKC8lKW1hdGhzaXplR0ZSLUZnbjY5USIyRihGam5GXW9GX29GYW9GY29GZW9GZ29GaW9GW3BGXXBGX3BGYXBGZHBGZnBGaHBGanBGXHFGXnFGYHFGYnFGZHFGZ3EvJS5saW5ldGhpY2tuZXNzR1EiMUYoLyUrZGVub21hbGlnbkdRJ2NlbnRlckYoLyUpbnVtYWxpZ25HRmFyLyUpYmV2ZWxsZWRHRjVGU0ZWLUYtNjNRKCZtaW51cztGKC9GMVEmaW5maXhGKEYzRjYvRjlRMG1lZGl1bW1hdGhzcGFjZUYoL0Y8RlxzRj1GP0ZBRkRGR0ZJRktGTUZQRlNGVi1GJDYlLUZaNiotRmduNjlGRkZqbkZdb0Zfb0Zhb0Zjb0Zlb0Znb0Zpb0ZbcEZdcEZfcEZhcEZkcEZmcEZocEZqcEZccUZecUZgcUZicUZkcUZncUZpcUZcckZfckZickZkckZTRlYtRi02M1ExJkludmlzaWJsZVRpbWVzO0YoRmlyRjNGNkY4RjtGPUY/RkFGREZHRklGS0ZNRlBGU0ZWLUYkNiMtSSZtc3FydEdGJTYlLUZnbjY5USM0MUYoRmpuRl1vRl9vRmFvRmNvRmVvRmdvRmlvRltwRl1wRl9wRmFwRmRwRmZwRmhwRmpwRlxxRl5xRmBxRmJxRmRxRmdxRlNGVjcjLCYjIiImIiIjISIiKiYjIiIiRmN0Rmd0KSIjVEZmdEZndEZkdA==
The eigenvectors:
[evA[1][3][1], evA[2][3][1]];
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtb0dGJTYzUSJbRigvJSVmb3JtR1EncHJlZml4RigvJSZmZW5jZUdRJXRydWVGKC8lKnNlcGFyYXRvckdRJmZhbHNlRigvJSdsc3BhY2VHUS50aGlubWF0aHNwYWNlRigvJSdyc3BhY2VHRjsvJSlzdHJldGNoeUdGNS8lKnN5bW1ldHJpY0dGOC8lKG1heHNpemVHUSlpbmZpbml0eUYoLyUobWluc2l6ZUdRIjFGKC8lKGxhcmdlb3BHRjgvJS5tb3ZhYmxlbGltaXRzR0Y4LyUnYWNjZW50R0Y4LyUwZm9udF9zdHlsZV9uYW1lR1EqMkR+T3V0cHV0RigvJSVzaXplR1EjMTJGKC8lK2ZvcmVncm91bmRHUSpbMCwwLDI1NV1GKC8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdRigtRiQ2JS1GJDYlRiwtRiQ2Iy1JJ210YWJsZUdGJTYkLUkkbXRyR0YlNiMtSSRtdGRHRiU2Iy1JJm1mcmFjR0YlNiotRiQ2Iy1JI21uR0YlNjlRIjJGKC8lJ2ZhbWlseUdRMFRpbWVzfk5ld35Sb21hbkYoLyUlc2l6ZUdGUy8lJWJvbGRHRjgvJSdpdGFsaWNHRjgvJSp1bmRlcmxpbmVHRjgvJSpzdWJzY3JpcHRHRjgvJSxzdXBlcnNjcmlwdEdGOC8lK2ZvcmVncm91bmRHRlYvJStiYWNrZ3JvdW5kR0ZZLyUnb3BhcXVlR0Y4LyUrZXhlY3V0YWJsZUdGOC8lKXJlYWRvbmx5R0Y1LyUpY29tcG9zZWRHRjgvJSpjb252ZXJ0ZWRHRjgvJStpbXNlbGVjdGVkR0Y4LyUscGxhY2Vob2xkZXJHRjgvJTBmb250X3N0eWxlX25hbWVHRlAvJSptYXRoY29sb3JHRlYvJS9tYXRoYmFja2dyb3VuZEdGWS8lK2ZvbnRmYW1pbHlHRl5wLyUsbWF0aHZhcmlhbnRHUSdub3JtYWxGKC8lKW1hdGhzaXplR0ZTLUYkNiMtRiQ2Ji1GLTYzUSomdW1pbnVzMDtGKEYwL0Y0RjhGNi9GOlEkMGVtRigvRj1GY3MvRj9GOEZARkJGRUZIRkpGTEZORlFGVEZXLUZkbzYqLUZpbzY5RkdGXHBGX3BGYXBGY3BGZXBGZ3BGaXBGW3FGXXFGX3FGYXFGY3FGZXFGZ3FGaXFGW3JGXXJGX3JGYXJGY3JGZXJGaHJGaG8vJS5saW5ldGhpY2tuZXNzR1EiMUYoLyUrZGVub21hbGlnbkdRJ2NlbnRlckYoLyUpbnVtYWxpZ25HRl90LyUpYmV2ZWxsZWRHRjhGVEZXLUYtNjNRIitGKC9GMVEmaW5maXhGKEZhc0Y2L0Y6UTBtZWRpdW1tYXRoc3BhY2VGKC9GPUZqdEZlc0ZARkJGRUZIRkpGTEZORlFGVEZXLUYkNiVGZnMtRi02M1ExJkludmlzaWJsZVRpbWVzO0YoRmd0RmFzRjZGYnNGZHNGZXNGQEZCRkVGSEZKRkxGTkZRRlRGVy1GJDYjLUkmbXNxcnRHRiU2JS1GaW82OVEjNDFGKEZccEZfcEZhcEZjcEZlcEZncEZpcEZbcUZdcUZfcUZhcUZjcUZlcUZncUZpcUZbckZdckZfckZhckZjckZlckZockZURldGanNGXXRGYHRGYnRGVEZXLUZebzYjLUZhbzYjRmhzLUYtNjNRIl1GKC9GMVEocG9zdGZpeEYoRjNGNkY5L0Y9UTJ2ZXJ5dGhpbm1hdGhzcGFjZUYoRj5GQEZCRkVGSEZKRkxGTkZRRlRGVy1GLTYzUSIsRihGZ3RGYXMvRjdGNUZicy9GPVEzdmVyeXRoaWNrbWF0aHNwYWNlRihGZXNGQEZCRkVGSEZKRkxGTkZRRlRGVy1GJDYlRiwtRiQ2Iy1GW282JC1GXm82Iy1GYW82Iy1GZG82KkZmby1GJDYjLUYkNiZGXnNGZnMtRi02M1EoJm1pbnVzO0YoRmd0RmFzRjZGaXRGW3VGZXNGQEZCRkVGSEZKRkxGTkZRRlRGV0ZcdUZqc0ZddEZgdEZidEZURldGaXVGXXZGXXY3Iy1fRilJLG1wcmludHNsYXNoR0YoNiQ3IzckLUknUlRBQkxFR0YoNiUiKDMlcEMtSSdNQVRSSVhHRig2IzckNyMsJComIiIjIiIiLCYjRmB5Rl95ISIiKiZGYnlGYHkpIiNURmJ5RmB5RmB5RmN5RmB5NyNGYHkmSSdWZWN0b3JHNiQlKnByb3RlY3RlZEdGKjYjSSdjb2x1bW5HRigtRmV4NiUiKCdIMiMpLUZpeDYjNyQ3IywkKiZGX3lGYHksJkZieUZjeUZkeUZjeUZjeUZgeUZneUZoeTcjNyQtRmh5NiMvSSQlaWRHRihGZ3gtRmh5NiMvRl1bbEZgeg==
Visualizing eigenvectors and eigenvalues in R2
We can use eigenvectors and eigenvalues to give a description of the action of a transformation from R2 to R2 . We break this into three cases, when the eigenvalues are real and distinct, when the eigenvalues are complex (must be conjugate in this case), and when the eigenvalues repeat (must be real in this case).
Case 1 - Distinct real eigenvalues
To produce distinct real eigenvalues, we start with a matrix that is symmetric
A := RandomMatrix(2,2,generator=rand(-5..5),
outputoptions=[shape=symmetric]);
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5USJBRigvJSdmYW1pbHlHUTBUaW1lc35OZXd+Um9tYW5GKC8lJXNpemVHUSMxMkYoLyUlYm9sZEdRJmZhbHNlRigvJSdpdGFsaWNHUSV0cnVlRigvJSp1bmRlcmxpbmVHRjgvJSpzdWJzY3JpcHRHRjgvJSxzdXBlcnNjcmlwdEdGOC8lK2ZvcmVncm91bmRHUSpbMCwwLDI1NV1GKC8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdRigvJSdvcGFxdWVHRjgvJStleGVjdXRhYmxlR0Y4LyUpcmVhZG9ubHlHRjsvJSljb21wb3NlZEdGOC8lKmNvbnZlcnRlZEdGOC8lK2ltc2VsZWN0ZWRHRjgvJSxwbGFjZWhvbGRlckdGOC8lMGZvbnRfc3R5bGVfbmFtZUdRKjJEfk91dHB1dEYoLyUqbWF0aGNvbG9yR0ZELyUvbWF0aGJhY2tncm91bmRHRkcvJStmb250ZmFtaWx5R0YyLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGKC8lKW1hdGhzaXplR0Y1LUkjbW9HRiU2M1EjOj1GKC8lJWZvcm1HUSZpbmZpeEYoLyUmZmVuY2VHRjgvJSpzZXBhcmF0b3JHRjgvJSdsc3BhY2VHUS90aGlja21hdGhzcGFjZUYoLyUncnNwYWNlR0ZbcC8lKXN0cmV0Y2h5R0Y4LyUqc3ltbWV0cmljR0Y4LyUobWF4c2l6ZUdRKWluZmluaXR5RigvJShtaW5zaXplR1EiMUYoLyUobGFyZ2VvcEdGOC8lLm1vdmFibGVsaW1pdHNHRjgvJSdhY2NlbnRHRjgvJTBmb250X3N0eWxlX25hbWVHRlgvJSVzaXplR0Y1LyUrZm9yZWdyb3VuZEdGRC8lK2JhY2tncm91bmRHRkctRiQ2JS1GX282M1EiW0YoL0Zjb1EncHJlZml4RigvRmZvRjtGZ28vRmpvUS50aGlubWF0aHNwYWNlRigvRl1wRl9yL0ZfcEY7RmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiMtSSdtdGFibGVHRiU2JC1JJG10ckdGJTYkLUkkbXRkR0YlNiMtSSNtbkdGJTY5USI1RihGMEYzRjYvRjpGOEY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ24vRmpuUSdub3JtYWxGKEZcby1GW3M2Iy1GXnM2OVEiMkYoRjBGM0Y2RmFzRjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZic0Zcby1GaHI2JEZkcy1GW3M2Iy1GXnM2OVEiNEYoRjBGM0Y2RmFzRjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZic0Zcby1GX282M1EiXUYoL0Zjb1EocG9zdGZpeEYoRl1yRmdvRl5yL0ZdcFEydmVyeXRoaW5tYXRoc3BhY2VGKEZhckZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcTcjLV9GKUksbXByaW50c2xhc2hHRig2JDcjPkkiQUdGKC1JJ1JUQUJMRUdGKDYlIilzJzQhPS1JJ01BVFJJWEdGKDYjNyQ3JCIiJiIiIzckRml1IiIlSSdNYXRyaXhHNiQlKnByb3RlY3RlZEdGKjcjLUZcdjYjL0kkJWlkR0YoRmJ1
We use the Eigenvector command to find the eigenvalues and corresponding eigenvectors.
EigenA := [Eigenvectors(A)]:
EigenVecA := [Column( EigenA[2],1),Column( EigenA[2],2)];
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5USpFaWdlblZlY0FGKC8lJ2ZhbWlseUdRMFRpbWVzfk5ld35Sb21hbkYoLyUlc2l6ZUdRIzEyRigvJSVib2xkR1EmZmFsc2VGKC8lJ2l0YWxpY0dRJXRydWVGKC8lKnVuZGVybGluZUdGOC8lKnN1YnNjcmlwdEdGOC8lLHN1cGVyc2NyaXB0R0Y4LyUrZm9yZWdyb3VuZEdRKlswLDAsMjU1XUYoLyUrYmFja2dyb3VuZEdRLlsyNTUsMjU1LDI1NV1GKC8lJ29wYXF1ZUdGOC8lK2V4ZWN1dGFibGVHRjgvJSlyZWFkb25seUdGOy8lKWNvbXBvc2VkR0Y4LyUqY29udmVydGVkR0Y4LyUraW1zZWxlY3RlZEdGOC8lLHBsYWNlaG9sZGVyR0Y4LyUwZm9udF9zdHlsZV9uYW1lR1EqMkR+T3V0cHV0RigvJSptYXRoY29sb3JHRkQvJS9tYXRoYmFja2dyb3VuZEdGRy8lK2ZvbnRmYW1pbHlHRjIvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YoLyUpbWF0aHNpemVHRjUtSSNtb0dGJTYzUSM6PUYoLyUlZm9ybUdRJmluZml4RigvJSZmZW5jZUdGOC8lKnNlcGFyYXRvckdGOC8lJ2xzcGFjZUdRL3RoaWNrbWF0aHNwYWNlRigvJSdyc3BhY2VHRltwLyUpc3RyZXRjaHlHRjgvJSpzeW1tZXRyaWNHRjgvJShtYXhzaXplR1EpaW5maW5pdHlGKC8lKG1pbnNpemVHUSIxRigvJShsYXJnZW9wR0Y4LyUubW92YWJsZWxpbWl0c0dGOC8lJ2FjY2VudEdGOC8lMGZvbnRfc3R5bGVfbmFtZUdGWC8lJXNpemVHRjUvJStmb3JlZ3JvdW5kR0ZELyUrYmFja2dyb3VuZEdGRy1GJDYlLUZfbzYzUSJbRigvRmNvUSdwcmVmaXhGKC9GZm9GO0Znby9Gam9RLnRoaW5tYXRoc3BhY2VGKC9GXXBGX3IvRl9wRjtGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtRiQ2JS1GJDYlRmhxLUYkNiMtSSdtdGFibGVHRiU2JC1JJG10ckdGJTYjLUkkbXRkR0YlNiMtSSZtZnJhY0dGJTYqLUYkNiMtSSNtbkdGJTY5USIyRihGMEYzRjYvRjpGOEY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ24vRmpuUSdub3JtYWxGKEZcby1GJDYjLUYkNiYtRl9vNjNRKiZ1bWludXMwO0YoRltyRmVvRmdvL0Zqb1EkMGVtRigvRl1wRmV0Rl5wRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUZiczYqLUZnczY5RmdwRjBGM0Y2RmpzRjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZbdEZcb0Zmcy8lLmxpbmV0aGlja25lc3NHUSIxRigvJStkZW5vbWFsaWduR1EnY2VudGVyRigvJSludW1hbGlnbkdGYHUvJSliZXZlbGxlZEdGOEZicUZkcS1GX282M1EiK0YoRmJvRmVvRmdvL0Zqb1EwbWVkaXVtbWF0aHNwYWNlRigvRl1wRml1Rl5wRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiVGZ3QtRl9vNjNRMSZJbnZpc2libGVUaW1lcztGKEZib0Zlb0Znb0ZkdEZmdEZecEZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1GJDYjLUkmbXNxcnRHRiU2JS1GZ3M2OVEjMTdGKEYwRjNGNkZqc0Y8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GW3RGXG9GYnFGZHFGW3VGXnVGYXVGY3VGYnFGZHEtRlxzNiMtRl9zNiNGaXQtRl9vNjNRIl1GKC9GY29RKHBvc3RmaXhGKEZdckZnb0Zeci9GXXBRMnZlcnl0aGlubWF0aHNwYWNlRihGYXJGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtRl9vNjNRIixGKEZib0Zlby9GaG9GO0ZkdC9GXXBRM3Zlcnl0aGlja21hdGhzcGFjZUYoRl5wRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiVGaHEtRiQ2Iy1GaXI2JC1GXHM2Iy1GX3M2Iy1GYnM2KkZkcy1GJDYjLUYkNiZGYXRGZ3QtRl9vNjNRKCZtaW51cztGKEZib0Zlb0Znb0ZodUZqdUZecEZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcUZbdkZbdUZedUZhdUZjdUZicUZkcUZodkZcd0ZcdzcjLV9GKUksbXByaW50c2xhc2hHRig2JDcjPkkqRWlnZW5WZWNBR0YoNyQtSSdSVEFCTEVHRig2JSIpK1daPS1JJ01BVFJJWEdGKDYjNyQ3IywkKiYiIiMiIiIsJiNGYXpGYHohIiIqJkZjekZheikiIzxGY3pGYXpGYXpGZHpGYXo3I0ZheiZJJ1ZlY3Rvckc2JCUqcHJvdGVjdGVkR0YqNiNJJ2NvbHVtbkdGKC1GZnk2JSIpWyQpWz0tRmp5NiM3JDcjLCQqJkZgekZheiwmRmN6RmR6RmV6RmR6RmR6RmF6Rmh6Rml6NyM3JC1GaXo2Iy9JJCVpZEdGKEZoeS1GaXo2Iy9GXlxsRmFbbA==
We plot the eigenvectors along with their images. The original vectors are graphed with double arrow heads and the images with single sided arrow heads.
v1 := arrow(EigenVecA[1],color=blue, shape=arrow,thickness=5):
v2 := arrow(EigenVecA[2],color=red, shape=arrow,thickness=5):
Av1 := arrow(A.EigenVecA[1],color=navy, shape=harpoon,thickness=3):
Av2 := arrow(A.EigenVecA[2],color=green, shape=harpoon,thickness=3):
display({v1, v2, Av1, Av2},scaling=constrained);
NictSSdDVVJWRVNHNiQlKnByb3RlY3RlZEdJKF9zeXNsaWJHNiI2JzckNyQkIiIhRi1GLDckJCEra1N3MnkhIzUkIiIiRi03JTckJCErXjdAWW5GMSQiKyl6Nic0d0YxRi43JCQhK143QFlkRjEkIisvIylRIVIpRjEtSSZTVFlMRUdGJTYjSSxQQVRDSE5PR1JJREdGJS1JJ0NPTE9VUkdGJTYmSSRSR0JHRiUkIiorKysrIiEiKUYsRiwtSSpUSElDS05FU1NHRiU2IyIiJi1GJDYmNyVGKzckJCErSz8pUSE+ISIqJCIrKD1aJVFDRlQ3JCQhKyUpei1YO0ZUJCIrW09jYj1GVEY/LUZENiZGRkYsRkdGLC1GSzYjIiIkLUYkNiY3JUYrNyQkIitKPylRUylGVCQiKzdHYmhsRlQ3JCQiKyUpei0mUidGVCQiK11qVnBjRlRGPy1GRDYmRkYkIikhXERQIkZJRmpvJCIpdmlvYkZJRmhuLUYkNic3JEYrNyQkIisxa3ghRyJGVEYyNyU3JCQiK1o3QFkoKkYxJCIrLSMpUVMnKUYxRmFwNyQkIitENml1NUZUJCIrJ3o2J2Z0RjFGPy1GRDYmRkZGLEYsRkdGSi1JKFNDQUxJTkdHRig2I0ksQ09OU1RSQUlORURHRiU=
Notice that the images of the eigenvectors are scalar multiples of the eigenvectors.
Case 3 - Repeated eigenvectors
The third case is repeated real eigenvalues. Sometimes, in this case, Maple does not produce a basis of eigenvectors. We will come back to this case in later sections.
A := <<2,0> | <2,2>>;
Eigenvectors(A);
Eigenvectors(A,output='list');
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5USJBRigvJSdmYW1pbHlHUTBUaW1lc35OZXd+Um9tYW5GKC8lJXNpemVHUSMxMkYoLyUlYm9sZEdRJmZhbHNlRigvJSdpdGFsaWNHUSV0cnVlRigvJSp1bmRlcmxpbmVHRjgvJSpzdWJzY3JpcHRHRjgvJSxzdXBlcnNjcmlwdEdGOC8lK2ZvcmVncm91bmRHUSpbMCwwLDI1NV1GKC8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdRigvJSdvcGFxdWVHRjgvJStleGVjdXRhYmxlR0Y4LyUpcmVhZG9ubHlHRjsvJSljb21wb3NlZEdGOC8lKmNvbnZlcnRlZEdGOC8lK2ltc2VsZWN0ZWRHRjgvJSxwbGFjZWhvbGRlckdGOC8lMGZvbnRfc3R5bGVfbmFtZUdRKjJEfk91dHB1dEYoLyUqbWF0aGNvbG9yR0ZELyUvbWF0aGJhY2tncm91bmRHRkcvJStmb250ZmFtaWx5R0YyLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGKC8lKW1hdGhzaXplR0Y1LUkjbW9HRiU2M1EjOj1GKC8lJWZvcm1HUSZpbmZpeEYoLyUmZmVuY2VHRjgvJSpzZXBhcmF0b3JHRjgvJSdsc3BhY2VHUS90aGlja21hdGhzcGFjZUYoLyUncnNwYWNlR0ZbcC8lKXN0cmV0Y2h5R0Y4LyUqc3ltbWV0cmljR0Y4LyUobWF4c2l6ZUdRKWluZmluaXR5RigvJShtaW5zaXplR1EiMUYoLyUobGFyZ2VvcEdGOC8lLm1vdmFibGVsaW1pdHNHRjgvJSdhY2NlbnRHRjgvJTBmb250X3N0eWxlX25hbWVHRlgvJSVzaXplR0Y1LyUrZm9yZWdyb3VuZEdGRC8lK2JhY2tncm91bmRHRkctRiQ2JS1GX282M1EiW0YoL0Zjb1EncHJlZml4RigvRmZvRjtGZ28vRmpvUS50aGlubWF0aHNwYWNlRigvRl1wRl9yL0ZfcEY7RmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiMtSSdtdGFibGVHRiU2JC1JJG10ckdGJTYkLUkkbXRkR0YlNiMtSSNtbkdGJTY5USIyRihGMEYzRjYvRjpGOEY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ24vRmpuUSdub3JtYWxGKEZcb0Zqci1GaHI2JC1GW3M2Iy1GXnM2OVEiMEYoRjBGM0Y2RmFzRjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZic0Zcb0Zqci1GX282M1EiXUYoL0Zjb1EocG9zdGZpeEYoRl1yRmdvRl5yL0ZdcFEydmVyeXRoaW5tYXRoc3BhY2VGKEZhckZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcTcjLV9GKUksbXByaW50c2xhc2hHRig2JDcjPkkiQUdGKC1JJ1JUQUJMRUdGKDYlIikhZT8jPS1JJ01BVFJJWEdGKDYjNyQ3JCIiI0ZjdTckIiIhRmN1SSdNYXRyaXhHNiQlKnByb3RlY3RlZEdGKjcjLUZmdTYjL0kkJWlkR0YoRl11
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtRiQ2JS1JI21vR0YlNjNRIltGKC8lJWZvcm1HUSdwcmVmaXhGKC8lJmZlbmNlR1EldHJ1ZUYoLyUqc2VwYXJhdG9yR1EmZmFsc2VGKC8lJ2xzcGFjZUdRLnRoaW5tYXRoc3BhY2VGKC8lJ3JzcGFjZUdGPS8lKXN0cmV0Y2h5R0Y3LyUqc3ltbWV0cmljR0Y6LyUobWF4c2l6ZUdRKWluZmluaXR5RigvJShtaW5zaXplR1EiMUYoLyUobGFyZ2VvcEdGOi8lLm1vdmFibGVsaW1pdHNHRjovJSdhY2NlbnRHRjovJTBmb250X3N0eWxlX25hbWVHUSoyRH5PdXRwdXRGKC8lJXNpemVHUSMxMkYoLyUrZm9yZWdyb3VuZEdRKlswLDAsMjU1XUYoLyUrYmFja2dyb3VuZEdRLlsyNTUsMjU1LDI1NV1GKC1GJDYjLUknbXRhYmxlR0YlNiQtSSRtdHJHRiU2Iy1JJG10ZEdGJTYjLUkjbW5HRiU2OVEiMkYoLyUnZmFtaWx5R1EwVGltZXN+TmV3flJvbWFuRigvJSVzaXplR0ZVLyUlYm9sZEdGOi8lJ2l0YWxpY0dGOi8lKnVuZGVybGluZUdGOi8lKnN1YnNjcmlwdEdGOi8lLHN1cGVyc2NyaXB0R0Y6LyUrZm9yZWdyb3VuZEdGWC8lK2JhY2tncm91bmRHRmVuLyUnb3BhcXVlR0Y6LyUrZXhlY3V0YWJsZUdGOi8lKXJlYWRvbmx5R0Y3LyUpY29tcG9zZWRHRjovJSpjb252ZXJ0ZWRHRjovJStpbXNlbGVjdGVkR0Y6LyUscGxhY2Vob2xkZXJHRjovJTBmb250X3N0eWxlX25hbWVHRlIvJSptYXRoY29sb3JHRlgvJS9tYXRoYmFja2dyb3VuZEdGZW4vJStmb250ZmFtaWx5R0Znby8lLG1hdGh2YXJpYW50R1Enbm9ybWFsRigvJSltYXRoc2l6ZUdGVUZbby1GLzYzUSJdRigvRjNRKHBvc3RmaXhGKEY1RjhGOy9GP1EydmVyeXRoaW5tYXRoc3BhY2VGKEZARkJGREZHRkpGTEZORlBGU0ZWRlktRi82M1EiLEYoL0YzUSZpbmZpeEYoL0Y2RjovRjlGNy9GPFEkMGVtRigvRj9RM3Zlcnl0aGlja21hdGhzcGFjZUYoL0ZBRjpGQkZERkdGSkZMRk5GUEZTRlZGWS1GJDYlRi4tRiQ2Iy1GaW42JC1GXG82JC1GX282Iy1GYm82OUZJRmVvRmhvRmpvRlxwRl5wRmBwRmJwRmRwRmZwRmhwRmpwRlxxRl5xRmBxRmJxRmRxRmZxRmhxRmpxRlxyRl5yRmFyLUZfbzYjLUZibzY5USIwRihGZW9GaG9Gam9GXHBGXnBGYHBGYnBGZHBGZnBGaHBGanBGXHFGXnFGYHFGYnFGZHFGZnFGaHFGanFGXHJGXnJGYXItRlxvNiRGYnRGYnRGY3I3Iy1fRilJLG1wcmludHNsYXNoR0YoNiQ3JC1JJ1JUQUJMRUdGKDYlIilfYz88LUknTUFUUklYR0YoNiM3JDcjIiIjRmd1JkknVmVjdG9yRzYkJSpwcm90ZWN0ZWRHRio2I0knY29sdW1uR0YoLUZgdTYlIilrcyZ6Ii1GZHU2IzckNyQiIiIiIiE3JEZndkZndkknTWF0cml4R0ZbdjckLUZpdTYjL0kkJWlkR0YoRmJ1LUZpdjYjL0Zed0Zhdg==
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtb0dGJTYzUSJbRigvJSVmb3JtR1EncHJlZml4RigvJSZmZW5jZUdRJXRydWVGKC8lKnNlcGFyYXRvckdRJmZhbHNlRigvJSdsc3BhY2VHUS50aGlubWF0aHNwYWNlRigvJSdyc3BhY2VHRjsvJSlzdHJldGNoeUdGNS8lKnN5bW1ldHJpY0dGOC8lKG1heHNpemVHUSlpbmZpbml0eUYoLyUobWluc2l6ZUdRIjFGKC8lKGxhcmdlb3BHRjgvJS5tb3ZhYmxlbGltaXRzR0Y4LyUnYWNjZW50R0Y4LyUwZm9udF9zdHlsZV9uYW1lR1EqMkR+T3V0cHV0RigvJSVzaXplR1EjMTJGKC8lK2ZvcmVncm91bmRHUSpbMCwwLDI1NV1GKC8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdRigtRiQ2Iy1GJDYlRiwtRiQ2Jy1JI21uR0YlNjlRIjJGKC8lJ2ZhbWlseUdRMFRpbWVzfk5ld35Sb21hbkYoLyUlc2l6ZUdGUy8lJWJvbGRHRjgvJSdpdGFsaWNHRjgvJSp1bmRlcmxpbmVHRjgvJSpzdWJzY3JpcHRHRjgvJSxzdXBlcnNjcmlwdEdGOC8lK2ZvcmVncm91bmRHRlYvJStiYWNrZ3JvdW5kR0ZZLyUnb3BhcXVlR0Y4LyUrZXhlY3V0YWJsZUdGOC8lKXJlYWRvbmx5R0Y1LyUpY29tcG9zZWRHRjgvJSpjb252ZXJ0ZWRHRjgvJStpbXNlbGVjdGVkR0Y4LyUscGxhY2Vob2xkZXJHRjgvJTBmb250X3N0eWxlX25hbWVHRlAvJSptYXRoY29sb3JHRlYvJS9tYXRoYmFja2dyb3VuZEdGWS8lK2ZvbnRmYW1pbHlHRmBvLyUsbWF0aHZhcmlhbnRHUSdub3JtYWxGKC8lKW1hdGhzaXplR0ZTLUYtNjNRIixGKC9GMVEmaW5maXhGKC9GNEY4L0Y3RjUvRjpRJDBlbUYoL0Y9UTN2ZXJ5dGhpY2ttYXRoc3BhY2VGKC9GP0Y4RkBGQkZFRkhGSkZMRk5GUUZURldGam5GXHItRiQ2JS1GLTYzUSJ8ZnJGKEYwRjNGNkY5RjxGPkZARkJGRUZIRkpGTEZORlFGVEZXLUYkNiMtRiQ2JUYsLUYkNiMtSSdtdGFibGVHRiU2JC1JJG10ckdGJTYjLUkkbXRkR0YlNiMtRltvNjlGR0Zeb0Zhb0Zjb0Zlb0Znb0Zpb0ZbcEZdcEZfcEZhcEZjcEZlcEZncEZpcEZbcUZdcUZfcUZhcUZjcUZlcUZncUZqcS1GZ3M2Iy1GanM2Iy1GW282OVEiMEYoRl5vRmFvRmNvRmVvRmdvRmlvRltwRl1wRl9wRmFwRmNwRmVwRmdwRmlwRltxRl1xRl9xRmFxRmNxRmVxRmdxRmpxLUYtNjNRIl1GKC9GMVEocG9zdGZpeEYoRjNGNkY5L0Y9UTJ2ZXJ5dGhpbm1hdGhzcGFjZUYoRj5GQEZCRkVGSEZKRkxGTkZRRlRGVy1GLTYzUSJ8aHJGKEZodEYzRjZGOUZqdEY+RkBGQkZFRkhGSkZMRk5GUUZURldGZXRGZXQ3Iy1fRilJLG1wcmludHNsYXNoR0YoNiQ3IzcjNyUiIiNGZ3U8Iy1JJ1JUQUJMRUdGKDYlIilzbD88LUknTUFUUklYR0YoNiM3JDcjIiIiNyMiIiEmSSdWZWN0b3JHNiQlKnByb3RlY3RlZEdGKjYjSSdjb2x1bW5HRig3IzcjNyVGZ3VGZ3U8Iy1GZXY2Iy9JJCVpZEdGKEZcdg==
Student[LinearAlgebra][EigenPlot](A, showunitvectors=true,
scaling=constrained, axes=normal, numvectors=8);
NlAtSSdDVVJWRVNHNiQlKnByb3RlY3RlZEdJKF9zeXNsaWJHNiI2JTdVNyQkIiM1ISIiJCIiIUYwNyQkIi9bOThxOUAqKiEjOSQiL0lrTkJMYDdGNDckJCIvakc2OyRlbypGNCQiLyZbOygpKSpvWyNGNDckJCIvRCkpZVt3KEgqRjQkIi9vJW9fWDdvJEY0NyQkIi8oUS8hbzFqKClGNCQiL3IsVG5gPFtGNDckJCIvJlxQJSpwLDQpRjQkIi9aI0hfX3koZUY0NyQkIi9TQHVpbypHKEY0JCIvcUdmNVpYb0Y0NyQkIi9wWygqKVJVUCdGNCQiL3p2RkM4MHhGNDckJCIvLnpcekVlYEY0JCIvKy1iI3pLVylGNDckJCIvNWw6SHpkVUY0JCIvLG1DMEZbISpGNDckJCIvJ1xQJSpwLDQkRjQkIi86Jkg7bDBeKkY0NyQkIi90JmU5OFEoPUY0JCIvcEcyRChHIykqRjQ3JCQiLzkkSCY+MHppISM6JCIvRkclR24tKSoqRjQ3JCQhL0AkSCY+MHppRmZvRmdvNyQkIS91JmU5OFEoPUY0RmFvNyQkIS8oXFAlKnAsNCRGNEZcbzckJCEvLGw6SHpkVUY0JCIvMG1DMEZbISpGNDckJCEvJip5XHpFZWBGNCQiLzAtYiN6S1cpRjQ3JCQhL2xbKCopUlVQJ0Y0JCIvI2V4VUtecShGNDckJCEvUUB1aW8qRyhGNCQiL3NHZjVaWG9GNDckJCEvJFxQJSpwLDQpRjQkIi9dI0hfX3koZUY0NyQkIS8mUS8hbzFqKClGNCQiL3UsVG5gPFtGNDckJCEvQykpZVt3KEgqRjQkIi9wJW9fWDdvJEY0NyQkIS9qRzY7JGVvKkY0JCIvJ1s7KCkpKm9bI0Y0NyQkIS9bOThxOUAqKkY0JCIvSmtOQkxgN0Y0NyQkISM1Ri4kIS9xaGN0YGhuISNHNyRGW3MkIS9La05CTGA3RjQ3JEZmciQhLylbOygpKSpvWyNGNDckRmFyJCEvciVvX1g3byRGNDckJCEvJVEvIW8xaigpRjQkIS92LFRuYDxbRjQ3JCQhLyNcUCUqcCw0KUY0JCEvXiNIX195KGVGNDckJCEvUEB1aW8qRyhGNCQhL3RHZjVaWG9GNDckJCEva1soKilSVVAnRjQkIS8kZXhVS15xKEY0NyQkIS8lKnlcekVlYEY0JCEvMC1iI3pLVylGNDckJCEvK2w6SHpkVUY0JCEvMG1DMEZbISpGNDckJCEvJ1tQJSpwLDQkRjQkIS89Jkg7bDBeKkY0NyQkIS9qJmU5OFEoPUY0JCEvckcyRChHIykqRjQ3JCQhLzMjSCY+MHppRmZvJCEvR0clR24tKSoqRjQ3JCQiL0cjSCY+MHppRmZvRml2NyQkIi9sJmU5OFEoPUY0JCEvcUcyRChHIykqRjQ3JCQiLylbUCUqcCw0JEY0Rl92NyQkIi8tbDpIemRVRjRGanU3JCQiLyYqeVx6RWVgRjQkIS8vLWIjektXKUY0NyQkIi9tWygqKVJVUCdGNCQhLyNleFVLXnEoRjQ3JCQiL1JAdWlvKkcoRjQkIS9yR2Y1WlhvRjQ3JCQiLyRcUCUqcCw0KUY0JCEvXSNIX195KGVGNDckJCIvJlEvIW8xaigpRjQkIS90LFRuYDxbRjQ3JEY9JCEvcCVvX1g3byRGNDckRjgkIS8nWzsoKSkqb1sjRjQ3JEYyJCEvSWtOQkxgN0Y0NyRGLCQiL01LcnVJXzghI0YtSSdDT0xPVVJHRiU2JkkkUkdCR0YlRjBGMEYwLUkqVEhJQ0tORVNTR0YlNiMiIiItSSlQT0xZR09OU0dGJTYmNyY3JCRGZnpGMCQhLysrKysrKz9GZm83JEZcW2wkIi8rKysrKys/RmZvNyQkIi8rKysrKytFISM4RmBbbDckRmNbbEZdW2w3JTckRmNbbCQiLysrKysrK1NGZm83JCQiIiRGMEYvNyRGY1tsJCEvKysrKysrU0Zmby1JJlNUWUxFR0YlNiNJLFBBVENITk9HUklER0YlLUZgejYmRmJ6JCIqKysrKyIhIilGL0YvLUYkNiM3KkZbW2xGX1tsRmJbbEZoW2xGW1xsRl5cbEZmW2xGW1tsLUZoejYmNyY3JCRGLkYwRmBbbDckRmFdbEZdW2w3JCQhLysrKysrK0VGZVtsRl1bbDckRmRdbEZgW2w3JTckRmRdbEZfXGw3JCQhIiRGMCQhIiFGMDckRmRdbEZpW2xGYVxsRmVcbC1GJDYjNypGYF1sRmJdbEZjXWxGaF1sRmldbEZeXmxGZl1sRmBdbC1GaHo2JkZfXWw3JUZoXWw3JEZqXWwkIi80dloiUThtJyEjSEZeXmxGYVxsLUZgejYmRmJ6Ri9GL0ZnXGwtRiQ2IzcqRmBdbEZiXWxGY11sRmhdbEZlXmxGXl5sRmZdbEZgXWwtRmh6NiZGano3JUZoW2w3JEZcXGwkIS80dloiUThtJ0ZoXmxGXlxsRmFcbEZpXmwtRiQ2IzcqRltbbEZfW2xGYltsRmhbbEZhX2xGXlxsRmZbbEZbW2wtRiQ2JTckNyRGXFtsRi83JEZhXWxGXF5sRmVcbEZjei1GJDYlNyQ3JEZhXWwkIi8uRFxnVz9BRmhebDckRlxbbCQhLy5EXGdXP0FGaF5sRmlebEZjei1GaHo2JjcmNyQkIi93Q0kwXmdyRjQkIi9WUnVCPSMqb0Y0NyQkIi9JLiM0RDspcEY0JCIvcHldSyYqXHNGNDckJCIvLSs0VSE0J0hGZVtsJCIvSzR2aE9jPUZlW2w3JCQiLzwjR3YjenlIRmVbbCQiL1JYKDMqZT89RmVbbDclNyQkIi8mKjNQKmY+JkhGZVtsJCIvRyIqPVpEdT1GZVtsNyQkIi8tJ2UhUmBOTkZlW2wkIi9zeFYuS0BARmVbbDckJCIvQ3RDcXQoKUhGZVtsJCIvVmpWMHEtPUZlW2xGYVxsLUZgejYmRmJ6Ri9GZ1xsRi8tRiQ2IzcqRmhgbEZdYWxGYmFsRl1ibEZiYmxGZ2JsRmdhbEZoYGwtRmh6NiY3JjckJCIveitVTkA5OUZmbyQiL3N4Vid5JmUpKkY0NyQkIS8mR0llOFVUIkZmbyQiL0JpTkA5OTVGZVtsNyQkIi87JVEneSZlZSJGZVtsJCIvQmlOQDk5RUZlW2w3JCQiLz40TkA5OTtGZVtsJCIveFBreSZlZSNGZVtsNyU3JCQiL2tAR2RycjpGZVtsJCIvWUNyVUdHRUZlW2w3JCQiL1pRKioqKioqKio+RmVbbEZcXGw3JCQiL3JycVVHRztGZVtsJCIvYXZHZHJyREZlW2xGYVxsRlxjbC1GJDYjNypGZGNsRmljbEZeZGxGaWRsRl5lbEZhZWxGY2RsRmRjbC1GaHo2JjcmNyQkIS80Lzl5MXJvRjQkIi9KIik0eTFycUY0NyQkIS80Lzl5MXJzRjQkIi9wYzR5MXJxRjQ3JCQhL0IrQHkxcnNGNCQiL3JdSXdaUT1GZVtsNyQkIS9CK0B5MXJvRjQkIi88YEl3WlE9RmVbbDclNyQkIS9CK0B5MXJ1RjQkIi9bXEl3WlE9RmVbbDckJCEvRnVBeTFycUY0JCIvcSFITT84NyNGZVtsNyQkIS9CK0B5MXJtRjQkIi9TYUl3WlE9RmVbbEZhXGxGXGNsLUYkNiM3KkZcZmxGYWZsRmZmbEZhZ2xGZmdsRltobEZbZ2xGXGZsLUZoejYmNyY3JCQhLyMzKysrKysiRmVbbCQiLyR6KmUqKioqKio+RmZvNyQkIS8hPSoqKioqKioqKioqRjQkIS8yLVQrKys/RmZvNyQkIS8iWzErKytnI0ZlW2wkIS9RbDEsKys/RmZvNyQkIS9YbSsrKytFRmVbbCQiL2lNJCopKioqKio+RmZvNyU3JCQhLypSMSsrK2cjRmVbbCQhL1FsMSwrK1NGZm83JCQhLy8jMysrKyskRmVbbCQhL0BoJUc/MUIiISNBNyQkIS9GbisrKytFRmVbbCQiL2lNJCopKioqKipSRmZvRmFcbEZcY2wtRiQ2IzcqRmZobEZbaWxGYGlsRltqbEZgamxGZmpsRmVpbEZmaGwtRmh6NiY3JjckJCEva1BGMF5nckY0JCEvJzR0UCM9IypvRjQ3JCQhL0k1KjNEOylwRjQkIS9HbmBLJipcc0Y0NyQkIS9zcTNVITQnSEZlW2wkIS9lJWU8bWomPUZlW2w3JCQhL1lgX0Z6eUhGZVtsJCEvJjQjKTMqZT89RmVbbDclNyQkIS9Pek8qZj4mSEZlW2wkIS9TbT5aRHU9RmVbbDckJCEvLWQwUmBOTkZlW2wkIS90a1cuS0BARmVbbDckJCEvJFtXLVB4KUhGZVtsJCEvOVJXMHEtPUZlW2xGYVxsRlxjbC1GJDYjNypGYVttRmZbbUZbXG1GZlxtRltdbUZgXW1GYFxtRmFbbS1GaHo2JjcmNyQkIS9zRixOQDk5RmZvJCEvaSFRa3kmZSkqRjQ3JCQiLyNSVmk4VVQiRmZvJCEvJT5jOFVULCJGZVtsNyQkIS9AeGl5JmVlIkZlW2wkIS8lPmM4VVRoI0ZlW2w3JCQhLyRHUzhVVGgiRmVbbCQhLzFRa3kmZWUjRmVbbDclNyQkIS9UOUZkcnI6RmVbbCQhLylRN0YlR0dFRmVbbDckJCEvVDopKioqKioqKj5GZVtsRmpdbDckJCEva2xwVUdHO0ZlW2wkIS83d0dkcnJERmVbbEZhXGxGXGNsLUYkNiM3KkZbXm1GYF5tRmVebUZgX21GZV9tRmhfbUZqXm1GW15tLUZoejYmNyY3JCQiLzolcCJ5MXJvRjQkIS9sMjJ5MXJxRjQ3JCQiLzolcCJ5MXJzRjQkIS9CXTF5MXJxRjQ3JCQiL1w9THkxcnNGNCQhL210SHdaUT1GZVtsNyQkIi9cPUx5MXJvRjQkIS9Tekh3WlE9RmVbbDclNyQkIi9cPUx5MXJ1RjQkIS95cUh3WlE9RmVbbDckJCIvZENQeTFycUY0JCEvby5VLktAQEZlW2w3JCQiL1w9THkxcm1GNCQhL0YjKUh3WlE9RmVbbEZhXGxGXGNsLUYkNiM3KkZjYG1GaGBtRl1hbUZoYW1GXWJtRmJibUZiYW1GY2BtLUZoejYmNyY3JCQiL2ssKysrKzVGZVtsJCEvJ2V6IioqKioqKj5GZm83JCQiL2YkKSoqKioqKioqKipGNCQiLzkvIzMrKysjRmZvNyQkIi9pSCwrKytFRmVbbCQiL3dJOC0rKz9GZm83JCQiLyFIOCsrK2cjRmVbbCQhL0NwJ3kqKioqKj5GZm83JTckJCIvKXo3KysrZyNGZVtsJCIvd0k4LSsrU0ZmbzckJCIvM2ssKysrSUZlW2wkIi9UQXAwQ2hDRmVqbDckJCIvYU0sKysrRUZlW2wkIS9DcCd5KioqKipSRmZvRmFcbEZcY2wtRiQ2IzcqRl1jbUZiY21GZ2NtRmJkbUZnZG1GXGVtRlxkbUZdY20tRmh6NiY3JjckJCIvIz1EYzhVVCJGZm8kIS8iR0FjOFVUIkZmbzckJCEvIz1EYzhVVCJGZm8kIi8iR0FjOFVUIkZmbzckJCIvL20jKkdWOmJGNCQiL2BIMWNGKXomRjQ3JCQiL1M7MGNGKXomRjQkIi8oXFEqR1Y6YkY0NyU3JCQiLycza2A2U1AmRjQkIi8iPUQncHBSZkY0NyQkIi8uOTZ5MXJxRjQkIi8xZjd5MXJxRjQ3JCQiL2VUaHBwUmZGNCQiL3BpUDosdWBGNEZhXGwtRmB6NiZGYnpGZ1xsRmdcbEYvLUYkNiM3KkZnZW1GXGZtRmFmbUZcZ21GYWdtRmZnbUZmZm1GZ2VtLUZoejYmNyY3JEZgW2wkIi9PUDp3MS1UISNENyRGXVtsJCEvT1A6dzEtVEZmaG03JCQhLyMzaywrKysjRmZvJCIvIWYqKioqKioqKip6RjQ3JCQiLz1mJCkqKioqKio+RmZvJCIvNS8rKysrISlGNDclNyQkIS8jM2ssKysrJUZmbyQiLyE9KioqKioqKioqekY0NyQkIS9vbzJRLl4/ISNCRlxbbDckJCIvPWYkKSoqKioqKlJGZm8kIi8/MysrKyshKUY0RmFcbEZbaG0tRiQ2IzcqRmNobUZnaG1GamhtRmVpbUZqaW1GXmptRl9pbUZjaG0tRmh6NiY3JjckJCIvIVE+YzhVVCJGZm8kIi8jM0djOFVUIkZmbzckJCEvIVE+YzhVVCJGZm8kIS8jM0djOFVUIkZmbzckJCEvbFUyY0YpeiZGNCQiLzdaIipHVjpiRjQ3JCQhLypRXSpHVjpiRjQkIi9HLi9jRil6JkY0NyU3JCQhLy5panBwUmZGNCQiLy8+TjosdWBGNDckJCEvNC85eTFycUY0JCIvK3A0eTFycUY0NyQkIS9eJSlROix1YEY0JCIvT0pncHBSZkY0RmFcbEZbaG0tRiQ2IzcqRmlqbUZeW25GY1tuRl5cbkZjXG5GaFxuRmhbbkZpam0tRmh6NiY3JjckJCEvcXVJXzgvIylGZmhtRmBbbDckJCIvcXVJXzgvIylGZmhtRl1bbDckJCEvIT0qKioqKioqKip6RjQkIS9tIkcuKysrI0ZmbzckJCEvPzMrKysrISlGNCQiL009bioqKioqKj5GZm83JTckJCEvZiQpKioqKioqKip6RjQkIS9tIkcuKysrJUZmbzckRmFdbCRGaWhtRl1qbTckJCEvVDsrKysrISlGNCQiL009bioqKioqKlJGZm9GYVxsRltobS1GJDYjNypGY11uRmZdbkZpXW5GZF5uRmlebkZbX25GXl5uRmNdbi1GaHo2JjcmNyQkIS8jKTRqTkA5OUZmbyQiLyFbO2M4VVQiRmZvNyQkIi8jKTRqTkA5OUZmbyQhLyFbO2M4VVQiRmZvNyQkIS8+RyEqR1Y6YkY0JCEveGIzY0YpeiZGNDckJCEvOiFIZ3YjKXomRjQkIS8iR2kqR1Y6YkY0NyU3JCQhL0AoUmA2U1AmRjQkIS9Ec2twcFJmRjQ3JCQhLyhSIzN5MXJxRjQkIS83XDp5MXJxRjQ3JCQhLzhAZnBwUmZGNCQhL0wxUzosdWBGNEZhXGxGW2htLUYkNiM3KkZmX25GW2BuRmBgbkZbYW5GYGFuRmVhbkZlYG5GZl9uLUZoejYmNyY3JEZdW2wkIS8/aCVHPzFCIiEjQzckRmBbbCQiLz9oJUc/MUIiRmNibjckJCIvW0FcKysrP0ZmbyQhL3AoKSoqKioqKioqekY0NyQkIS9feF0qKioqKio+RmZvJCEvSjcrKysrISlGNDclNyQkIi9bQVwrKytTRmZvJCEvUnYqKioqKioqKnpGNDckJCIvLjFCOTVgaEZdam1GYV1sNyQkIS9feF0qKioqKipSRmZvJCEvaEMrKysrISlGNEZhXGxGW2htLUYkNiM3KkZgYm5GZGJuRmdibkZiY25GZ2NuRmpjbkZcY25GYGJuLUZoejYmNyY3JCQhL3pOaE5AOTlGZm8kIS8kKVFqTkA5OUZmbzckJCIvek5oTkA5OUZmbyQiLyQpUWpOQDk5RmZvNyQkIi8hKm80Y0YpeiZGNCQhL0Y0KilHVjpiRjQ3JCQiL3VUKCpHVjpiRjQkIS8ueCxjRil6JkY0NyU3JCQiL1sjZSdwcFJmRjQkIS9Rdks6LHVgRjQ3JCQiLzolcCJ5MXJxRjQkIS8lKnkxeTFycUY0NyQkIi87R1Q6LHVgRjQkIS8jNCJlcHBSZkY0RmFcbEZbaG0tRiQ2IzcqRmVkbkZqZG5GX2VuRmplbkZfZm5GZGZuRmRlbkZlZG4tRmh6NiY3JjckJCIvJVxoL0YzayJGY2JuRl1bbDckJCEvJVxoL0YzayJGY2JuRmBbbDckJCIvZiQpKioqKioqKip6RjQkIi9KamwrKys/RmZvNyQkIi9UOysrKyshKUY0JCEvcE9NKioqKioqPkZmbzclNyQkIi89bioqKioqKioqekY0JCIvSmpsKysrU0ZmbzckRlxbbCQiL3J1SV84LyMpRl1qbTckJCIvI0cuKysrKylGNCQhL3BPTSoqKioqKlJGZm9GYVxsRltobS1GJDYjNypGX2duRmJnbkZlZ25GYGhuRmVobkZoaG5GamduRl9nbi1JKFNDQUxJTkdHRig2I0ksQ09OU1RSQUlORURHRiUtSSpBWEVTU1RZTEVHRiU2I0knTk9STUFMR0YlLUkmVElUTEVHRiU2I1FMVGhlfkltYWdlc35vZn5Vbml0flZlY3RvcnN+YW5kfkVpZ2VudmVjdG9yc0Yo
Eigenvectors(A,output='list');
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtb0dGJTYzUSJbRigvJSVmb3JtR1EncHJlZml4RigvJSZmZW5jZUdRJXRydWVGKC8lKnNlcGFyYXRvckdRJmZhbHNlRigvJSdsc3BhY2VHUS50aGlubWF0aHNwYWNlRigvJSdyc3BhY2VHRjsvJSlzdHJldGNoeUdGNS8lKnN5bW1ldHJpY0dGOC8lKG1heHNpemVHUSlpbmZpbml0eUYoLyUobWluc2l6ZUdRIjFGKC8lKGxhcmdlb3BHRjgvJS5tb3ZhYmxlbGltaXRzR0Y4LyUnYWNjZW50R0Y4LyUwZm9udF9zdHlsZV9uYW1lR1EqMkR+T3V0cHV0RigvJSVzaXplR1EjMTJGKC8lK2ZvcmVncm91bmRHUSpbMCwwLDI1NV1GKC8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdRigtRiQ2Iy1GJDYlRiwtRiQ2Jy1JI21uR0YlNjlRIjJGKC8lJ2ZhbWlseUdRMFRpbWVzfk5ld35Sb21hbkYoLyUlc2l6ZUdGUy8lJWJvbGRHRjgvJSdpdGFsaWNHRjgvJSp1bmRlcmxpbmVHRjgvJSpzdWJzY3JpcHRHRjgvJSxzdXBlcnNjcmlwdEdGOC8lK2ZvcmVncm91bmRHRlYvJStiYWNrZ3JvdW5kR0ZZLyUnb3BhcXVlR0Y4LyUrZXhlY3V0YWJsZUdGOC8lKXJlYWRvbmx5R0Y1LyUpY29tcG9zZWRHRjgvJSpjb252ZXJ0ZWRHRjgvJStpbXNlbGVjdGVkR0Y4LyUscGxhY2Vob2xkZXJHRjgvJTBmb250X3N0eWxlX25hbWVHRlAvJSptYXRoY29sb3JHRlYvJS9tYXRoYmFja2dyb3VuZEdGWS8lK2ZvbnRmYW1pbHlHRmBvLyUsbWF0aHZhcmlhbnRHUSdub3JtYWxGKC8lKW1hdGhzaXplR0ZTLUYtNjNRIixGKC9GMVEmaW5maXhGKC9GNEY4L0Y3RjUvRjpRJDBlbUYoL0Y9UTN2ZXJ5dGhpY2ttYXRoc3BhY2VGKC9GP0Y4RkBGQkZFRkhGSkZMRk5GUUZURldGam5GXHItRiQ2JS1GLTYzUSJ8ZnJGKEYwRjNGNkY5RjxGPkZARkJGRUZIRkpGTEZORlFGVEZXLUYkNiMtRiQ2JUYsLUYkNiMtSSdtdGFibGVHRiU2JC1JJG10ckdGJTYjLUkkbXRkR0YlNiMtRltvNjlGR0Zeb0Zhb0Zjb0Zlb0Znb0Zpb0ZbcEZdcEZfcEZhcEZjcEZlcEZncEZpcEZbcUZdcUZfcUZhcUZjcUZlcUZncUZqcS1GZ3M2Iy1GanM2Iy1GW282OVEiMEYoRl5vRmFvRmNvRmVvRmdvRmlvRltwRl1wRl9wRmFwRmNwRmVwRmdwRmlwRltxRl1xRl9xRmFxRmNxRmVxRmdxRmpxLUYtNjNRIl1GKC9GMVEocG9zdGZpeEYoRjNGNkY5L0Y9UTJ2ZXJ5dGhpbm1hdGhzcGFjZUYoRj5GQEZCRkVGSEZKRkxGTkZRRlRGVy1GLTYzUSJ8aHJGKEZodEYzRjZGOUZqdEY+RkBGQkZFRkhGSkZMRk5GUUZURldGZXRGZXQ3Iy1fRilJLG1wcmludHNsYXNoR0YoNiQ3IzcjNyUiIiNGZ3U8Iy1JJ1JUQUJMRUdGKDYlIilrKVJzIi1JJ01BVFJJWEdGKDYjNyQ3IyIiIjcjIiIhJkknVmVjdG9yRzYkJSpwcm90ZWN0ZWRHRio2I0knY29sdW1uR0YoNyM3IzclRmd1Rmd1PCMtRmV2NiMvSSQlaWRHRihGXHY=
Visualizing eigenvectors and eigenvalues in R3
Case 1 - Distinct real eigenvalues
We can repeat the process for 3 by 3 matrices, but we need some computational trickery. Since solving cubics involves using complex numbers in a way that does not easily simplify, we evaluate the decimal approximation of the matrix of eigenvectors. The imaginary part is small enough to be roundoff error, so we eliminate it before graphing.
A := RandomMatrix(3,3,generator=rand(-5..5),
outputoptions=[shape=symmetric]);
EigenA := [Eigenvectors(A)]:
EigenA1 := evalf(EigenA[2]);
EigenA2 := map(Re,evalf(EigenA[2]));
EigenVecA := ([Column(EigenA2,1),Column(EigenA2,2),Column(EigenA2,3)]):
v1 := arrow(EigenVecA[1],color=blue, shape=arrow,thickness=3):
v2 := arrow(EigenVecA[2],color=red, shape=arrow,thickness=3):
v3 := arrow(EigenVecA[3],color=brown, shape=arrow,thickness=3):
Av1 := arrow(A.EigenVecA[1],color=navy, shape=harpoon,thickness=3):
Av2 := arrow(A.EigenVecA[2],color=green, shape=harpoon,thickness=3):
Av3 := arrow(A.EigenVecA[3],color=yellow, shape=harpoon,thickness=3):
display({v1, v2,v3, Av1, Av2,Av3},scaling=constrained, axes=normal);
Computing eigenvalues and eigenvectors - Step by step
While it is useful to have a simple Maple command for Eigenvalues and Eigenvectors, we would also like to step through the computations, so that we can do this by hand if needed.
First define a matrix and the identity of the same size.
A := <<-1,0,5>|<0,4,0 >|<5,0,-1>>;
I3 := IdentityMatrix(3);
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5USJBRigvJSdmYW1pbHlHUTBUaW1lc35OZXd+Um9tYW5GKC8lJXNpemVHUSMxMkYoLyUlYm9sZEdRJmZhbHNlRigvJSdpdGFsaWNHUSV0cnVlRigvJSp1bmRlcmxpbmVHRjgvJSpzdWJzY3JpcHRHRjgvJSxzdXBlcnNjcmlwdEdGOC8lK2ZvcmVncm91bmRHUSpbMCwwLDI1NV1GKC8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdRigvJSdvcGFxdWVHRjgvJStleGVjdXRhYmxlR0Y4LyUpcmVhZG9ubHlHRjsvJSljb21wb3NlZEdGOC8lKmNvbnZlcnRlZEdGOC8lK2ltc2VsZWN0ZWRHRjgvJSxwbGFjZWhvbGRlckdGOC8lMGZvbnRfc3R5bGVfbmFtZUdRKjJEfk91dHB1dEYoLyUqbWF0aGNvbG9yR0ZELyUvbWF0aGJhY2tncm91bmRHRkcvJStmb250ZmFtaWx5R0YyLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGKC8lKW1hdGhzaXplR0Y1LUkjbW9HRiU2M1EjOj1GKC8lJWZvcm1HUSZpbmZpeEYoLyUmZmVuY2VHRjgvJSpzZXBhcmF0b3JHRjgvJSdsc3BhY2VHUS90aGlja21hdGhzcGFjZUYoLyUncnNwYWNlR0ZbcC8lKXN0cmV0Y2h5R0Y4LyUqc3ltbWV0cmljR0Y4LyUobWF4c2l6ZUdRKWluZmluaXR5RigvJShtaW5zaXplR1EiMUYoLyUobGFyZ2VvcEdGOC8lLm1vdmFibGVsaW1pdHNHRjgvJSdhY2NlbnRHRjgvJTBmb250X3N0eWxlX25hbWVHRlgvJSVzaXplR0Y1LyUrZm9yZWdyb3VuZEdGRC8lK2JhY2tncm91bmRHRkctRiQ2JS1GX282M1EiW0YoL0Zjb1EncHJlZml4RigvRmZvRjtGZ28vRmpvUS50aGlubWF0aHNwYWNlRigvRl1wRl9yL0ZfcEY7RmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiMtSSdtdGFibGVHRiU2JS1JJG10ckdGJTYlLUkkbXRkR0YlNiMtSSNtbkdGJTY5USkmbWludXM7MUYoRjBGM0Y2L0Y6RjhGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduL0ZqblEnbm9ybWFsRihGXG8tRltzNiMtRl5zNjlRIjBGKEYwRjNGNkZhc0Y8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GYnNGXG8tRltzNiMtRl5zNjlRIjVGKEYwRjNGNkZhc0Y8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GYnNGXG8tRmhyNiVGZHMtRltzNiMtRl5zNjlRIjRGKEYwRjNGNkZhc0Y8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GYnNGXG9GZHMtRmhyNiVGaXNGZHNGanItRl9vNjNRIl1GKC9GY29RKHBvc3RmaXhGKEZdckZnb0Zeci9GXXBRMnZlcnl0aGlubWF0aHNwYWNlRihGYXJGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHE3Iy1fRilJLG1wcmludHNsYXNoR0YoNiQ3Iz5JIkFHRigtSSdSVEFCTEVHRig2JSIpc2YtPC1JJ01BVFJJWEdGKDYjNyU3JSEiIiIiISIiJjclRmB2IiIlRmB2NyVGYXZGYHZGX3ZJJ01hdHJpeEc2JCUqcHJvdGVjdGVkR0YqNyMtRmV2NiMvSSQlaWRHRihGaXU=
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5USNJM0YoLyUnZmFtaWx5R1EwVGltZXN+TmV3flJvbWFuRigvJSVzaXplR1EjMTJGKC8lJWJvbGRHUSZmYWxzZUYoLyUnaXRhbGljR1EldHJ1ZUYoLyUqdW5kZXJsaW5lR0Y4LyUqc3Vic2NyaXB0R0Y4LyUsc3VwZXJzY3JpcHRHRjgvJStmb3JlZ3JvdW5kR1EqWzAsMCwyNTVdRigvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XUYoLyUnb3BhcXVlR0Y4LyUrZXhlY3V0YWJsZUdGOC8lKXJlYWRvbmx5R0Y7LyUpY29tcG9zZWRHRjgvJSpjb252ZXJ0ZWRHRjgvJStpbXNlbGVjdGVkR0Y4LyUscGxhY2Vob2xkZXJHRjgvJTBmb250X3N0eWxlX25hbWVHUSoyRH5PdXRwdXRGKC8lKm1hdGhjb2xvckdGRC8lL21hdGhiYWNrZ3JvdW5kR0ZHLyUrZm9udGZhbWlseUdGMi8lLG1hdGh2YXJpYW50R1EnaXRhbGljRigvJSltYXRoc2l6ZUdGNS1JI21vR0YlNjNRIzo9RigvJSVmb3JtR1EmaW5maXhGKC8lJmZlbmNlR0Y4LyUqc2VwYXJhdG9yR0Y4LyUnbHNwYWNlR1EvdGhpY2ttYXRoc3BhY2VGKC8lJ3JzcGFjZUdGW3AvJSlzdHJldGNoeUdGOC8lKnN5bW1ldHJpY0dGOC8lKG1heHNpemVHUSlpbmZpbml0eUYoLyUobWluc2l6ZUdRIjFGKC8lKGxhcmdlb3BHRjgvJS5tb3ZhYmxlbGltaXRzR0Y4LyUnYWNjZW50R0Y4LyUwZm9udF9zdHlsZV9uYW1lR0ZYLyUlc2l6ZUdGNS8lK2ZvcmVncm91bmRHRkQvJStiYWNrZ3JvdW5kR0ZHLUYkNiUtRl9vNjNRIltGKC9GY29RJ3ByZWZpeEYoL0Zmb0Y7RmdvL0Zqb1EudGhpbm1hdGhzcGFjZUYoL0ZdcEZfci9GX3BGO0ZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1GJDYjLUknbXRhYmxlR0YlNiUtSSRtdHJHRiU2JS1JJG10ZEdGJTYjLUkjbW5HRiU2OUZncEYwRjNGNi9GOkY4RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbi9Gam5RJ25vcm1hbEYoRlxvLUZbczYjLUZeczY5USIwRihGMEYzRjZGYHNGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmFzRlxvRmNzLUZocjYlRmNzRmpyRmNzLUZocjYlRmNzRmNzRmpyLUZfbzYzUSJdRigvRmNvUShwb3N0Zml4RihGXXJGZ29GXnIvRl1wUTJ2ZXJ5dGhpbm1hdGhzcGFjZUYoRmFyRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxNyMtX0YpSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSNJM0dGKC1JJ1JUQUJMRUdGKDYlIiklZSQpbyItSSdNQVRSSVhHRig2IzclNyUiIiIiIiFGZXU3JUZldUZkdUZldTclRmV1RmV1RmR1SSdNYXRyaXhHNiQlKnByb3RlY3RlZEdGKjcjLUZodTYjL0kkJWlkR0YoRl51
Compute the characteristic polynomial, which is det(A - x*I3)
chpol := Determinant(A - x*I3);
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5USZjaHBvbEYoLyUnZmFtaWx5R1EwVGltZXN+TmV3flJvbWFuRigvJSVzaXplR1EjMTJGKC8lJWJvbGRHUSZmYWxzZUYoLyUnaXRhbGljR1EldHJ1ZUYoLyUqdW5kZXJsaW5lR0Y4LyUqc3Vic2NyaXB0R0Y4LyUsc3VwZXJzY3JpcHRHRjgvJStmb3JlZ3JvdW5kR1EqWzAsMCwyNTVdRigvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XUYoLyUnb3BhcXVlR0Y4LyUrZXhlY3V0YWJsZUdGOC8lKXJlYWRvbmx5R0Y7LyUpY29tcG9zZWRHRjgvJSpjb252ZXJ0ZWRHRjgvJStpbXNlbGVjdGVkR0Y4LyUscGxhY2Vob2xkZXJHRjgvJTBmb250X3N0eWxlX25hbWVHUSoyRH5PdXRwdXRGKC8lKm1hdGhjb2xvckdGRC8lL21hdGhiYWNrZ3JvdW5kR0ZHLyUrZm9udGZhbWlseUdGMi8lLG1hdGh2YXJpYW50R1EnaXRhbGljRigvJSltYXRoc2l6ZUdGNS1JI21vR0YlNjNRIzo9RigvJSVmb3JtR1EmaW5maXhGKC8lJmZlbmNlR0Y4LyUqc2VwYXJhdG9yR0Y4LyUnbHNwYWNlR1EvdGhpY2ttYXRoc3BhY2VGKC8lJ3JzcGFjZUdGW3AvJSlzdHJldGNoeUdGOC8lKnN5bW1ldHJpY0dGOC8lKG1heHNpemVHUSlpbmZpbml0eUYoLyUobWluc2l6ZUdRIjFGKC8lKGxhcmdlb3BHRjgvJS5tb3ZhYmxlbGltaXRzR0Y4LyUnYWNjZW50R0Y4LyUwZm9udF9zdHlsZV9uYW1lR0ZYLyUlc2l6ZUdGNS8lK2ZvcmVncm91bmRHRkQvJStiYWNrZ3JvdW5kR0ZHLUYkNiotRl9vNjNRKiZ1bWludXMwO0YoL0Zjb1EncHJlZml4RihGZW9GZ28vRmpvUSQwZW1GKC9GXXBGXnJGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtSSNtbkdGJTY5USM5NkYoRjBGM0Y2L0Y6RjhGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduL0ZqblEnbm9ybWFsRihGXG8tRl9vNjNRIitGKEZib0Zlb0Znby9Gam9RMG1lZGl1bW1hdGhzcGFjZUYoL0ZdcEZbc0ZecEZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1GJDYlLUZhcjY5USMzMkYoRjBGM0Y2RmRyRjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZlckZcby1GX282M1ExJkludmlzaWJsZVRpbWVzO0YoRmJvRmVvRmdvRl1yRl9yRl5wRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYtNjlRInhGKEYwRjNGNkY5RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZpbkZcb0Znci1GJDYlLUZhcjY5USIyRihGMEYzRjZGZHJGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmVyRlxvRmJzLUklbXN1cEdGJTYlRmVzRmpzLyUxc3VwZXJzY3JpcHRzaGlmdEdRIjBGKC1GX282M1EoJm1pbnVzO0YoRmJvRmVvRmdvRmpyRlxzRl5wRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiMtRl50NiVGZXMtRmFyNjlRIjNGKEYwRjNGNkZkckY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GZXJGXG9GYHQ3Iy1fRilJLG1wcmludHNsYXNoR0YoNiQ3Iz5JJmNocG9sR0YoLCoiIycqISIiKiYiI0siIiJJInhHRihGanVGanUqJiIiI0ZqdSlGW3ZGXXZGanVGanUqJClGW3YiIiRGanVGZ3U3I0ZldQ==
Set equal to zero and solve.
AEvals := [solve(chpol = 0,x)];
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5USdBRXZhbHNGKC8lJ2ZhbWlseUdRMFRpbWVzfk5ld35Sb21hbkYoLyUlc2l6ZUdRIzEyRigvJSVib2xkR1EmZmFsc2VGKC8lJ2l0YWxpY0dRJXRydWVGKC8lKnVuZGVybGluZUdGOC8lKnN1YnNjcmlwdEdGOC8lLHN1cGVyc2NyaXB0R0Y4LyUrZm9yZWdyb3VuZEdRKlswLDAsMjU1XUYoLyUrYmFja2dyb3VuZEdRLlsyNTUsMjU1LDI1NV1GKC8lJ29wYXF1ZUdGOC8lK2V4ZWN1dGFibGVHRjgvJSlyZWFkb25seUdGOy8lKWNvbXBvc2VkR0Y4LyUqY29udmVydGVkR0Y4LyUraW1zZWxlY3RlZEdGOC8lLHBsYWNlaG9sZGVyR0Y4LyUwZm9udF9zdHlsZV9uYW1lR1EqMkR+T3V0cHV0RigvJSptYXRoY29sb3JHRkQvJS9tYXRoYmFja2dyb3VuZEdGRy8lK2ZvbnRmYW1pbHlHRjIvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YoLyUpbWF0aHNpemVHRjUtSSNtb0dGJTYzUSM6PUYoLyUlZm9ybUdRJmluZml4RigvJSZmZW5jZUdGOC8lKnNlcGFyYXRvckdGOC8lJ2xzcGFjZUdRL3RoaWNrbWF0aHNwYWNlRigvJSdyc3BhY2VHRltwLyUpc3RyZXRjaHlHRjgvJSpzeW1tZXRyaWNHRjgvJShtYXhzaXplR1EpaW5maW5pdHlGKC8lKG1pbnNpemVHUSIxRigvJShsYXJnZW9wR0Y4LyUubW92YWJsZWxpbWl0c0dGOC8lJ2FjY2VudEdGOC8lMGZvbnRfc3R5bGVfbmFtZUdGWC8lJXNpemVHRjUvJStmb3JlZ3JvdW5kR0ZELyUrYmFja2dyb3VuZEdGRy1GJDYlLUZfbzYzUSJbRigvRmNvUSdwcmVmaXhGKC9GZm9GO0Znby9Gam9RLnRoaW5tYXRoc3BhY2VGKC9GXXBGX3IvRl9wRjtGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtRiQ2Jy1JI21uR0YlNjlRKSZtaW51czs2RihGMEYzRjYvRjpGOEY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ24vRmpuUSdub3JtYWxGKEZcby1GX282M1EiLEYoRmJvRmVvL0Zob0Y7L0Zqb1EkMGVtRigvRl1wUTN2ZXJ5dGhpY2ttYXRoc3BhY2VGKEZecEZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1GZXI2OVEiNEYoRjBGM0Y2RmhyRjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZpckZcb0Zbc0Zjcy1GX282M1EiXUYoL0Zjb1EocG9zdGZpeEYoRl1yRmdvRl5yL0ZdcFEydmVyeXRoaW5tYXRoc3BhY2VGKEZhckZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcTcjLV9GKUksbXByaW50c2xhc2hHRig2JDcjPkknQUV2YWxzR0YoNyUhIiciIiVGZ3Q3I0ZldA==
Compute the basis of the first eigenspace by finding the basis of the corresponding nullspace.
AEbasis1 := NullSpace(A-AEvals[1]*I3);
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5USlBRWJhc2lzMUYoLyUnZmFtaWx5R1EwVGltZXN+TmV3flJvbWFuRigvJSVzaXplR1EjMTJGKC8lJWJvbGRHUSZmYWxzZUYoLyUnaXRhbGljR1EldHJ1ZUYoLyUqdW5kZXJsaW5lR0Y4LyUqc3Vic2NyaXB0R0Y4LyUsc3VwZXJzY3JpcHRHRjgvJStmb3JlZ3JvdW5kR1EqWzAsMCwyNTVdRigvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XUYoLyUnb3BhcXVlR0Y4LyUrZXhlY3V0YWJsZUdGOC8lKXJlYWRvbmx5R0Y7LyUpY29tcG9zZWRHRjgvJSpjb252ZXJ0ZWRHRjgvJStpbXNlbGVjdGVkR0Y4LyUscGxhY2Vob2xkZXJHRjgvJTBmb250X3N0eWxlX25hbWVHUSoyRH5PdXRwdXRGKC8lKm1hdGhjb2xvckdGRC8lL21hdGhiYWNrZ3JvdW5kR0ZHLyUrZm9udGZhbWlseUdGMi8lLG1hdGh2YXJpYW50R1EnaXRhbGljRigvJSltYXRoc2l6ZUdGNS1JI21vR0YlNjNRIzo9RigvJSVmb3JtR1EmaW5maXhGKC8lJmZlbmNlR0Y4LyUqc2VwYXJhdG9yR0Y4LyUnbHNwYWNlR1EvdGhpY2ttYXRoc3BhY2VGKC8lJ3JzcGFjZUdGW3AvJSlzdHJldGNoeUdGOC8lKnN5bW1ldHJpY0dGOC8lKG1heHNpemVHUSlpbmZpbml0eUYoLyUobWluc2l6ZUdRIjFGKC8lKGxhcmdlb3BHRjgvJS5tb3ZhYmxlbGltaXRzR0Y4LyUnYWNjZW50R0Y4LyUwZm9udF9zdHlsZV9uYW1lR0ZYLyUlc2l6ZUdGNS8lK2ZvcmVncm91bmRHRkQvJStiYWNrZ3JvdW5kR0ZHLUYkNiUtRl9vNjNRInxmckYoL0Zjb1EncHJlZml4RigvRmZvRjtGZ28vRmpvUS50aGlubWF0aHNwYWNlRigvRl1wRl9yL0ZfcEY7RmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiMtRiQ2JS1GX282M1EiW0YoRltyRl1yRmdvRl5yRmByRmFyRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiMtSSdtdGFibGVHRiU2JS1JJG10ckdGJTYjLUkkbXRkR0YlNiMtSSNtbkdGJTY5USkmbWludXM7MUYoRjBGM0Y2L0Y6RjhGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduL0ZqblEnbm9ybWFsRihGXG8tRl9zNiMtRmJzNiMtRmVzNjlRIjBGKEYwRjNGNkZoc0Y8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaXNGXG8tRl9zNiMtRmJzNiMtRmVzNjlGZ3BGMEYzRjZGaHNGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmlzRlxvLUZfbzYzUSJdRigvRmNvUShwb3N0Zml4RihGXXJGZ29GXnIvRl1wUTJ2ZXJ5dGhpbm1hdGhzcGFjZUYoRmFyRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUZfbzYzUSJ8aHJGKEZbdUZdckZnb0ZeckZddUZhckZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcTcjLV9GKUksbXByaW50c2xhc2hHRig2JDcjPkkpQUViYXNpczFHRig8Iy1JJ1JUQUJMRUdGKDYlIih3MCwmLUknTUFUUklYR0YoNiM3JTcjISIiNyMiIiE3IyIiIiZJJ1ZlY3Rvckc2JCUqcHJvdGVjdGVkR0YqNiNJJ2NvbHVtbkdGKDcjPCMtRml2NiMvSSQlaWRHRihGXnY=
Repeat for the second eigenvalue
AEbasis2 := NullSpace(A-AEvals[2]*I3);
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5USlBRWJhc2lzMkYoLyUnZmFtaWx5R1EwVGltZXN+TmV3flJvbWFuRigvJSVzaXplR1EjMTJGKC8lJWJvbGRHUSZmYWxzZUYoLyUnaXRhbGljR1EldHJ1ZUYoLyUqdW5kZXJsaW5lR0Y4LyUqc3Vic2NyaXB0R0Y4LyUsc3VwZXJzY3JpcHRHRjgvJStmb3JlZ3JvdW5kR1EqWzAsMCwyNTVdRigvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XUYoLyUnb3BhcXVlR0Y4LyUrZXhlY3V0YWJsZUdGOC8lKXJlYWRvbmx5R0Y7LyUpY29tcG9zZWRHRjgvJSpjb252ZXJ0ZWRHRjgvJStpbXNlbGVjdGVkR0Y4LyUscGxhY2Vob2xkZXJHRjgvJTBmb250X3N0eWxlX25hbWVHUSoyRH5PdXRwdXRGKC8lKm1hdGhjb2xvckdGRC8lL21hdGhiYWNrZ3JvdW5kR0ZHLyUrZm9udGZhbWlseUdGMi8lLG1hdGh2YXJpYW50R1EnaXRhbGljRigvJSltYXRoc2l6ZUdGNS1JI21vR0YlNjNRIzo9RigvJSVmb3JtR1EmaW5maXhGKC8lJmZlbmNlR0Y4LyUqc2VwYXJhdG9yR0Y4LyUnbHNwYWNlR1EvdGhpY2ttYXRoc3BhY2VGKC8lJ3JzcGFjZUdGW3AvJSlzdHJldGNoeUdGOC8lKnN5bW1ldHJpY0dGOC8lKG1heHNpemVHUSlpbmZpbml0eUYoLyUobWluc2l6ZUdRIjFGKC8lKGxhcmdlb3BHRjgvJS5tb3ZhYmxlbGltaXRzR0Y4LyUnYWNjZW50R0Y4LyUwZm9udF9zdHlsZV9uYW1lR0ZYLyUlc2l6ZUdGNS8lK2ZvcmVncm91bmRHRkQvJStiYWNrZ3JvdW5kR0ZHLUYkNiUtRl9vNjNRInxmckYoL0Zjb1EncHJlZml4RigvRmZvRjtGZ28vRmpvUS50aGlubWF0aHNwYWNlRigvRl1wRl9yL0ZfcEY7RmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiUtRiQ2JS1GX282M1EiW0YoRltyRl1yRmdvRl5yRmByRmFyRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiMtSSdtdGFibGVHRiU2JS1JJG10ckdGJTYjLUkkbXRkR0YlNiMtSSNtbkdGJTY5RmdwRjBGM0Y2L0Y6RjhGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduL0ZqblEnbm9ybWFsRihGXG8tRl9zNiMtRmJzNiMtRmVzNjlRIjBGKEYwRjNGNkZnc0Y8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaHNGXG9GXnMtRl9vNjNRIl1GKC9GY29RKHBvc3RmaXhGKEZdckZnb0Zeci9GXXBRMnZlcnl0aGlubWF0aHNwYWNlRihGYXJGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtRl9vNjNRIixGKEZib0Zlby9GaG9GOy9Gam9RJDBlbUYoL0ZdcFEzdmVyeXRoaWNrbWF0aHNwYWNlRihGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtRiQ2JUZmci1GJDYjLUZcczYlRmpzRl5zRmpzRmF0LUZfbzYzUSJ8aHJGKEZkdEZdckZnb0ZeckZmdEZhckZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcTcjLV9GKUksbXByaW50c2xhc2hHRig2JDcjPkkpQUViYXNpczJHRig8JC1JJ1JUQUJMRUdGKDYlIihPNCwmLUknTUFUUklYR0YoNiM3JTcjIiIiNyMiIiFGanYmSSdWZWN0b3JHNiQlKnByb3RlY3RlZEdGKjYjSSdjb2x1bW5HRigtRmN2NiUiKDs1LCYtRmd2NiM3JUZcd0ZqdkZcd0ZedzcjPCQtRl53NiMvSSQlaWRHRihGZXYtRl53NiMvRl94RmZ3
Check that these eigenspace basis vectors are in fact eigenvectors of A
[ AEbasis1[1], AEvals[1], A . AEbasis1[1]];
[ AEbasis2[1], AEvals[2], A . AEbasis2[1]];
[ AEbasis2[2], AEvals[3], A . AEbasis2[2]];
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtb0dGJTYzUSJbRigvJSVmb3JtR1EncHJlZml4RigvJSZmZW5jZUdRJXRydWVGKC8lKnNlcGFyYXRvckdRJmZhbHNlRigvJSdsc3BhY2VHUS50aGlubWF0aHNwYWNlRigvJSdyc3BhY2VHRjsvJSlzdHJldGNoeUdGNS8lKnN5bW1ldHJpY0dGOC8lKG1heHNpemVHUSlpbmZpbml0eUYoLyUobWluc2l6ZUdRIjFGKC8lKGxhcmdlb3BHRjgvJS5tb3ZhYmxlbGltaXRzR0Y4LyUnYWNjZW50R0Y4LyUwZm9udF9zdHlsZV9uYW1lR1EqMkR+T3V0cHV0RigvJSVzaXplR1EjMTJGKC8lK2ZvcmVncm91bmRHUSpbMCwwLDI1NV1GKC8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdRigtRiQ2Jy1GJDYlRiwtRiQ2Iy1JJ210YWJsZUdGJTYlLUkkbXRyR0YlNiMtSSRtdGRHRiU2Iy1JI21uR0YlNjlRKSZtaW51czsxRigvJSdmYW1pbHlHUTBUaW1lc35OZXd+Um9tYW5GKC8lJXNpemVHRlMvJSVib2xkR0Y4LyUnaXRhbGljR0Y4LyUqdW5kZXJsaW5lR0Y4LyUqc3Vic2NyaXB0R0Y4LyUsc3VwZXJzY3JpcHRHRjgvJStmb3JlZ3JvdW5kR0ZWLyUrYmFja2dyb3VuZEdGWS8lJ29wYXF1ZUdGOC8lK2V4ZWN1dGFibGVHRjgvJSlyZWFkb25seUdGNS8lKWNvbXBvc2VkR0Y4LyUqY29udmVydGVkR0Y4LyUraW1zZWxlY3RlZEdGOC8lLHBsYWNlaG9sZGVyR0Y4LyUwZm9udF9zdHlsZV9uYW1lR0ZQLyUqbWF0aGNvbG9yR0ZWLyUvbWF0aGJhY2tncm91bmRHRlkvJStmb250ZmFtaWx5R0Zpby8lLG1hdGh2YXJpYW50R1Enbm9ybWFsRigvJSltYXRoc2l6ZUdGUy1GXm82Iy1GYW82Iy1GZG82OVEiMEYoRmdvRmpvRlxwRl5wRmBwRmJwRmRwRmZwRmhwRmpwRlxxRl5xRmBxRmJxRmRxRmZxRmhxRmpxRlxyRl5yRmByRmNyLUZebzYjLUZhbzYjLUZkbzY5RkdGZ29Gam9GXHBGXnBGYHBGYnBGZHBGZnBGaHBGanBGXHFGXnFGYHFGYnFGZHFGZnFGaHFGanFGXHJGXnJGYHJGY3ItRi02M1EiXUYoL0YxUShwb3N0Zml4RihGM0Y2RjkvRj1RMnZlcnl0aGlubWF0aHNwYWNlRihGPkZARkJGRUZIRkpGTEZORlFGVEZXLUYtNjNRIixGKC9GMVEmaW5maXhGKC9GNEY4L0Y3RjUvRjpRJDBlbUYoL0Y9UTN2ZXJ5dGhpY2ttYXRoc3BhY2VGKC9GP0Y4RkBGQkZFRkhGSkZMRk5GUUZURlctRmRvNjlRKSZtaW51czs2RihGZ29Gam9GXHBGXnBGYHBGYnBGZHBGZnBGaHBGanBGXHFGXnFGYHFGYnFGZHFGZnFGaHFGanFGXHJGXnJGYHJGY3JGaXMtRiQ2JUYsLUYkNiMtRltvNiUtRl5vNiMtRmFvNiMtRmRvNjlRIjZGKEZnb0Zqb0ZccEZecEZgcEZicEZkcEZmcEZocEZqcEZccUZecUZgcUZicUZkcUZmcUZocUZqcUZcckZeckZgckZjckZlci1GXm82Iy1GYW82I0ZldEZic0ZiczcjLV9GKUksbXByaW50c2xhc2hHRig2JDcjNyUtSSdSVEFCTEVHRig2JSIodzAsJi1JJ01BVFJJWEdGKDYjNyU3IyEiIjcjIiIhNyMiIiImSSdWZWN0b3JHNiQlKnByb3RlY3RlZEdGKjYjSSdjb2x1bW5HRighIictRmF2NiUiKSNcQ3IiLUZldjYjNyU3IyIiJ0ZqdjcjRmR3Rl53NyM3JS1GXnc2Iy9JJCVpZEdGKEZjdkZkdy1GXnc2Iy9GY3hGZ3c=
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtb0dGJTYzUSJbRigvJSVmb3JtR1EncHJlZml4RigvJSZmZW5jZUdRJXRydWVGKC8lKnNlcGFyYXRvckdRJmZhbHNlRigvJSdsc3BhY2VHUS50aGlubWF0aHNwYWNlRigvJSdyc3BhY2VHRjsvJSlzdHJldGNoeUdGNS8lKnN5bW1ldHJpY0dGOC8lKG1heHNpemVHUSlpbmZpbml0eUYoLyUobWluc2l6ZUdRIjFGKC8lKGxhcmdlb3BHRjgvJS5tb3ZhYmxlbGltaXRzR0Y4LyUnYWNjZW50R0Y4LyUwZm9udF9zdHlsZV9uYW1lR1EqMkR+T3V0cHV0RigvJSVzaXplR1EjMTJGKC8lK2ZvcmVncm91bmRHUSpbMCwwLDI1NV1GKC8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdRigtRiQ2Jy1GJDYlRiwtRiQ2Iy1JJ210YWJsZUdGJTYlLUkkbXRyR0YlNiMtSSRtdGRHRiU2Iy1JI21uR0YlNjlGRy8lJ2ZhbWlseUdRMFRpbWVzfk5ld35Sb21hbkYoLyUlc2l6ZUdGUy8lJWJvbGRHRjgvJSdpdGFsaWNHRjgvJSp1bmRlcmxpbmVHRjgvJSpzdWJzY3JpcHRHRjgvJSxzdXBlcnNjcmlwdEdGOC8lK2ZvcmVncm91bmRHRlYvJStiYWNrZ3JvdW5kR0ZZLyUnb3BhcXVlR0Y4LyUrZXhlY3V0YWJsZUdGOC8lKXJlYWRvbmx5R0Y1LyUpY29tcG9zZWRHRjgvJSpjb252ZXJ0ZWRHRjgvJStpbXNlbGVjdGVkR0Y4LyUscGxhY2Vob2xkZXJHRjgvJTBmb250X3N0eWxlX25hbWVHRlAvJSptYXRoY29sb3JHRlYvJS9tYXRoYmFja2dyb3VuZEdGWS8lK2ZvbnRmYW1pbHlHRmhvLyUsbWF0aHZhcmlhbnRHUSdub3JtYWxGKC8lKW1hdGhzaXplR0ZTLUZebzYjLUZhbzYjLUZkbzY5USIwRihGZm9GaW9GW3BGXXBGX3BGYXBGY3BGZXBGZ3BGaXBGW3FGXXFGX3FGYXFGY3FGZXFGZ3FGaXFGW3JGXXJGX3JGYnJGXW8tRi02M1EiXUYoL0YxUShwb3N0Zml4RihGM0Y2RjkvRj1RMnZlcnl0aGlubWF0aHNwYWNlRihGPkZARkJGRUZIRkpGTEZORlFGVEZXLUYtNjNRIixGKC9GMVEmaW5maXhGKC9GNEY4L0Y3RjUvRjpRJDBlbUYoL0Y9UTN2ZXJ5dGhpY2ttYXRoc3BhY2VGKC9GP0Y4RkBGQkZFRkhGSkZMRk5GUUZURlctRmRvNjlRIjRGKEZmb0Zpb0ZbcEZdcEZfcEZhcEZjcEZlcEZncEZpcEZbcUZdcUZfcUZhcUZjcUZlcUZncUZpcUZbckZdckZfckZickZicy1GJDYlRiwtRiQ2Iy1GW282JS1GXm82Iy1GYW82I0ZedEZkckZndEZbc0ZbczcjLV9GKUksbXByaW50c2xhc2hHRig2JDcjNyUtSSdSVEFCTEVHRig2JSIoTzQsJi1JJ01BVFJJWEdGKDYjNyU3IyIiIjcjIiIhRmp1JkknVmVjdG9yRzYkJSpwcm90ZWN0ZWRHRio2I0knY29sdW1uR0YoIiIlLUZjdTYlIihPP0MpLUZndTYjNyU3I0ZkdkZcdkZbd0ZedjcjNyUtRl52NiMvSSQlaWRHRihGZXVGZHYtRl52NiMvRmF3Rmd2
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtb0dGJTYzUSJbRigvJSVmb3JtR1EncHJlZml4RigvJSZmZW5jZUdRJXRydWVGKC8lKnNlcGFyYXRvckdRJmZhbHNlRigvJSdsc3BhY2VHUS50aGlubWF0aHNwYWNlRigvJSdyc3BhY2VHRjsvJSlzdHJldGNoeUdGNS8lKnN5bW1ldHJpY0dGOC8lKG1heHNpemVHUSlpbmZpbml0eUYoLyUobWluc2l6ZUdRIjFGKC8lKGxhcmdlb3BHRjgvJS5tb3ZhYmxlbGltaXRzR0Y4LyUnYWNjZW50R0Y4LyUwZm9udF9zdHlsZV9uYW1lR1EqMkR+T3V0cHV0RigvJSVzaXplR1EjMTJGKC8lK2ZvcmVncm91bmRHUSpbMCwwLDI1NV1GKC8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdRigtRiQ2Jy1GJDYlRiwtRiQ2Iy1JJ210YWJsZUdGJTYlLUkkbXRyR0YlNiMtSSRtdGRHRiU2Iy1JI21uR0YlNjlRIjBGKC8lJ2ZhbWlseUdRMFRpbWVzfk5ld35Sb21hbkYoLyUlc2l6ZUdGUy8lJWJvbGRHRjgvJSdpdGFsaWNHRjgvJSp1bmRlcmxpbmVHRjgvJSpzdWJzY3JpcHRHRjgvJSxzdXBlcnNjcmlwdEdGOC8lK2ZvcmVncm91bmRHRlYvJStiYWNrZ3JvdW5kR0ZZLyUnb3BhcXVlR0Y4LyUrZXhlY3V0YWJsZUdGOC8lKXJlYWRvbmx5R0Y1LyUpY29tcG9zZWRHRjgvJSpjb252ZXJ0ZWRHRjgvJStpbXNlbGVjdGVkR0Y4LyUscGxhY2Vob2xkZXJHRjgvJTBmb250X3N0eWxlX25hbWVHRlAvJSptYXRoY29sb3JHRlYvJS9tYXRoYmFja2dyb3VuZEdGWS8lK2ZvbnRmYW1pbHlHRmlvLyUsbWF0aHZhcmlhbnRHUSdub3JtYWxGKC8lKW1hdGhzaXplR0ZTLUZebzYjLUZhbzYjLUZkbzY5RkdGZ29Gam9GXHBGXnBGYHBGYnBGZHBGZnBGaHBGanBGXHFGXnFGYHFGYnFGZHFGZnFGaHFGanFGXHJGXnJGYHJGY3JGXW8tRi02M1EiXUYoL0YxUShwb3N0Zml4RihGM0Y2RjkvRj1RMnZlcnl0aGlubWF0aHNwYWNlRihGPkZARkJGRUZIRkpGTEZORlFGVEZXLUYtNjNRIixGKC9GMVEmaW5maXhGKC9GNEY4L0Y3RjUvRjpRJDBlbUYoL0Y9UTN2ZXJ5dGhpY2ttYXRoc3BhY2VGKC9GP0Y4RkBGQkZFRkhGSkZMRk5GUUZURlctRmRvNjlRIjRGKEZnb0Zqb0ZccEZecEZgcEZicEZkcEZmcEZocEZqcEZccUZecUZgcUZicUZkcUZmcUZocUZqcUZcckZeckZgckZjckZicy1GJDYlRiwtRiQ2Iy1GW282JUZdby1GXm82Iy1GYW82I0ZedEZdb0Zbc0ZbczcjLV9GKUksbXByaW50c2xhc2hHRig2JDcjNyUtSSdSVEFCTEVHRig2JSIoOzUsJi1JJ01BVFJJWEdGKDYjNyU3IyIiITcjIiIiRmp1JkknVmVjdG9yRzYkJSpwcm90ZWN0ZWRHRio2I0knY29sdW1uR0YoIiIlLUZjdTYlIihzU1YnLUZndTYjNyVGanU3I0ZkdkZqdUZedjcjNyUtRl52NiMvSSQlaWRHRihGZXVGZHYtRl52NiMvRmF3Rmd2
Note that these vectors are the correct multiples of the eigenvectors.
Exercises
6) Follow the recipe above to find the eigenvalues and eigenvectors of the matrix B.
B := Matrix(3,3,[[0,1,1],[0,2,0],[-2,1,3]]);
7) Repeat for the matrix C.
C := Matrix(3,3,[[4,-3,1],[4,-1,0],[1,7,-4]]);
When you are done: what can you say about the matrix C?
Using eigenvalues and eigenvectors to ease computation
With A defined above, we want to compute the product of a power of A times a vector.
We'll easily compute LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2OVEhRicvJSdmYW1pbHlHUTBUaW1lc35OZXd+Um9tYW5GJy8lJXNpemVHUSMxMkYnLyUlYm9sZEdRJmZhbHNlRicvJSdpdGFsaWNHUSV0cnVlRicvJSp1bmRlcmxpbmVHRjcvJSpzdWJzY3JpcHRHRjcvJSxzdXBlcnNjcmlwdEdGNy8lK2ZvcmVncm91bmRHUShbMCwwLDBdRicvJStiYWNrZ3JvdW5kR0ZDLyUnb3BhcXVlR0Y3LyUrZXhlY3V0YWJsZUdGNy8lKXJlYWRvbmx5R0Y3LyUpY29tcG9zZWRHRjcvJSpjb252ZXJ0ZWRHRjcvJStpbXNlbGVjdGVkR0Y3LyUscGxhY2Vob2xkZXJHRjcvJTBmb250X3N0eWxlX25hbWVHUSsyRH5Db21tZW50RicvJSptYXRoY29sb3JHRkMvJS9tYXRoYmFja2dyb3VuZEdGQy8lK2ZvbnRmYW1pbHlHRjEvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLyUpbWF0aHNpemVHRjQtRiM2JUYrLUklbXN1cEdGJDYlLUYsNjlRIkFGJ0YvRjJGNUY4RjtGPUY/RkFGREZGRkhGSkZMRk5GUEZSRlRGV0ZZRmVuRmduRmpuLUkjbW5HRiQ2OVEjMjVGJ0YvRjJGNS9GOUY3RjtGPUY/RkFGREZGRkhGSkZMRk5GUEZSRlRGV0ZZRmVuL0ZoblEnbm9ybWFsRidGam4vJTFzdXBlcnNjcmlwdHNoaWZ0R1EiMEYnRitGKw== * X where X = LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2OVEhRicvJSdmYW1pbHlHUTBUaW1lc35OZXd+Um9tYW5GJy8lJXNpemVHUSMxMkYnLyUlYm9sZEdRJmZhbHNlRicvJSdpdGFsaWNHUSV0cnVlRicvJSp1bmRlcmxpbmVHRjcvJSpzdWJzY3JpcHRHRjcvJSxzdXBlcnNjcmlwdEdGNy8lK2ZvcmVncm91bmRHUShbMCwwLDBdRicvJStiYWNrZ3JvdW5kR0ZDLyUnb3BhcXVlR0Y3LyUrZXhlY3V0YWJsZUdGNy8lKXJlYWRvbmx5R0Y3LyUpY29tcG9zZWRHRjcvJSpjb252ZXJ0ZWRHRjcvJStpbXNlbGVjdGVkR0Y3LyUscGxhY2Vob2xkZXJHRjcvJTBmb250X3N0eWxlX25hbWVHUSsyRH5Db21tZW50RicvJSptYXRoY29sb3JHRkMvJS9tYXRoYmFja2dyb3VuZEdGQy8lK2ZvbnRmYW1pbHlHRjEvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLyUpbWF0aHNpemVHRjQtSSVtc3VwR0YkNiUtRiM2JS1JI21vR0YkNjNRIltGJy8lJWZvcm1HUSdwcmVmaXhGJy8lJmZlbmNlR0Y6LyUqc2VwYXJhdG9yR0Y3LyUnbHNwYWNlR1EudGhpbm1hdGhzcGFjZUYnLyUncnNwYWNlR0ZecC8lKXN0cmV0Y2h5R0Y6LyUqc3ltbWV0cmljR0Y3LyUobWF4c2l6ZUdRKWluZmluaXR5RicvJShtaW5zaXplR1EiMUYnLyUobGFyZ2VvcEdGNy8lLm1vdmFibGVsaW1pdHNHRjcvJSdhY2NlbnRHRjcvJTBmb250X3N0eWxlX25hbWVHRlYvJSVzaXplR0Y0LyUrZm9yZWdyb3VuZEdGQy8lK2JhY2tncm91bmRHRkMtRiM2KC1JI21uR0YkNjlRIjJGJ0YvRjJGNS9GOUY3RjtGPUY/RkFGREZGRkhGSkZMRk5GUEZSRlRGV0ZZRmVuL0ZoblEnbm9ybWFsRidGam4tRmJvNjNRIixGJy9GZm9RJmluZml4RicvRmlvRjcvRltwRjovRl1wUSQwZW1GJy9GYHBRM3Zlcnl0aGlja21hdGhzcGFjZUYnL0ZicEY3RmNwRmVwRmhwRltxRl1xRl9xRmFxRmNxRmVxRmdxLUZccjY5USIzRidGL0YyRjVGX3JGO0Y9Rj9GQUZERkZGSEZKRkxGTkZQRlJGVEZXRllGZW5GYHJGam5GYnItRiM2JC1GYm82M1EqJnVtaW51czA7RidGZW9GZ3JGam9GaXIvRmBwRmpyRl1zRmNwRmVwRmhwRltxRl1xRl9xRmFxRmNxRmVxRmdxLUZccjY5RmpwRi9GMkY1Rl9yRjtGPUY/RkFGREZGRkhGSkZMRk5GUEZSRlRGV0ZZRmVuRmByRmpuRistRmJvNjNRIl1GJy9GZm9RKHBvc3RmaXhGJ0Zob0Zqb0ZccC9GYHBRMnZlcnl0aGlubWF0aHNwYWNlRidGYXBGY3BGZXBGaHBGW3FGXXFGX3FGYXFGY3FGZXFGZ3EtRiw2OVEidEYnRi9GMkY1RjhGO0Y9Rj9GQUZERkZGSEZKRkxGTkZQRlJGVEZXRllGZW5GZ25Gam4vJTFzdXBlcnNjcmlwdHNoaWZ0R1EiMEYnRis=
First note that the three eigenvectors which we computed earlier form a basis for LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2OVEhRicvJSdmYW1pbHlHUTBUaW1lc35OZXd+Um9tYW5GJy8lJXNpemVHUSMxMkYnLyUlYm9sZEdRJmZhbHNlRicvJSdpdGFsaWNHUSV0cnVlRicvJSp1bmRlcmxpbmVHRjcvJSpzdWJzY3JpcHRHRjcvJSxzdXBlcnNjcmlwdEdGNy8lK2ZvcmVncm91bmRHUShbMCwwLDBdRicvJStiYWNrZ3JvdW5kR0ZDLyUnb3BhcXVlR0Y3LyUrZXhlY3V0YWJsZUdGNy8lKXJlYWRvbmx5R0Y3LyUpY29tcG9zZWRHRjcvJSpjb252ZXJ0ZWRHRjcvJStpbXNlbGVjdGVkR0Y3LyUscGxhY2Vob2xkZXJHRjcvJTBmb250X3N0eWxlX25hbWVHUSsyRH5Db21tZW50RicvJSptYXRoY29sb3JHRkMvJS9tYXRoYmFja2dyb3VuZEdGQy8lK2ZvbnRmYW1pbHlHRjEvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLyUpbWF0aHNpemVHRjQtRiM2JUYrLUklbXN1cEdGJDYlLUYsNjlRIlJGJ0YvRjJGNUY4RjtGPUY/RkFGREZGRkhGSkZMRk5GUEZSRlRGV0ZZRmVuRmduRmpuLUkjbW5HRiQ2OVEiM0YnRi9GMkY1L0Y5RjdGO0Y9Rj9GQUZERkZGSEZKRkxGTkZQRlJGVEZXRllGZW4vRmhuUSdub3JtYWxGJ0Zqbi8lMXN1cGVyc2NyaXB0c2hpZnRHUSIwRidGK0Yr
AEMat := <AEbasis1[1] | AEbasis2[1] | AEbasis2[2]>;
Rank(AEMat);
Now find the coefficients of X with respect to the basis above.
X := <2,3,-1>;
Xcvec := LinearSolve(AEMat,X);
Now LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2OVEhRicvJSdmYW1pbHlHUTBUaW1lc35OZXd+Um9tYW5GJy8lJXNpemVHUSMxMkYnLyUlYm9sZEdRJmZhbHNlRicvJSdpdGFsaWNHUSV0cnVlRicvJSp1bmRlcmxpbmVHRjcvJSpzdWJzY3JpcHRHRjcvJSxzdXBlcnNjcmlwdEdGNy8lK2ZvcmVncm91bmRHUShbMCwwLDBdRicvJStiYWNrZ3JvdW5kR0ZDLyUnb3BhcXVlR0Y3LyUrZXhlY3V0YWJsZUdGNy8lKXJlYWRvbmx5R0Y3LyUpY29tcG9zZWRHRjcvJSpjb252ZXJ0ZWRHRjcvJStpbXNlbGVjdGVkR0Y3LyUscGxhY2Vob2xkZXJHRjcvJTBmb250X3N0eWxlX25hbWVHUSsyRH5Db21tZW50RicvJSptYXRoY29sb3JHRkMvJS9tYXRoYmFja2dyb3VuZEdGQy8lK2ZvbnRmYW1pbHlHRjEvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLyUpbWF0aHNpemVHRjQtRiM2JUYrLUklbXN1cEdGJDYlLUYsNjlRIkFGJ0YvRjJGNUY4RjtGPUY/RkFGREZGRkhGSkZMRk5GUEZSRlRGV0ZZRmVuRmduRmpuLUkjbW5HRiQ2OVEjMjVGJ0YvRjJGNS9GOUY3RjtGPUY/RkFGREZGRkhGSkZMRk5GUEZSRlRGV0ZZRmVuL0ZoblEnbm9ybWFsRidGam4vJTFzdXBlcnNjcmlwdHNoaWZ0R1EiMEYnRitGKw== * X is computed as follows:
Xcvec[1]*AEvals[1]^25*AEbasis1[1]
+ Xcvec[2]*AEvals[2]^25*AEbasis2[1]
+ Xcvec[3]*AEvals[2]^25*AEbasis2[2];
Note that this calculation involved no matrix multiplications. What advantage might this have?
Maple will actually compute directly, so let's check.
A^25 . X;
Exercise:
8) Given the matrix B defined below, compute LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2OVEhRicvJSdmYW1pbHlHUTBUaW1lc35OZXd+Um9tYW5GJy8lJXNpemVHUSMxMkYnLyUlYm9sZEdRJmZhbHNlRicvJSdpdGFsaWNHUSV0cnVlRicvJSp1bmRlcmxpbmVHRjcvJSpzdWJzY3JpcHRHRjcvJSxzdXBlcnNjcmlwdEdGNy8lK2ZvcmVncm91bmRHUShbMCwwLDBdRicvJStiYWNrZ3JvdW5kR0ZDLyUnb3BhcXVlR0Y3LyUrZXhlY3V0YWJsZUdGNy8lKXJlYWRvbmx5R0Y3LyUpY29tcG9zZWRHRjcvJSpjb252ZXJ0ZWRHRjcvJStpbXNlbGVjdGVkR0Y3LyUscGxhY2Vob2xkZXJHRjcvJTBmb250X3N0eWxlX25hbWVHUSsyRH5Db21tZW50RicvJSptYXRoY29sb3JHRkMvJS9tYXRoYmFja2dyb3VuZEdGQy8lK2ZvbnRmYW1pbHlHRjEvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLyUpbWF0aHNpemVHRjQtRiM2JUYrLUklbXN1cEdGJDYlLUYsNjlRIkJGJ0YvRjJGNUY4RjtGPUY/RkFGREZGRkhGSkZMRk5GUEZSRlRGV0ZZRmVuRmduRmpuLUkjbW5HRiQ2OVEjMTNGJ0YvRjJGNS9GOUY3RjtGPUY/RkFGREZGRkhGSkZMRk5GUEZSRlRGV0ZZRmVuL0ZoblEnbm9ybWFsRidGam4vJTFzdXBlcnNjcmlwdHNoaWZ0R1EiMEYnRitGKw== * X (X as above) using eigenvalues and eigenvectors. Check by direct computation. (Note we cannot use this technique for powers of C, since we can't find a basis of eigenvectors of C.)
B := Matrix(3,3,[[0,1,1],[0,2,0],[-2,1,3]]);