function r8bto_mv_test ( )

%*****************************************************************************80
%
%% R8BTO_MV_TEST tests R8BTO_MV.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    14 March 2009
%
%  Author:
%
%    John Burkardt
%
  l = 3;
  m = 2;
  n = m * l;
%
%  I don't know how MATLAB wants me to enter a triply-dimensioned array.
%  I've tried several different approaches with no success, so 
%  here's the clodhopper's method:
%
  a = zeros ( m, m, 2*l-1 );
  
  a(1,1,1) = 1.0;
  a(2,1,1) = 5.0; 
  a(1,2,1) = 2.0;
  a(2,2,1) = 5.0;
  a(1,1,2) = 3.0;
  a(2,1,2) = 6.0; 
  a(1,2,2) = 4.0;
  a(2,2,2) = 6.0;
  a(1,1,3) = 5.0;
  a(2,1,3) = 7.0; 
  a(1,2,3) = 6.0;
  a(2,2,3) = 7.0;
  a(1,1,4) = 7.0;
  a(2,1,4) = 8.0; 
  a(1,2,4) = 8.0;
  a(2,2,4) = 8.0;
  a(1,1,5) = 9.0;
  a(2,1,5) = 9.0; 
  a(1,2,5) = 0.0;
  a(2,2,5) = 9.0;
  
  fprintf ( 1, '\n' );
  fprintf ( 1, 'R8BTO_MV_TEST\n' );
  fprintf ( 1, '  R8BTO_MV computes A * x.\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 A:' );

  xmat = r8ge_indicator ( m, l );
   
  x = reshape ( xmat, n, 1 );

  r8vec_print ( n, x, '  The vector x:' );

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

  r8vec_print ( n, b, '  The product A*x:' );

  return
end
