G = tf([1], [1, 11, 39, 29])

%

bode(G)

%

K = 1 / db2mag(-41)

%

margin(K * G)

%

Gcl = feedback(K * G, 1)

%

step(Gcl)

%

stepinfo(Gcl)

%

pzmap(K * G)