Hi Robert! Oh! Its not a homework problem... I read the Golub book, it tells me what an orthogonal complement is, however, I cannot understand how I can code it. I know about svd from numpy's mlab, but I what I want to know is how can I compute an orthogonal complement, using SVD or otherwise. Thanks for the prompt reply :)