
A=[.70000, .70711; .70001, .70711]
[Q,R]=cgs(A);
CGS=norm(Q'*Q-eye(2))

[Q,R]=mgs(A);
MGS=norm(Q'*Q-eye(2))

[Q,R]=qr(A);
HH=norm(Q'*Q-eye(2))

