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