fprintf('\n');
f0='%4.0f %5.1f %8.4f %23.12e %5.0f %8.2f\n';
disp('  n   bet     gam              I               N      can')
dig=40;
eps0=.5e-12;
%for n=[0 6 12 20 60 100]
for n=[1 7 13 21 61 101]
  for bet=[.1 .5 1 5 10 25 50 100]
%  for bet=0
    tic
%    for gam=[1.1 1.5 2 2.5 5 7.5 10]
%    for gam=[1 .5 .1 .05 .01 .001]
    for gam=[.5 .25 .1 .05 .01 .001 .0001]
      [I,N,can]=integral(dig,n,bet,gam,eps0); 
      fprintf(f0,n,bet,gam,I,N,can)
    end
%    fprintf('\n')
    toc
  end
  fprintf('\n')
end
