%RUN_CHECK_POS

global a b ab ableg
N=40; M=ceil((N+1)/2);
ab=r_jacobi(N,1,0); ableg=r_jacobi(M);
for n=2:40
  [min_int,min_anm1]=check_pos(n);
  if (min_int==0)|(min_anm1==0), error('zero'); end
  [n min_int min_anm1]
  fprintf('\n') 
end
