As written, SVD::solve, in the overloaded version for matrix right-hand sides b (book page 71), requires (and tests) that the matrix A be square. This is not exactly a bug, since the subsection is named "SVD of a Square Matrix". However, it is clear in later subsections that the intent was that this routine also work for non-square A, which it doesn't.

User take2rec has kindly contributed a fixed version, at http://www.nr.com/forum/showthread.php?t=1355.

We will use this fixed version in reprintings of the book and CDROM.