QR decomp finds R correctly (i.e. matches octave), solves square system; u_k (which will make Q) don't match octave