%
% Make plots of the eigenvalues of the matrices in question
%
close all
clear all
nilmats
e = zeros(7*16,1);
for k=1:16
  e(1+7*(k-1):7+7*(k-1)) = eig(nilg{k});
end
figure(1)
plot( real(e), imag(e), 'k.')
axis([-0.01,0.01,-0.01,0.01])
axis('square')
hold on
badmats
e5 = zeros(7*32,1);
for k=1:32
  e5(1+7*(k-1):7+7*(k-1)) = eig(fiveg{k});
end
plot( real(e5), imag(e5), 'k.' )
e3 = zeros(7*368,1);
thirdOrderZeros
for k=1:368
  e3(1+7*(k-1):7+7*(k-1)) = eig(threeg{k});
end
figure(2)
plot( real(e3), imag(e3), 'k.')
small = 2.0e-5;
axis( small*[-1,1,-1,1] )
axis('square')