import numpy as np
M = np.mat('4 11 14; 8 7 -2')
print(M)
U, Svd, V = np.linalg.svd(M, full_matrices=False)
print(U)
print(Svd)
print(V)
print(U * U.T)
print(V * V.T)
S = np.diag(Svd)
print(S)
print(U * S * V)
import numpy as np
M = np.mat('4 11 14; 8 7 -2')
print(M)
U, Svd, V = np.linalg.svd(M, full_matrices=False)
print(U)
print(Svd)
print(V)
print(U * U.T)
print(V * V.T)
S = np.diag(Svd)
print(S)
print(U * S * V)