qr( A, Q, R)  computes the QR decomposition of A.  For an m-by-n matrix A,
	Q*R = A, where Q is m-by-m orthogonal and R is m-by-n upper-triangular.