function [Q, R] = QRGivens(A)

% Use Givens rotations to compute a QR factorization for A

return;

end