a1 = {0, -1,-0.5, 0  }
b1 = {  0.1, 0.1, 0.1}
G1 = ResponseFunction.New( {a1,b1,mu=0,type="ListOfPoles", name="G1"} )

a2 = {0, 0.5,  1,  1.5}
b2 = {   0.1, 0.1, 0.1}
G2 = ResponseFunction.New( {a2,b2,mu=0,type="ListOfPoles", name="G2"} )

G = G1 + G2

print("The response function of G is")
print(G)
omega = 1.1
Gamma = 0.001
print()
print("Evaluated at omega =",omega," and Gamma =",Gamma," yields\nG (",G(omega,Gamma),") = G1 (",G1(omega,Gamma),") + G2(",G2(omega,Gamma),")")