global c alpha beta
f0='%8.0f %6.3f %6.0f\n';
disp('      n     c     Ncap')
eps0=.5e-12;
for n=[10 55 100 200 300]
  for c=[.1 .5 .9 .99 .999]
    Ncapmax=0;
    for alpha=[-.9:.2:-.1 0 .1:.2:.9 1.5 3 6]
      for beta=[-.9:.2:-.1 0 .1:.2:.9 1.5 3 6]
        if beta>=alpha
          [ab,Ncap]=r_subjacobi(n,eps0,c,alpha,beta);
%          [ab,Ncap]=r_subjacobi0(n,eps0,c,alpha,beta);
          if Ncap>Ncapmax, Ncapmax=Ncap; end
        end
      end
    end
    fprintf(f0,n,c,Ncapmax)
  end
  fprintf('\n')
end
