% Pure QR algorithm

function A = pure_qr(A,N)

for k=1:N
    
    [QM,RM] = mgs(A);
    A = RM*QM;
    
end