  diary make_table.dia

  vn0 = 2.^[5 : 10] - 1;
  venorminf = [  3.2189644400852924e-03 ...
                 8.0357767939931435e-04 ...
                 2.0082180979019704e-04 ...
                 5.0200916270171092e-05 ...
                 1.2549946873496864e-05 ...
                 3.1374743667544891e-06];

  vh = 1./(vn0+1);
  vratio = zeros(size(venorminf));
  vratio(2:end) = venorminf(1:end-1) ./ venorminf(2:end);
  vlog2ratio = zeros(size(venorminf));
  vlog2ratio(2:end) = log10(vratio(2:end)) / log10(2.0);
  vconst = venorminf ./ vh.^2;

  % ... make table row by row with explicit formatting:
  lenvn0 = length(vn0);
  fprintf(['     N           h       enorm       ratio ', ...
                   '  log2ratio   enorm/h^2\n']);
  for i = 1 : lenvn0
    fprintf('%6d %11.4e %11.4e %11.4f %11.4f %11.4f\n', ...
      vn0(i), vh(i), venorminf(i), vratio(i), vlog2ratio(i), vconst(i));
  end;

  diary off;

