a = 4.49; % Länge der konventionellen Einheitszelle
C = 200; % Federkonstante in N/m
m = 0.344*10^-24; % Masse in kg

kvaln = 1000000; %10^6 Werte für die Dispersionsrelation
%Berechnung dauert sehr lange

% Position des betrachteten Atoms und seiner Nachbarn
% FCC:
an = a/2 * [0,0,0; 1,1,0; 1,0,1; 0,1,1; -1,-1,0; -1,0,-1; 0,-1,-1; 1,-1,0; 1,0,-1; 0,1,-1; -1,1,0; -1,0,1; 0,-1,1];


% BCC:
% an = a/2 * [0,0,0; 1,1,1; 1,-1,1; 1,-1,-1; 1,1,-1; -1,1,1; -1,-1,1; -1,-1,-1; -1,1,-1];
% SC:
% an = a/2 * [0,0,0; 1,0,0; -1,0,0; 0,1,0; 0,-1,0; 0,0,1; 0,0,-1];


k = kvals(a,kvaln); %ruft die functino kvals auf, welche zufällige k-Werte liefert
%a...Gitterkonstante
%N=kvaln...Anzahl der Punkte

w = NaN(3,size(k,2)); %erzeugen des Vektors für die berechneten Werte

for i = 1:size(k,2)
    disp(i); %displays the current number
    k_aktuell = [k(1,i) k(2,i) k(3,i)]; %Erstellung des betrachteten k-Vektors
    w(:,i) = disp_solver(an, m, C, k_aktuell);
end

