Returns solution of linear system matrix * solution = b-vector. Employs the
singular value decomposition method. See the discussion in Numerical Recipes in
C, section 2.6, especially as to the semantics of `threshold.'