function zipf_probability_test ( )

%*****************************************************************************80
%
%% ZIPF_PROBABILITY_TEST tests ZIPF_PROBABILITY.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    21 February 2016
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'ZIPF_PROBABILITY_TEST\n' );
  fprintf ( 1, '  ZIPF_PROBABILITY sets up a probablity vector X of N+2\n' );
  fprintf ( 1, '  elements containing in X[1:N] the probabilities of\n' );
  fprintf ( 1, '  outcomes 1 through Nin a Zipf distribution with\n' );
  fprintf ( 1, '  parameter P.\n' );

  n = 5;
  p = 1.0;
  x = zipf_probability ( n, p );
  r8vec_print ( n + 2, x, '  X for N = 5, P = 1.0' );

  n = 5;
  p = 2.0;
  x = zipf_probability ( n, p );
  r8vec_print ( n + 2, x, '  X for N = 5, P = 2.0' );

  n = 10;
  p = 2.0;
  x = zipf_probability ( n, p );
  r8vec_print ( n + 2, x, '  X for N = 10, P = 2.0' );
%
%  Terminate.
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'ZIPF_PROBABILITY_TEST\n' );
  fprintf ( 1, '  Normal end of execution.\n' );

  return
end

