%RUN_SUCONJ_RESTR_JAC

tic
global sa sb  sab sableg
dig=32; digits(dig);
sa=vpa(51/50); sb=vpa(51/50); N=45;
M=ceil(N/2);
sab=sr_jacobi(dig,N,sa,sb); sableg=sr_jacobi(dig,M);
for n=30
  [U,d]=sUconj_restr_jac(dig,n); 
  d=double(d)
  plot(d,'k*','MarkerSize',16);set(gca,'FontSize',24) 
%  semilogx(d,'k*','MarkerSize',16);set(gca,'FontSize',24)
%  axis([0 .055 -.03 .03])
  for i=1:n
    if real(d(i))<=0, error('conjecture is false'); end
  end
%  n
  disp('conjecture is true')
end
toc
