% clf

set(0,'defaultaxesfontsize',18);
x=[0.1:0.1:4.0];
B=TNVandBD(x);
A=TNExpand(B); A=A(end:-1:1,:);
z=SREigenvalues(B);
y=eig(A);
y=sort(abs(y)); y=y(end:-1:1);
z=abs(z);
semilogy(z,'o','linewidth',1,'markersize',9); hold on
semilogy(y,'r+','linewidth',1,'markersize',9)
semilogy(1:length(x),z(1)/10^16*ones(1,length(x)),'k--','linewidth',1,'markersize',12);
axis([0,length(x)+1,z(end)/100,z(1)*100]);
xlabel('Eigenvalue index');
ylabel('|\lambda_i|');