[Snippet]Python SVD

a = numpy.asarray([[1,2,3],[4,5,6],[7,8,9]])
U,S,VT = numpy.linalg.svd(a)
print U.dot(numpy.diag(S)).dot(VT)
S[2] = 0
print U.dot(numpy.diag(S).dot(VT))