tf1 = TransferFunction(1, s**2 + 2*s + 1, s)
tf2 = TransferFunction(1, s**2 - 2*s + 1, s)
plot_nichols(tf1, tf2, xlim=(-360, 360))