---

w_freq = 0:0.001:2; % уменьшаем частоту, чтобы степени не улетали в миллионы
p_freq = 1j * w_freq;
GM = polyval(den_zam, p_freq);

figure;
plot(real(GM), imag(GM), 'LineWidth', 1.5);
grid on;
xlim([-500 500]); % ограничиваем масштаб по оси Re
ylim([-500 500]); % ограничиваем масштаб по оси Im
xlabel('Re(D)'); ylabel('Im(D)');
title('Годограф Михайлова');