function [Q, R] = QRHouseholder(A)

% Use Householder reflections to compute a QR factorization for A

return;

end