function ki=qrit(A,n)

ki=A;
for i=1:n
    [Q,R]=qr(ki);
    ki=R*Q;
end
ki=Q*R;