clc
clear
n=8;
A=rand(n,n);
[Q,R]=qrgv(A); % calling for the QR givens function
test_norm = norm(A-Q*R)
matrixplot(Q*Q');