using Plots , ControlSystemsBase , Plots.Measures
begin
num = [2000, 1000]
den = [1, 60, 500, 0]
pp = bodeplot(sys, xminorticks=true, yminorticks=true, layout=(1,2), size=(800,300), margins=5mm, label=false, xlabel="Frequency [rad/s]")
begin
numerator = [1e-2, 1e-4, 1e-2]
denominator = [1/64, 2e-2/8, 1, 0, 0]
gg = bodeplot(sys2, xminorticks=true, yminorticks=true, layout=(1,2), size=(800,300), margins=5mm, label=false, xlabel="Frequency [rad/s]")
end