Thanks Robert! I was using mlab's svd function, which returns an mxn matrix for u, and hence was unable to see how to compute the orthogonal complement! I realize that numpy's svd gives the mxm matrix! Thanks again. -Smriti