%RUN_UCONJ_RESTR_LAG

%Checks the restricted Stenger U-conjecture for a generalized
%Laguerre weight function

tic
global a ab ableg
ableg=r_jacobi01(30);
for a=0
  ab=r_laguerre(30,a);
  for n=30
    if n>30, error('n exceeds 30'); end
    dig=32; [U,d]=Uconj_restr_lag(dig,n); 
    d
    for k=1:n
      if imag(d(k))==0
        d(k)=d(k)+(1e-20)*i; break
      end
    end
  end
%  plot(d,'k*','MarkerSize',16);set(gca,'FontSize',24) 
  semilogx(d,'k*','MarkerSize',16);set(gca,'FontSize',24)
  axis([.6 1.5 -.45 .45])
  for i=1:n
    if real(d(i))<=0, error('conjecture is false'); end
  end
  disp('conjecture is true')
end
