 function hardcase(B,delta)
 %function hardcase(B,delta)
% generate a hard case with plots - singleton smallest eig of B is
% assumed
[V,D]=eig(B);
n=size(B,2);
g=V(:,round(n/2));
lambdas=plotrs(B,g,delta);
