function ft = f(at)
    global kvec c epsilon
    u=fft(at);
    u2=u.*u/2;
    b=ifft(u2);
    ft=-kvec.*kvec*epsilon.*at+1i*kvec.*b;
end
