function r83s_mtv_test ( )

%*****************************************************************************80
%
%% R83S_MTV_TEST tests R83S_MTV.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    05 September 2015
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'R83S_MTV_TEST\n' );
  fprintf ( 1, '  R83S_MV computes b=A''*x, where A is an R83S matrix.\n' );
  fprintf ( 1, '  We check three cases, M<N, M=N, M>N.\n' );

  for i = 1 : 3

    if ( i == 1 )
      m = 3;
      n = 5;
    elseif ( i == 2 )
      m = 5;
      n = 5;
    elseif ( i == 3 )
      m = 5;
      n = 3;
    end

    seed = 123456789;
    [ a_83s, seed ] = r83s_random ( m, n, seed );
    x = r8vec_indicator1 ( m );
    ax_83s = r83s_mtv ( m, n, a_83s, x );
    a_ge = r83s_to_r8ge ( m, n, a_83s );
    ax_ge = r8ge_mtv ( m, n, a_ge, x );
    r8vec2_print ( n, ax_83s, ax_ge, '  Product comparison:' );

  end

  return
end
