function f = f1g(v,k)
v_0 = ifft(v);
v_x = ifft(i*k.*v);
v_c = ifft(conj(fliplr(v)));
f =  2*(v_x.^2).*v_c./(1+abs(v_0).^2);