%RUN_UCONJ_RESTR_LOG2

%Checks the restricted Stenger U-conjecture for an algebraic/
%square-logarithmic weight function

tic
global ab ableg
%abal=loadvpa('coeff_lalg0log2.txt',9,32,200,2);
%abal=loadvpa('coeff_lalgm1halflog2.txt',9,32,100,2);
%abal=loadvpa('coeff_lalg1halflog2.txt',9,32,100,2);
%abal=loadvpa('coeff_lalg1log2.txt',9,32,100,2);
abal=loadvpa('coeff_lalg2log2.txt',9,32,100,2);
%abal=loadvpa('coeff_lalg5log2.txt',9,32,100,2);
ab=double(abal(1:30,:));
ableg=r_jacobi01(30);
for n=15
  if n>30, error('n exceeds 30'); end
  dig=32; [U,d]=Uconj_restr_log2(dig,n); 
  d
  for k=1:n
    if imag(d(k))==0
      d(k)=d(k)+(1e-20)*i; break
    end
  end
%  plot(d,'k*','MarkerSize',16);set(gca,'FontSize',24) 
  semilogx(d,'k*','MarkerSize',16);set(gca,'FontSize',24)
%  axis([10^(-3) 8*10^(-2) -.025 .025])
  for i=1:n
    if real(d(i))<=0, error('conjecture is false'); end
  end
  disp('conjecture is true')
end
