function disk01_integrands_test02 ( )

%*****************************************************************************80
%
%% disk01_integrands_test02(): integral of R^EX in the unit disk.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    11 October 2022
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'disk01_integrands_test02():\n' );
  fprintf ( 1, '  Use a simple Monte Carlo approach to estimate the\n' );
  fprintf ( 1, '  integral of R^EX over the disk of radius 1\n' );
  fprintf ( 1, '  centered at the origin.\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '           N   E       Exact Approximate   Error\n' );

  for ex = 1 : 2 : 5
    fprintf ( 1, '\n' );
    p02_parameters ( ex );
    n = 1;
    for n_log2 = 0 : 20
      xy = disk_unit_sample ( n );
      f = p02_f ( n, xy );
      q = disk_unit_area ( ) * sum ( f ) / n;
      t = p02_exact ( );
      fprintf ( 1, '  %10d  %2d  %10.4f  %10.4f  %10.4e\n', n, ex, t, q, abs ( t - q ) );
      n = n * 2;
    end
  end

  return
end
