---

w_full = 0:0.01:12; % увеличили частоту, чтобы дойти до 10.87
GM_full = polyval(den_zam, 1j*w_full);
figure; plot(real(GM_full), imag(GM_full), 'b', 'LineWidth', 1.5);
grid on; xlabel('Re'); ylabel('Im');
title('Полный годограф Михайлова');

w_small = 0:0.001:0.5;
GM_small = polyval(den_zam, 1j*w_small);
figure; plot(real(GM_small), imag(GM_small), 'r', 'LineWidth', 2);
grid on; xlabel('Re'); ylabel('Im');
title('Начало годографа: квадранты I -> II -> III');