function rankdef
  format long
  A=randn(6,3);
  c=randn(3,1);
  A=[A A*c];
  b=randn(6,1);
%
  x1=A\b
  res1=norm(b-A*x1)
  disp(' ')
%
  [Q,R]=qr(A,0);
  x2=R\(Q'*b)
  res2=norm(b-A*x2)
  disp(' ')
%
  [U,S,V]=svd(A,0);
  x3=V*(S\(U'*b))
  res3=norm(b-A*x3)
  disp(' ')
%
  R=chol(A'*A);
  x4=R\(R'\(A'*b))
  res4=norm(b-A*x4)