No BSD License
a = rand(3) [u w v]=svd(a) [s m] = eig(a*a') [s1 m1] = eig(a'*a) g = sort(diag(w.*w))
Contact us at files@mathworks.com