 
 
 
qr takes as argument a numeric
square matrix A of size n.
qr factorizes numerically 
this matrix as Q*R where
Q is an orthogonal matrix (tQ*Q=I) and R is an upper triangular 
matrix. 
qr(A) returns only R, run Q=A*inv(R) to get Q.
Input :
Output is the matrix R :
Input :
Output is the matrix R :
 
 
