%RUN_UCONJ_JAC

tic
global a b  ab0 ab00
a=0; b=0;
M=ceil((75+a)/2);
ab0=r_jacobi(75,a,b); ab00=r_jacobi(M,0,b);
for n=15
  dig=16; [U,d1]=Uconj_jac(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_jac(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 .08 -.05 .05])
  for i=1:n
    if real(d(i))<=0, error('conjecture is false'); end
  end
%  n
  disp('conjecture is true')
end
end
end
toc
