A =

     1     3     1     2
     0     0     0     0
     1     0     0     0
     1     3     1     2

[U S V]=svd(A)

U =

    0.7063   -0.0344    0.7071    0.0000
         0         0   -0.0000    1.0000
    0.0486    0.9988    0.0000    0.0000
    0.7063   -0.0344   -0.7071   -0.0000

S =

    5.4835         0         0         0
         0    0.9650         0         0
         0         0    0.0000         0
         0         0         0    0.0000

V =

    0.2665    0.9638         0         0
    0.7728   -0.2136    0.5976    0.0000
    0.2576   -0.0712   -0.3586    0.8944
    0.5152   -0.1424   -0.7171   -0.4472

V2 = V(:,3:4)

[Q R e]=qr(V2','vector')

Q =

   -0.3721    0.9282
    0.9282    0.3721


R =

    0.9636   -0.1482         0   -0.2224
         0   -0.8321         0    0.5547

e =

     3     4     1     2