function vdc_sequence_test ( )

%*****************************************************************************80
%
%% vdc_sequence_test() tests vdc_sequence().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    09 August 2016
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'vdc_sequence_test():\n' );
  fprintf ( 1, '  VDC_SEQUENCE returns elements I1 through I2 of \n' );
  fprintf ( 1, '  a van der Corput sequence.\n' );

  i1 = 7;
  i2 = 7;
  n = abs ( i2 - i1 ) + 1;
  r = vdc_sequence ( i1, i2 );
  fprintf ( 1, '\n' );
  r8vec_transpose_print ( n, r, '  R=VDC_SEQUENCE(  7,  7):' );

  i1 = 0;
  i2 = 8;
  n = abs ( i2 - i1 ) + 1;
  r = vdc_sequence ( i1, i2 );
  fprintf ( 1, '\n' );
  r8vec_transpose_print ( n, r, '  R=VDC_SEQUENCE(  0,  8):' );

  i1 = 8;
  i2 = 0;
  n = abs ( i2 - i1 ) + 1;
  r = vdc_sequence ( i1, i2 );
  fprintf ( 1, '\n' );
  r8vec_transpose_print ( n, r, '  R=VDC_SEQUENCE(  8,  0):' );

  i1 = -3;
  i2 = +5;
  n = abs ( i2 - i1 ) + 1;
  r = vdc_sequence ( i1, i2 );
  fprintf ( 1, '\n' );
  r8vec_transpose_print ( n, r, '  R=VDC_SEQUENCE( -3,  5):' );

  i1 = 100;
  i2 = 105;
  n = abs ( i2 - i1 ) + 1;
  r = vdc_sequence ( i1, i2 );
  fprintf ( 1, '\n' );
  r8vec_transpose_print ( n, r, '  R=VDC_SEQUENCE(100,105):' );

  return
end
