function value = r8_epsilon_compute ( )

%*****************************************************************************80
%
%% r8_epsilon_compute() returns the R8 roundoff unit.
%
%  Discussion:
%
%    The roundoff unit is a number R which is a power of 2 with the 
%    property that, to the precision of the computer's arithmetic,
%      1 < 1 + R
%    but 
%      1 = ( 1 + R / 2 )
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    11 August 2010
%
%  Author:
%
%    John Burkardt
%
%  Output:
%
%    real VALUE, the unit roundoff.
%
  one = 1;

  value = one;
  temp = value / 2.0;
  test = r8_add ( one, temp );

  while ( one < test )  
    value = temp;
    temp = value / 2.0;
    test = r8_add ( one, temp );
  end

  return
end
