function i4_binomial_pdf_test ( )

%*****************************************************************************80
%
%% I4_BINOMIAL_PDF_TEST calls I4_BINOMIAL_PDF.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license. 
%
%  Modified:
%
%    24 January 2018
%
%  Author:
%
%    John Burkardt
%
  initialize ( );

  fprintf ( 1, '\n' );
  fprintf ( 1, 'I4_BINOMIAL_PDF_TEST\n' );
  fprintf ( 1, '  I4_BINOMIAL_PDF reports\n' );
  fprintf ( 1, '  PROB, the probability that\n' );
  fprintf ( 1, '  N trials, with\n' );
  fprintf ( 1, '  P probability of success result in\n' );
  fprintf ( 1, '  K successes.\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '   N         P   K        PROB\n' );
  fprintf ( 1, '\n' );

  n = 5;
  p = 0.25;
 
  for k = 0 : n
    prob = i4_binomial_pdf ( n, p, k );
    fprintf ( 1, '  %2d  %8f  %2d  %14.6g\n', n, p, k, prob );
  end

  return
end