function r8blt_sl_test ( )

%*****************************************************************************80
%
%% R8BLT_SL_TEST tests R8BLT_SL.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    06 October 2015
%
%  Author:
%
%    John Burkardt
%
  ml = 3;
  n = 10;
  seed = 123456789;

  fprintf ( 1, '\n' );
  fprintf ( 1, 'R8BLT_SL_TEST\n' );
  fprintf ( 1, '  R8BLT_SL solves A*x=b, where A is an R8BLT matrix;\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  Matrix order N = %d\n', n );
  fprintf ( 1, '  Lower bandwidth ML = %d\n', ml );

  [ a, seed ] = r8blt_random ( n, ml, seed );
  r8blt_print ( n, ml, a, '  The R8BLT matrix:' );
%
%  Set the desired solution.
%
  x = r8vec_indicator1 ( n );
  r8vec_print ( n, x, '  x:' );
%
%  Compute the corresponding right hand side.
%
  b = r8blt_mv ( n, ml, a, x );
  r8vec_print ( n, b, '  The right hand side:' );
%
%  Solve the linear system.
%
  x = r8blt_sl ( n, ml, a, b );
  r8vec_print ( n, x, '  Solution x:' );

  return
end
