function disk_integrals_test01 ( n )

%*****************************************************************************80
%
%% disk_integrals_test01() compares exact and Monte Carlo integral estimates.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    04 November 2002
%
%  Author:
%
%    John Burkardt
%
%  Input:
%
%    integer n: the number of sample points to use.
%
  m = 2;
  test_num = 20;

  fprintf ( 1, '\n' );
  fprintf ( 1, 'disk_integrals_test01():\n' );
  fprintf ( 1, '  Estimate monomial integrals using Monte Carlo\n' );
  fprintf ( 1, '  over the interior of the unit disk in 2D.\n' );
%
%  Get sample points.
%
  r = 2.0;
  x = disk_sample ( r, n );

  fprintf ( 1, '\n' );
  fprintf ( 1, '  Disk radius is %g\n', r );
  fprintf ( 1, '  Number of sample points is %d\n', n );
%
%  Randomly choose X,Y exponents between 0 and 8.
%
  fprintf ( 1, '\n' );
  fprintf ( 1, '  If any exponent is odd, the integral is zero.\n' );
  fprintf ( 1, '  We will restrict this test to randomly chosen even exponents.\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  Ex  Ey     MC-Estimate           Exact      Error\n' );
  fprintf ( 1, '\n' );

  for e1 = 0 : 2 : 8

    for e2 = 0 : 2 : 8

      e = [ e1, e2 ];

      value = monomial_value ( m, n, e, x );

      result = disk_area ( r ) * sum ( value(1:n) ) / n;
      exact = disk_monomial_integral ( r, e );
      error = abs ( result - exact );

      fprintf ( 1, '  %2d  %2d  %14.6g  %14.6g  %10.2e\n', ...
        e(1:m), result, exact, error );

    end
  end

  return
end
