function r8utt_mm_test ( )

%*****************************************************************************80
%
%% R8UTT_MM_TEST tests R8UTT_MM.
%
%  Licensing:
%
%    This code is distributed under the MIT license. 
%
%  Modified:
%
%    31 October 2015
%
%  Author:
%
%    John Burkardt
%
  n = 5;
 
  fprintf ( 1, '\n' );
  fprintf ( 1, 'R8UTT_MM_TEST\n' );
  fprintf ( 1, '  R8UTT_MM computes C = A * B for R8UTT matrices.\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  Matrix order N = %d\n', n );

  seed = 123456789;
  [ a, seed ] = r8utt_random ( n, seed );
  r8utt_print ( n, a, '  Factor A:' );
  [ b, seed ] = r8utt_random ( n, seed );
  r8utt_print ( n, b, '  Factor B:' );
  c = r8utt_mm ( n, a, b );
  r8utt_print ( n, c, '  The product C = A * B' );

  a_ge = r8utt_to_r8ge ( n, a );
  b_ge = r8utt_to_r8ge ( n, b );
  c_ge = a_ge * b_ge;
  r8ge_print ( n, n, c_ge, '  The R8GE product C:' )

  return
end
