%RUN_UCONJ_EXT_HERM

tic
global mu ab sabhrh abjac
dig=32;
sabhrh=loadvpa('coeffhrghermite0.txt',9,32,100,2);
for mu=0:1/2:25
  mu
  ab=r_hermite(30,mu);
  abjac=r_jacobi01(100,0,2*mu);
  for n=15
    [U,d]=Uconj_ext_herm(dig,n);
%    d
    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([10^(-3) 10^(0) -.25 .25])
    for i=1:n
      if real(d(i))<=0, error('conjecture is false'); end
    end
    disp('conjecture is true')
  end
end
toc
