  average_old = 0.0;
  for logn = 1 : 6
    n = 10^logn;
    average = vector_code ( n );
    fprintf ( 1, '%7d  %14.6f  %14.6f\n', n, average, average - average_old );
    average_old = average;
  end

