function uu_product_integral_test ( )

%*****************************************************************************80
%
%% uu_product_integral_test() tests uu_product_integral().
%
%  Discussion:
%
%    This process should match the UU_MASS_MATRIX computation.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    13 July 2015
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'UU_PRODUCT_INTEGRAL_TEST:\n' );
  fprintf ( 1, '  UU_PRODUCT_INTEGRAL computes the product integral\n' );
  fprintf ( 1, '  of a pair of Chebyshev U polynomials U(i,x) and U(j,x).\n' );
  fprintf ( 1, '  A(I,J) = integral ( -1 <=x <= +1 ) U(i,x) U(j,x) sqrt ( 1 - x^2 ) dx\n' );
  fprintf ( 1, '  0    if i is not equal to j;\n' );
  fprintf ( 1, '  pi/2 if i = j\n' );

  n = 4;

  a = zeros ( n + 1, n + 1 );

  for i = 0 : n
    for j = 0 : n
      a(i+1,j+1) = uu_product_integral ( i, j );
    end
  end

  r8mat_print ( n + 1, n + 1, a, '  U(i,x)*U(j,x) integral matrix:' );

  return
end
