tf = TransferFunction(1, s**3 + 2*s**2 + s, s)
p1 = plot_bode_phase(
    tf, unwrap=False, show=False, title="unwrap=False")
p2 = plot_bode_phase(
    tf, unwrap=True, show=False, title="unwrap=True")
plotgrid(p1, p2)