T=1.5;
k=0.2;
a=[1 3 3 3+k*T k];
[x]=roots(a);
y1=real(x);
y2=imag(x);
my=min(abs(y2./y1));
eta=max(y1);
%далее корневой годограф
T=1.5;
r=0.2;
w=tf([0 0 0 3+k*T k],[1 3 3 3.3 0.2]);
rlocus(w); sgrid;
%Карта нулей и полюсов
w=tf([0 1 1 2 1],[1 3 3 3 1]);
[p]=pole(w);
[z]=zero(w);
pzmap(w);