NEDOS   = 10
HalfBW = 1
dE =HalfBW/NEDOS

a = {0}
b = {}

for i=1,NEDOS do
  a[#a+1] = (i-0.5) * dE -- energy axis from 0 to 1
  b[#b+1] = 0.5 / NEDOS -- flat density of states
end

G0 = ResponseFunction.New( {a,b,mu=0,type="ListOfPoles", name="G0"} )
G0_inv = ResponseFunction.InvertEnergy(G0)

print(G0)
print(G0_inv)