%RUN_UCONJ_RESTR_BLOCKHAHN

%Checks the restricted Stenger conjecture for the (N+1)-block Hahn
%weight function w(x;N+1,0,0) 

tic
global ab ab0 
a=0; b=0;
for N=10
  dig=32;
  ab=loadvpa('coeffblockhahn10_0_0.txt',9,32,100,2);
  ab=double(ab);
  ab0=r_jacobi01(100);
  for n=25
    n
    [U,d,xn]=Uconj_restr_blockhahn(dig,n); 
    d
    plot(d,'k*','MarkerSize',16);set(gca,'FontSize',24) 
%    semilogx(d,'k*','MarkerSize',16);set(gca,'FontSize',24)
    axis([-.2 1.1 -.9 .9])
    for i=1:n
      conject=1;
      if real(d(i))<=0 
        disp('conjecture is false');
        [N n], d
        conject=0;
        break
      end
    end
%    if conject==0, break; end
  end
end
toc
