U^2 + V^2 where U=cos(2*pi*f + p) - cos(p), V=cos(2*pi*f + p - 2*pi*f/N) - cos(p - 2*pi*f/N)