m = 5;
n = 3;
a = [ 1, 1, 1; 1, 2, 4; 1, 3, 9; 1, 4, 16; 1, 5, 25 ];
b = [ 1.0; 2.3; 4.6; 3.1; 1.2];
lda = m;

tol = eps / max ( max ( abs ( a(1:m,1:n) ) ) ); 
%
%  Factor the matrix.
%
[ kr, jpvt, qraux ] = dqrank ( a, lda, m, n, tol )

