function t_mass_matrix_test ( )

%*****************************************************************************80
%
%% t_mass_matrix_test() tests t_mass_matrix().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    15 July 2015
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 't_mass_matrix_test():\n' );
  fprintf ( 1, '  t_mass_matrix() computes the mass matrix for the\n' );
  fprintf ( 1, '  Chebyshev T polynomials T(i,x).\n' );
  fprintf ( 1, '  A(I,J) = integral ( -1 <=x <= +1 ) T(i,x) T(j,x) / sqrt ( 1 - x^2 ) dx\n' );
  fprintf ( 1, '  0    if i is not equal to j;\n' );
  fprintf ( 1, '  pi   if i = j = 0;\n' );
  fprintf ( 1, '  pi/2 if i = j =/= 0.\n' );

  n = 3;

  a = t_mass_matrix ( n );

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

  return
end
