A = [4,-1,-2;1,2,-3;0,0,4]
AtA = A'*A
eigv = eig(AtA)
normA = max(sqrt(eigv))
fprintf('For comparison: norm(A,2)=%g', norm(A,2))
