function u_mass_matrix_test ( )

%*****************************************************************************80
%
%% u_mass_matrix_test() tests u_mass_matrix().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    15 July 2015
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'u_mass_matrix_test():\n' );
  fprintf ( 1, '  u_mass_matrix() computes the mass matrix for the\n' );
  fprintf ( 1, '  Chebyshev U polynomials U(i,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 = 3;

  a = u_mass_matrix ( n );

  r8mat_print ( n + 1, n + 1, a, '  U mass matrix:' );

  return
end
