function [varargout] = qr(varargin)
%QR     Orthogonal-triangular decomposition.
%   [Q,R] = QR(A), where A is m-by-n, produces an m-by-n upper triangular
%   matrix R and an m-by-m unitary matrix Q so that A = Q*R.
%
%   [Q,R] = QR(A,0) produces the "economy size" decomposition.
%   If m>n, only the first n columns of Q and the first n rows of R are
%   computed. If m<=n, this is the same as [Q,R] = QR(A).
%
%   If A is full:
%
%   [Q,R,E] = QR(A) produces unitary Q, upper triangular R and a
%   permutation matrix E so that A*E = Q*R. The column permutation E is
%   chosen so that ABS(DIAG(R)) is decreasing.
%
%   [Q,R,E] = QR(A,0) produces an "economy size" decomposition in which E
%   is a permutation vector, so that A(:,E) = Q*R.
%
%   X = QR(A) and X = QR(A,0) return the output of LAPACK's *GEQRF routine.
%   TRIU(X) is the upper triangular factor R.
%
%   If A is sparse:
%
%   R = QR(A) computes a "Q-less QR decomposition" and returns the upper
%   triangular factor R. Note that R = CHOL(A'*A). Since Q is often nearly
%   full, this is preferred to [Q,R] = QR(A).
%
%   R = QR(A,0) produces "economy size" R. If m>n, R has only n rows. If
%   m<=n, this is the same as R = QR(A).
%
%   [C,R] = QR(A,B), where B has as many rows as A, returns C = Q'*B.
%   The least-squares solution to A*X = B is X = R\C.
%
%   [C,R] = QR(A,B,0) produces "economy size" results. If m>n, C and R have
%   only n rows. If m<=n, this is the same as [C,R] = QR(A,B).
%
%   Example: The least squares approximate solution to A*x = b can be found
%   with the Q-less QR decomposition and one step of iterative refinement:
%
%         if issparse(A), R = qr(A); else R = triu(qr(A)); end
%         x = R\(R'\(A'*b));
%         r = b - A*x;
%         e = R\(R'\(A'*r));
%         x = x + e;
%
%   See also LU, NULL, ORTH, QRDELETE, QRINSERT, QRUPDATE.

%   Copyright 1984-2004 The MathWorks, Inc.
%   $Revision: 5.13.4.3 $  $Date: 2004/06/25 18:52:14 $
%   Built-in function.

if nargout == 0
  builtin('qr', varargin{:});
else
  [varargout{1:nargout}] = builtin('qr', varargin{:});
end