For any m×n matrix A, with linearly independent columns: A=QR Q is m x n its columns are an orthonormal basis for Col A R is n x n upper triangular positive diagonal ∣∣rj∣∣=∣∣aj∣∣ where r and a are columns of R and A