function [z] = mcqr(A, y)
    [m, n] = size(A);
    [Q, R] = qr(A3);
    e = Q' * y
    Rt = R(1:n, :)
    z = solsup(Rt, e(1:n,:) 
endfunction