G = zpk([-1], [j, -j], 10)

%

nyquist(G)

%

margin(G)

%

bode(G, {10, 100})

%

mag2db(evalfr(G, 0))

%

bode(G, {1.5, 5})