%RUN_UCONJ_EXT_BLOCKHAHN

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

tic
global N ab ab0 w
a=0; b=0;
for N=10
  w=zeros(N+1,1);
  p=1:N;
  w(1)=prod(b+p)/prod(p); w(N+1)=prod(a+p)/prod(p);
  for k=2:N
    p=vpa(1:k-1); q=vpa(1:N-k+1);
    w(k)=prod(a+p)*prod(b+q)/(prod(p)*prod(q));
  end
  dig=32;
  ab=loadvpa('coeffblockhahn10_0_0.txt',9,32,100,2);
  ab=double(ab);
  ab0=r_jacobi01(100);
  for n=2:30
    n
    [U,d,xn,flxn]=Uconj_ext_blockhahn(dig,n); 
%    d
%    plot(d,'k*','MarkerSize',16);set(gca,'FontSize',24) 
%    semilogx(d,'k*','MarkerSize',16);set(gca,'FontSize',24)
%    axis([-.3 1.2 -.8 .8])
    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
