gram schmidt numpy code example

Example: python qr decomposition

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