function normal_comparison ( )

  fprintf ( 1, '\n' );
  fprintf ( 1, '           N  RANDN  sum(RAND)\n' );
  fprintf ( 1, '\n' );

  n = 1;
  for logn = 0 : 6
    tic;
    r1 = randn ( n, 1 );
    t1 = toc;
    tic;
    r2 = sum ( rand ( n, 10 ), 2 );
    t2 = toc;
    fprintf ( 1, '  %10d  %f  %f\n', n, t1, t2 );
    n = n * 10;
  end

  fprintf ( 1, '\n' );
  fprintf ( 1, '  Repeat, but compute 1 by 1.\n' );
  fprintf ( 1, '\n' );

  n = 1;
  for logn = 0 : 6
    r = zeros(n,1);
    tic;
    for i = 1 : n
      r(i) = randn ( 1, 1 );
    end
    t1 = toc;
    tic
    for i = 1 : n
      r(i) = sum ( rand ( 1, 10 ), 2 );
    end
    t2 = toc;  
    fprintf ( 1, '  %10d  %f  %f\n', n, t1, t2 );  
    n = n * 10;
  end

  return
end
