schur(mat, details)

Schur decomposes matrix A to QR, where Q is unitary, and R is triangular.

matThe matrix to be decomposed.
detailsIf Q is to be calculated and returned. (Defaults to false, otherwise the return type is structure.)
return valueThe triangular R, or the structure {r:R, q:Q}, depending on the parameter 'details'.

a = [444,42,827;804,8,12]'; s = qr(a,true) s.q * s.q' s.q * s.r