% plot IDOS
% N_norm=1;
% figure
% plot(idos(1,:),idos(2,:)/N_norm,'b')

% plot smoothed DOS
figure
plot(ens,sum(map(:,:),2),'b')