function [Q,R] = QRDecomp(A)
%QRDECOMP QR factorization of A.
%   [Q,R]= QRDECOMP(A) returns the QR factorization of A.
    [Q,R] = qr(A);
