from sympy.abc import a, b, c, d, e, f, s
from sympy.physics.control.lti import TransferFunction
from spb import *

tf1 = TransferFunction(a*s**2 + b*s + c, d*s**4 + e*s**3 + f*s**2, s)
plot_bode(
    tf1, initial_exp=-2, final_exp=2,
    params={
        a: (0.5, -10, 10),
        b: (0.1, -1, 1),
        c: (8, -10, 10),
        d: (10, -10, 10),
        e: (0.1, -1, 1),
        f: (1, -10, 10),
    },
    imodule="panel", ncols=3
)