from sympy import symbols
from spb import *
a, s = symbols("a, s")
G = (s**2 + a) / (s**3 + 2*s**2 + 3*s + 4)
params={a: (-0.5, -4, 4)}
plot_root_locus(G, params=params, xlim=(-4, 1))