function p = disk_positive_sample ( )

%*****************************************************************************80
%
%% disk_positive_sample() samples the unit positive disk.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    11 August 2023
%
%  Author:
%
%    John Burkardt
%
%  Output:
%
%    real p(1,2): a point selected uniformly at random from
%    the positive disk of radius 1 and center (0,0).
%
  theta = 2.0 * pi * rand ( 1, 1 );
  r = sqrt ( rand ( 1, 1 ) );
  x = abs ( r * cos ( theta ) );
  y = abs ( r * sin ( theta ) );
  p = [ x, y ];

  return
end

