אולי
נעזרתי במטלב כדי לפתור את המטריצה. מצאתי ערכים עצמיים: 0, 9 ו- 9. בגלל סימטריות המטריצה ניתן לקבל גם פירוק אורתוגונלי עבור הוקטורים העצמיים:
MN = UEU' = sum(k=1..3)[e(k)*u(k)*u'(k)] = sum(k=1..2)[9*u(k)*u'(k)] = (3*[u1 u2])*(3*[u1 u2])'
כאשר E מטריצה אלכסונית של הערכים העצמיים, e מסמל ערך עצמי, U מטריצת וקטורים עצמיים. u וקטור עצמי (עמודה). ו- ' מסמן את פעולת ה- transpose. כלומר מצאתי את M ו- N (עד כדי נרמול בקבוע):
==> NM = 9*[u1 u2]'*[u1 u2] = 9*I = [9 0; 0 9]