function r8bto_sl_test ( )

%*****************************************************************************80
%
%% R8BTO_SL_TEST tests R8BTO_SL.
%
%  Licensing:
%
%    This code is distributed under the MIT license. 
%
%  Modified:
%
%    08 July 2016
%
%  Author:
%
%    John Burkardt
%
  m = 2;
  l = 3;
  n = m * l;
  p = 2 * l - 1;
%
%  Hey, this is how MATLAB itself tells me you have to initialize
%  a 3D array...
%
  a        = [ 9.0, 1.0; 2.0, 8.0 ];
  a(:,:,2) = [ 3.0, 4.0; 6.0, 6.0 ];
  a(:,:,3) = [ 5.0, 6.0; 7.0, 7.0 ];
  a(:,:,4) = [ 7.0, 8.0; 8.0, 8.0 ];
  a(:,:,5) = [ 9.0, 0.0; 9.0, 9.0 ];

  fprintf ( 1, '\n' );
  fprintf ( 1, 'R8BTO_SL_TEST\n' );
  fprintf ( 1, '  R8BTO_SL solves a block Toeplitz system.\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  Block order M =  %d\n', m );
  fprintf ( 1, '  Block number L = %d\n', l );
  fprintf ( 1, '  Matrix order N = %d\n', n );

  r8bto_print ( m, l, a, '  The block Toeplitz matrix:' )
%
%  Set the desired solution.
%
  x = r8vec_indicator1 ( n );
%
%  Compute the right hand side.
%
  b = r8bto_mv ( m, l, a, x );

  r8vec_print ( n, b, '  The right hand side B:' );

  x = r8bto_sl ( m, l, a, b );

  r8vec_print ( n, x, '  The computed solution X:' );

  return
end
