function disk01_integrands_test01 ( )

%*****************************************************************************80
%
%% disk01_integrands_test01(): integral of X^EX in the unit disk.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    19 April 2019
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'disk01_integrands_test01():\n' );
  fprintf ( 1, '  Use a simple Monte Carlo approach to estimate the\n' );
  fprintf ( 1, '  integral of X^EX inside the disk of radius 1\n' );
  fprintf ( 1, '  centered at the origin.\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '           N   EX      Exact Approximate   Error\n' );

  for ex = 2 : 2 : 6
    fprintf ( 1, '\n' );
    p01_parameters ( ex );
    n = 1;
    for n_log2 = 0 : 20
      xy = disk_unit_sample ( n );
      f = p01_f ( n, xy );
      q = disk_unit_area ( ) * sum ( f ) / n;
      t = p01_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
