function r8_error_f_test ( )

%*****************************************************************************80
%
%% R8_ERROR_F_TEST tests R8_ERROR_F.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    02 April 2016
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'R8_ERROR_F_TEST\n' );
  fprintf ( 1, '  R8_ERROR_F evaluates the error function ERF.\n' );

  fprintf ( 1, '\n' );
  fprintf ( 1, 'X   -> Y = R8_ERROR_F(X) -> Z = R8_ERROR_F_INVERSE(Y)\n' );
  fprintf ( 1, '\n' );

  seed = 123456789;

  for i = 1 : 20

    [ x, seed ] = normal_01_sample ( seed );
    y = r8_error_f ( x );
    z = r8_error_f_inverse ( y );
    fprintf ( 1, '  %14g  %14g  %14g\n', x, y, z );

  end

  return
end
