from sympy.abc import s
from sympy.physics.control.lti import TransferFunction
from spb import plot_nichols
tf = TransferFunction(50*s**2 - 20*s + 15, -10*s**2 + 40*s + 30, s)
plot_nichols(tf)