function [ k ] = kvals( a,N)
%liefert zufaellige k-Werte
%a...Gitterkonstante
%N...Anzahl der Punkte
kmax = 2*pi/a;
k=kmax*rand(3,N);
end

