function [v] = normalize(u)
v = u / sqrt(sum(sum(u.^2)));
