%RUN_VCONJ_RESTR_LOG1

%Checks the restricted Stenger V-conjecture for an algebraic/
%logarithmic weight function

tic
global ab ableg
%abal=loadvpa('coeff_lalg0log1.txt',9,32,200,2);
%abal=loadvpa('coeff_lalgm1halflog1.txt',9,32,200,2);
%abal=loadvpa('coeff_lalg1halflog1.txt',9,32,200,2);
%abal=loadvpa('coeff_lalg1log1.txt',9,32,100,2);
%abal=loadvpa('coeff_lalg2log1.txt',9,32,100,2);
abal=loadvpa('coeff_lalg5log1.txt',9,32,100,2);
ab=double(abal(1:30,:));
ableg=r_jacobi01(30);
for n=30
  if n>30, error('n exceeds 30'); end
  dig=32; [V,d]=Vconj_restr_log1(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([2*10^(-6) 2*10^(-1) -.03 .03])
  for i=1:n
    if real(d(i))<=0, error('conjecture is false'); end
  end
  disp('conjecture is true')
end
