function uu_product_test ( )

%*****************************************************************************80
%
%% uu_product_test() tests uu_product().
%
%  Licensing:
%
%    This code is distributed under the MIT license. 
%
%  Modified:
%
%    13 July 2015
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'UU_PRODUCT_TEST:\n' );
  fprintf ( 1, '  UU_PRODUCT(I,J;X) = U(I,X) * U(J,X)\n' );

  fprintf ( 1, '\n' );
  fprintf ( 1, '   I   J      X               UI              UJ              UI*UJ       UU_PRODUCT\n' ); 
  fprintf ( 1, '\n' ); 
  for test = 1 : 10
    x = - 1.0 + 2.0 * rand ( 1, 1 );
    i = randi ( [ 0, 6 ], 1 );
    ui = u_polynomial_value ( i, x );
    j = randi ( [ -1, 4 ], 1 );
    uj = u_polynomial_value ( j, x );
    uiuj = uu_product ( i, j, x );
    fprintf ( 1, '  %2d  %2d  %14.6g  %14.6g  %14.6g  %14.6g  %14.6g\n', ...
      i, j, x, ui, uj, ui * uj, uiuj );
  end

  return
end
