QRfact

Usage: QR(<matrix>)
QR Factorization of a MxN General Matrix. On exit, the input array is changed to the upper triangle and diagonal is the min(M,N) by N upper triangular matrix R. The lower triangle, together with the output vector, is the orthogonal matrix Q as a product of min(M,N) elementary reflectors.
Returns The scalar factors of the elementary reflectors.
Ex:QR([1,2;3,4])
Uses the LAPACK function geqrf. See:http://www.osl.iu.edu/research/mtl/doc/geqrfLapackMatAVectorT.html