%% TEST10 tests R8_NORMAL_01.
%
%  Modified:
%
%    28 June 2006
%
%  Author:
%
%    John Burkardt
%
  clear

  seed = 123456789;

  fprintf ( 1, '\n' );
  fprintf ( 1, 'TEST08\n' );
  fprintf ( 1, '  R8_NORMAL_01 computes a sequence of \n' );
  fprintf ( 1, '  normally distributed pseudorandom numbers.\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  Verify that we can change the seed\n' );
  fprintf ( 1, '  and get the desired results.\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  The initial seed is %d\n', seed );
  fprintf ( 1, '\n' );

  seed = seed_init;

  fprintf ( 1, '\n' );
  fprintf ( 1, '         I    Seed(in)   Seed(out)   R8_NORMAL_01\n' );
  fprintf ( 1, '\n' );

  for ( i = 1; i <= 5; i++ )
  {
    seed_input = seed;
    [ value, seed ] = r8_normal_01 ( seed );
    seed_output = seed;

    fprintf ( 1, '  %8d  %12d  %12d  %14f\n', ...
      i, seed_input, seed_output, value );
  }

  seed = seed_init;

  fprintf ( 1, '\n' );
  fprintf ( 1, '  Resetting seed to repeat, after an ODD number of steps.\n' );
  fprintf ( 1, '\n' );

  for ( i = 6; i <= 10; i++ )
  {
    seed_input = seed;
    [ value, seed ] = r8_normal_01 ( seed );
    seed_output = seed;

    fprintf ( 1, '  %8d  %12d  %12d  %14f\n', ...
      i, seed_input, seed_output, value );
  }

  seed = seed_init;

  fprintf ( 1, '\n' );
  fprintf ( 1, '  Resetting seed to repeat, after an EVEN number of steps.\n' );
  fprintf ( 1, '\n' );

  for ( i = 11; i <= 15; i++ )
  {
    seed_input = seed;
    [ value, seed ] = r8_normal_01 ( seed );
    seed_output = seed;

    fprintf ( 1, '  %8d  %12d  %12d  %14f\n', ...
      i, seed_input, seed_output, value );
  }