function y = normalize(f)

y = f/(norm(f, 'fro'));
