Definitely learn how to do QR factorization on your calculator. If you also wish to see it done on a computer screen, try the following example sequence of commands at live.sympy.org: A=Matrix([[1,1,1],[1,1,3],[2,3,4]]) A Q,R=A.QRdecomposition() Q R