The QR decomposition of a square matrix A is A=QR, where
Q is an orthogonal matrix (QTQ=I) and R is upper triangular.
The qr command finds the QR decomposition of a
matrix.
qr takes
A, a numeric square matrix.
qr(A) returns a list [Q,R] with Q and R from the QR
decomposition.