function normal_test12 ( )

%*****************************************************************************80
%
%% TEST12 tests R8_NORMAL_01 and R8VEC_NORMAL_01.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license. 
%
%  Modified:
%
%    17 July 2006
%
%  Author:
%
%    John Burkardt
%
  n = 10;
  seed_init = 123456789;

  fprintf ( 1, '\n' );
  fprintf ( 1, 'TEST12\n' );
  fprintf ( 1, '  R8_NORMAL_01 gets pseudeonormal values one at a time.\n' );
  fprintf ( 1, '  R8VEC_NORMAL_01 computes a vector of values.\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  For the same seed, the results should be identical,\n' );
  fprintf ( 1, '  but R8VEC_NORMAL_01 might be faster.\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  Initial seed is %d\n', seed_init );

  seed = seed_init;

  a = zeros ( n, 1 );

  for i = 1 : n
    [ a(i), seed ] = r8_normal_01 ( seed );
  end

  seed = seed_init;
  [ b, seed ] = r8vec_normal_01 ( n, seed );

  fprintf ( 1, '\n' );
  fprintf ( 1, '       I      A(I)          B(I)\n' );
  fprintf ( 1, '          (R8_NORMAL_01)  (R8VEC_NORMAL_01)\n' );
  fprintf ( 1, '\n' );

  for i = 1 : n
    fprintf ( 1, '  %6d  %12f  %12f\n', i, a(i), b(i) );
  end

  return
end