array = np.array([19, 14, 16, 13, 15,15,12,8,9,18,12,5 ]).reshape(4,3) q,r = LA.qr(array) r # if diagonal exponent is zero, then v_i should be removed