function fisher_pdf_test ( )

%*****************************************************************************80
%
%% FISHER_PDF_TEST tests FISHER_PDF.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    10 April 2016
%
%  Author:
%
%    John Burkardt
%
  clear

  n = 10;
  test_num = 3;

  fprintf ( 1, '\n' );
  fprintf ( 1, 'FISHER_PDF_TEST\n' );
  fprintf ( 1, '  FISHER_PDF evaluates the Fisher PDF.\n' );

  for test = 1 : test_num

    if ( test == 1 )
      kappa = 0.0;
      mu = [ 1.0, 0.0, 0.0 ];
    elseif ( test == 2 )
      kappa = 0.5;
      mu = [ 1.0, 0.0, 0.0 ];
    elseif ( test == 3 ) 
      kappa = 10.0;
      mu = [ 1.0, 0.0, 0.0 ];
    end

    fprintf ( 1, '\n' );
    fprintf ( 1, '  PDF parameters:\n' );
    fprintf ( 1, '    Concentration parameter KAPPA =      %f\n', kappa );
    fprintf ( 1, '    Direction MU(1:3) = %f  %f  %f\n', mu(1:3) );

    fprintf ( 1, '\n' );
    fprintf ( 1, '      X                         PDF\n' );
    fprintf ( 1, '\n' );

    seed = 123456789;

    for j = 1 : n

      [ x, seed ] = fisher_sample ( kappa, mu, n, seed );

      pdf = fisher_pdf ( x(1:3), kappa, mu );

      fprintf ( 1, '  %8f  %8f  %8f    %14g\n', x(1:3), pdf );

    end

  end

  return
end
