%RUN_UCONJ_TWOINT

tic
global xi p q  ab0 ab1
%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(38+p+q); ab1=r_twoint(75,xi,p,q);
      for n=5
%      for n=1:30
        dig=16; [U,d1]=Uconj_twoint(dig,n); errd=1;
        while errd > .5e-6
          dig=dig+4;
          d0=d1; d00=d0;
          for j=1:n
            if imag(d00(j))<0, d00(j)=conj(d00(j)); end
          end
          [U,d1]=Uconj_twoint(dig,n); d11=d1;
          for j=1:n
            if imag(d11(j))<0, d11(j)=conj(d11(j)); end
          end
          err=abs((d00-d11)./d11);
          errd=max(err);
        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([0 .02 -.012 .012])
        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
