function r8vec_step_test ( )

%*****************************************************************************80
%
%% r8vec_step_test() tests r8vec_step().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    03 June 2012
%
%  Author:
%
%    John Burkardt
%
  n = 21;

  fprintf ( 1, '\n' );
  fprintf ( 1, 'R8VEC_STEP_TEST\n' );
  fprintf ( 1, '  R8VEC_STEP evaluates a step function;\n' );

  a = 10.0;
  b = 20.0;
  x = r8vec_linspace ( n, a, b );

  x0 = 17.25;
  f = r8vec_step ( x0, n, x );

  fprintf ( 1, '\n' );
  fprintf ( 1, '  Step function jump is at X=%g\n', x0 );

  r8vec2_print ( x, f, '  Table of step function values:' );

  return
end
