% ASSP - Students Project 
% Calculate DOS for Tight Binding
%
% Thomas Ganner, Klaus Lang
%-------------------------------------------------------------------------

function plotDOS


epsilon = [0,1,2,4];
t = [0.5,1,2,4];

% ... 1 -> SC, 2 -> FCC, 3 -> BCC
for lattice = 1 : 3
  for i=1:4
    switch(lattice)
      case 1 
        name = 'SC';         
      case 2 
        name = 'FCC';
      case 3 
        name = 'BCC';
    end

    Daten = load(['DOS_',name,'_Tight_Binding_1.txt']);
    dE  = Daten(:,1);
    DOS = Daten(:,2);

    figure(lattice)
    subplot(2,2,i)
    dE = epsilon(i)- t(i) .* dE;
    plot(dE,DOS,'r*');
    hold on;
    xlabel('E')
    ylabel('D(E)');
    ylim([0 1])
    title(['Density of States for ',name,' lattice: \epsilon = ',num2str(epsilon(i)),', t = ',num2str(t(i))]);
  end
end