from sympy import symbols, exp
s = symbols("s")
G1 = 1 / (s * (s + 1) * (s + 10))
G2 = G1 * exp(-5*s)
plot_bode(G1, G2, phase_unit="deg", n=1e04)