NEDOS   = 10
HalfBW = 1
dE =HalfBW/NEDOS

a0 = {0}
b0 = {}

for i=-NEDOS,NEDOS do
  a0[#a0+1] = i * dE -- energy axis from -1 to 1
  b0[#b0+1] = 0.5 -- flat density of states
end

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

print(ResponseFunction.ToTable(G0))
print(ResponseFunction.ToTable(G0_inv))