% INEQ Conjectured inequality for all zeros of Jacobi 
% polynomials P_n^{a,b}
%
function [k0,r0]=ineq(n,a,b)
ab=r_jacobi(n+1,a,b);
k0=0; r0=0;
for k=1:n
%
% Note that in gauss.m the zeros are ordered increasingly.
%
  xw=gauss(k,ab); xw1=gauss(k+1,ab);
%  for r=k
  for r=1:k
    th=acos(xw(k+1-r,1)); th1=acos(xw1(k+2-r,1));
    if k*th>=(k+1)*th1 
      k0=k; r0=r;
      break
    end
  end
  if k0>0
    break
  end
end