function plums_test ( )

%*****************************************************************************80
%
%% plums_test() tests plums().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    10 January 2022
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'plums_test():\n' );
  fprintf ( 1, '  plums() estimates the expected distance D from the surface\n' );
  fprintf ( 1, '  to the nearest of N plums randomly placed inside a unit\n' );
  fprintf ( 1, '  spherical pudding.\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  N  D(Estimate)  D(Exact)\n' );
  fprintf ( 1, '\n' );

  for n = 5 : 10
    s = plums ( n );
    s_exact = 1.0 / ( 3 * n + 1 );
    fprintf ( 1, '  %2d  %g  %g\n', n, s, s_exact );
  end

  return
end

