%RUN_UCONJ_EXT_TWOINT

tic
global xi p q ab ab0 
for xi=1/2
%for xi=1/10:2/10:9/10
  xi
  for p=0 
    for q=0
%  for p=0:5
%    for q=0:5
      ab0=r_jacobi(16+p+q); ab=r_twoint(30,xi,p,q);
      for n=15
        dig=16; [U,d1]=Uconj_ext_twoint(dig,n); err=1;
        while err > .5e-4
          dig=dig+4;
          d0=d1; 
          [U,d1]=Uconj_ext_twoint(dig,n); 
          err=errd(n,d1,d0);
        end
%        dig
        d=d1;
        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([5*10^(-5) 2*10^(-1) -.05 .05])
        for i=1:n
          if real(d(i))<=0 
            disp('conjecture is false')
            [xi p q n]
          end
        end
      end
    end
  end
end
disp('conjecture is true')
toc
