function uniform_in_annulus_sector_test ( )

%*****************************************************************************80
%
%% uniform_in_annulus_sector_test() tests uniform_in_annulus_sector().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    13 April 2022
%
%  Author:
%
%    John Burkardt
%
  n = 1000;
  pc = [ 10.0, 5.0 ];
  r1 = 1.0;
  r2 = 3.0;
  theta1 = 0.0;
  theta2 = pi / 2.0;

  fprintf ( 1, '\n' );
  fprintf ( 1, 'uniform_in_annulus_sector_test():\n' );
  fprintf ( 1, '  uniform_in_annulus_sector() generates uniform \n' );
  fprintf ( 1, '  points in an annular sector by mapping.\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  Number of points N =          %12d\n', n );
  fprintf ( 1, '  Center PC(1:2) =              %12f  %12f\n', pc(1:2) );
  fprintf ( 1, '  Inner radius is R1 =          %12f\n', r1 );
  fprintf ( 1, '  Outer radius is R2 =          %12f\n', r2 );
  fprintf ( 1, '  THETA1 =                      %12f\n', theta1 );
  fprintf ( 1, '  THETA2 =                      %12f\n', theta2 );
 
  p = uniform_in_annulus_sector ( n, pc, r1, r2, theta1, theta2 );
%
%  Make a plot.
%
  clf ( );
  plot ( p(:,1), p(:,2), 'b.' );
  grid ( 'on' );
  axis ( 'equal' );
  xlabel ( '<--- X --->' );
  ylabel ( '<--- Y --->' );
  title ( 'uniform\_in\_annulus\_sector()' );
  filename = 'uniform_in_annulus_sector.png';
  print ( '-dpng', filename );
  fprintf ( 1, '  Graphics saved as "%s".\n', filename );

  return
end
