load('Ed.mat')
nk = 500;
k = linspace(1,length(Ed_line),nk);
Ed = interp1(1:length(Ed_line),Ed_line,k,'spline');
nw = 1000;
w = linspace(-1,0.2,nw);
Gammar1 = 0.005;
Gammar2 = 0.005;
delta = 0.01;
Ak_d = zeros(nw,nk);
for i=1:nw
for j=1:nk
sigma_d = -1i*Gammar1+delta*delta/(w(i)+Ed(j)+1i*Gammar2);
Ak_d(i,j) = imag(sigma_d)/((w(i)-Ed(j)-real(sigma_d))^2+imag(sigma_d)^2);
end
end
[X,Y] = meshgrid(k,w);
mesh(X,Y,-Ak_d/pi);