gimse

import numpy as np
def qr(a):
    m=np.shape(a)[0]
    n=np.shape(a)[1]
    r=a.copy()

    ident=np.identity(m)
    q=np.identity(m)
    for i in range(0,n):
        x = r[i:, i].copy()