GH = tf(1, [1, 2, 5, 6])

%

bode(GH)

%

CL = feedback(GH, 1)

%

pole(CL)

%

margin(GH)

%

nyquist(GH)

%

nyquist(3.5 * GH)

%

bode(3.5 * GH)

%

margin(5 * GH)

%

nyquist(5 * GH)