% Example 3.12 and 3.16
% Matrix A below is rank deficient
A = [ -1  1 0
      -1  0 1
       0 -1 1 ];
[Q, R] = qr(A),

[U, S, V] = svd(A),


