function v_mass_matrix_test ( )

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

  n = 3;

  a = t_mass_matrix ( n );

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

  return
end
