function value = steinerberger_integral01 ( n )

%*****************************************************************************80
%
%% steinerberger_integral01() returns the Steinerberger integral from 0 to 1.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    16 June 2023
%
%  Author:
%
%    John Burkardt
%
%  Reference:
%
%    John D Cook,
%    Pushing numerical integration software to its limits,
%    https://www.johndcook.com/blog/2023/06/12/stressing-numerical-integration/
%    Posted 12 June 2023.
%
%    John D Cook,
%    Plotting a function with lots of local minima,
%    https://www.johndcook.com/blog/2023/06/12/lots-of-local-minima/
%    Posted 12 June 2023.
%
%    Stefan Steinerberger,
%    A amusing sequence of functions,
%    Mathematics Magazine,
%    Volume 91, Number 4, October 2018, pages 262-266.
%
%  Input:
%
%    integer N: the index of the function.
%
%  Output:
%
%    real VALUE: the value of the integral from 0 to 1.
%
  h = harmonic_number ( n );
  value = 2.0 * h / pi;

  return
end

