all (abs (svd ([1, 2; 2, 1]) - [3; 1]) < sqrt (eps))
