---

clc;
clear;
close all;

% Диапазон частот
w = linspace(0, 10, 2000);

% Комплексная переменная p = j*w
p = 1i*w;

% Характеристический полином
A = p.^4 + 3*p.^3 + 4*p.^2 + 2*p + 1;

% Действительная и мнимая части
Re = real(A);
Im = imag(A);

% Построение годографа Михайлова
figure;
plot(Re, Im, 'LineWidth', 2);
grid on;
axis equal;

xlabel('Re[A(j\omega)]');
ylabel('Im[A(j\omega)]');
title('Годограф Михайлова');

% Оси координат
hold on;
xline(0,'k--');
yline(0,'k--');

% Начальная и конечная точки
plot(Re(1), Im(1), 'go', 'MarkerSize', 8, 'LineWidth', 2);
plot(Re(end), Im(end), 'ro', 'MarkerSize', 8, 'LineWidth', 2);

legend('Годограф', 'Ось Y', 'Ось X', ...
       'Начало (\omega=0)', 'Конец');