Back to Gallery
  • Remix
  • Share
  • New Entry
on 9 Oct 2021
  • 6
  • 49
  • 0
  • 0
  • 216
m=4;
A=companion(4+randi(m));
I=eye(4+randi(m));
B=companion(4+randi(m));
K=kron(A,kron(I,B));
e=eig(K);
e=e-round(e);
plot(e,'.')
axis square
function C=companion(p)
a=1;
for j=1:p
a=conv(a,[1 0 -1]);
end
m=2*p;
C=[-a(2:m+1); eye(m-1,m)];
end
Remix Tree