function ft = f(at)
    global kvec c epsilon
    ft=(-kvec.*kvec*epsilon+1i*kvec*c).*at;
end

