a1 = {0, -1,   0,   5,   6}
b1 = {  0.2, 0.2, 0.3, 0.3}
G1 = ResponseFunction.New( {a1,b1,mu=0,type="ListOfPoles", name="G1"} )


a1 = {0, -0.5, 5.75}
b1 = {    0.2,  0.2}
G2 = ResponseFunction.New( {a1,b1,mu=0,type="ListOfPoles", name="G2"} )

G3 = G1-G2
print(G3)