function zeta_m1_test ( )

%*****************************************************************************80
%
%% ZETA_M1_TEST tests ZETA_M1.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    16 January 2017
%
%  Author:
%
%    John Burkardt
%
  tol = 1.0E-10;
  
  fprintf ( 1, '\n' );
  fprintf ( 1, 'ZETA_M1_TEST\n' );
  fprintf ( 1, '  ZETA_M1 evaluates the Riemann Zeta Minus One function.\n' );
  fprintf ( 1, '  Requested tolerance = %g\n', tol );
  fprintf ( 1, '\n' );
  fprintf ( 1, '       P              Zeta_M1               Zeta_M1\n' );
  fprintf ( 1, '                        exact              computed\n' );
  fprintf ( 1, '\n' );

  n_data = 0;

  while ( true )

    [ n_data, p, z1 ] = zeta_m1_values ( n_data );

    if ( n_data == 0 )
      break
    end

    z2 = zeta_m1 ( p, tol );

    fprintf ( 1, '  %8.2f  %20e  %20e\n', p, z1, z2 );

  end

  return
end
