[v,w] = eig(D) v = -0.2801 0.7408 -0.6106 -0.9537 -0.2873 0.0889 -0.1096 0.6072 0.7869 w = 0.0913 0 0 0 0.2266 0 0 0 0.4178 [x,y,z] = svd(D) x = -0.6106 -0.7408 -0.2801 0.0889 0.2873 -0.9537 0.7869 -0.6072 -0.1096 y = 0.4178 0 0 0 0.2266 0 0 0 0.0913 z = -0.6106 -0.7408 -0.2801 0.0889 0.2873 -0.9537 0.7869 -0.6072 -0.1096 v*w/v ans = 0.2873 -0.0465 -0.0960 -0.0465 0.1050 -0.0008 -0.0960 -0.0008 0.3434 x*y*z' ans = 0.2873 -0.0465 -0.0960 -0.0465 0.1050 -0.0008 -0.0960 -0.0008 0.3434