RQ decomposition
[R, Q]= sip_rq(X)
XFunction 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