%RUN_UCONJ_RESTR_HAHN

%Checks the restricted Stenger conjecture for (N+1)-point
%discrete Hahn weight functions

tic
global ab ab0  
dig=32;
%for b=[[-.9:.6:.9] [1.7:.7:3.8] [4.7:.9:7.4]]
for b=0
%for a=[[-.9:.6:.9] [1.7:.7:3.8] [4.7:.9:7.4]]
  for a=b
    [a b]
    for N=30
      ab0=r_jacobi01(N);
      ab=r_hahn(N,a,b);
      for n=30
        [U,d,xn]=Uconj_restr_hahn(dig,n); 
%        xn
%        plot(d,'k*','MarkerSize',16);set(gca,'FontSize',24) 
%       semilogx(d,'k*','MarkerSize',16);set(gca,'FontSize',24)
%        axis([-.5 1.5 -.4 .4])
        for i=1:n
          if real(d(i))<=0 
            disp('conjecture is false');
            [N n], d 
            break
          end
        end
      end
    end
  end
end
toc
