function f = f1g(v,k,x) v_0 = ifft(v); f = (-abs(v_0).^2-x.^2).*v_0;