function maple_area_estimate_qmc_test ( b, p_m, p_n )

%*****************************************************************************80
%
%% maple_area_estimate_qmc_test() tests maple_area_estimate_qmc().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    20 August 2016
%
%  Author:
%
%    John Burkardt
%
%  Input:
%
%    integer B(*,2), a list of pixels that form the boundary.
%
%    integer P_M, P_N, the number of rows and columns of pixels.
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'maple_area_estimate_qmc_test():\n' );
  fprintf ( 1, '  maple_area_estimate_qmc() is given a polygon defining the\n' );
  fprintf ( 1, '  boundary of the maple leaf, and uses SAMPLE_NUM Quasi Monte Carlo\n' );
  fprintf ( 1, '  points to estimate the relative area.\n' );
  fprintf ( 1, '  Hayes got a relative area of 0.4185 using a 1024x1024 grid.\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '   I  SAMPLE_NUM        Estimate\n' );
  fprintf ( 1, '\n' );

  n = 1;

  for log2n = 0 : 10

    sample_num = n * n;

    estimate = maple_area_estimate_qmc ( b, p_m, p_n, sample_num );

    fprintf ( 1, '  %2d    %8d  %14.6f\n', log2n, sample_num, estimate );
    n = 2 * n;

  end

  return
end

