hold on;
l=length(kvec);
for j=1:length(F)
    plot(kvec,ones(l,1)*F(j));
end