function hermite_polynomial_test13 ( p, e )

%*****************************************************************************80
%
%% hermite_polynomial_test13() tests hf_power_product().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    25 February 2012
%
%  Author:
%
%    John Burkardt
%
%  Input:
%
%    integer P, the maximum degree of the polynomial factors.
%
%    integer E, the exponent of X in the integrand.
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'hermite_polynomial_test13():\n' );
  fprintf ( 1, '  Compute a Hermite function power product table.\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  Tij = integral ( -oo < X < +oo ) X^E Hf(I,X) Hf(J,X) dx\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  where Hf(I,X) = Hermite function of "degree" I.\n' );

  fprintf ( 1, '\n' );
  fprintf ( 1, '  Maximum degree P = %d\n', p );
  fprintf ( 1, '  Exponent of X = %d\n', e );

  table = hf_power_product ( p, e );

  r8mat_print ( p + 1, p + 1, table, '  Power product table:' );

  return
end
