RQ decomposition
[R, Q]= sip_rq(X)
X
Function sip_rq
operates similarly to the qr function, but
in reverse. It produces an upper triangular matrix R
of the
same dimensions as X
and an orthogonal matrix
Q
so that X = Q*R
http://siptoolbox.sourceforge.net