function truncated_normal_b_pdf_test ( )

%*****************************************************************************80
%
%% TRUNCATED_NORMAL_B_PDF_TEST tests TRUNCATED_NORMAL_B_PDF;
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    14 September 2013
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'TRUNCATED_NORMAL_B_PDF_TEST\n' );
  fprintf ( 1, '  TRUNCATED_NORMAL_B_PDF evaluates the Upper Truncated Normal PDF.\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  The "parent" normal distribution has\n' );
  fprintf ( 1, '    mean = mu\n' );
  fprintf ( 1, '    standard deviation = sigma\n' );
  fprintf ( 1, '  The parent distribution is truncated to\n' );
  fprintf ( 1, '  the interval (-oo,b]\n' );

  fprintf ( 1, '\n' );
  fprintf ( 1, '                                                  Stored         Computed\n' );
  fprintf ( 1, '       X        Mu         S         B             PDF             PDF\n' );
  fprintf ( 1, '\n');

  n_data = 0;

  while ( true )

    [ n_data, mu, sigma, b, x, pdf1 ] ...
      = truncated_normal_b_pdf_values ( n_data );

    if ( n_data == 0 )
      break
    end

    pdf2 = truncated_normal_b_pdf ( x, mu, sigma, b );

    fprintf( 1, '  %8.1f  %8.1f  %8.1f  %8.1f  %14g  %14g\n', ...
      x, mu, sigma, b, pdf1, pdf2 );

  end

  return
end
