function x = p16_root ( i )

%*****************************************************************************80
%
%% p16_root returns a root for problem 16.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    10 October 2020
%
%  Author:
%
%    John Burkardt
%
%  Input:
%
%    integer I, the index of the requested root.
%
%  Output:
%
%    real X, the value of the root.
%
  if ( i == 1 )
    x = 22.65657866956775;
  else
    fprintf ( 1, '\n' );
    fprintf ( 1, 'P16_root - Fatal error!\n' );
    fprintf ( 1, '  Illegal root index = %d\n', i );
    error ( 'P16_root - Fatal error!' );
  end

  return
end
