function elliptic_inc_pik_test ( )

%*****************************************************************************80
%
%% ELLIPTIC_INC_PIK_TEST tests ELLIPTIC_INC_PIK.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    22 June 2018
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'ELLIPTIC_INC_PIK_TEST:\n' );
  fprintf ( 1, '  ELLIPTIC_INC_PIK returns values of\n' );
  fprintf ( 1, '  the incomplete elliptic integral of the\n' );
  fprintf ( 1, '  third kind, with parameters PHI, N, K.\n' );
  fprintf ( 1, '  Compare with MATLAB ellipticPi(n,phi,k^2)\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '          Phi             N               K                ')
  fprintf ( 1, 'ellipticPi(n,phi,k^2)  elliptic_inc_pik(phi,n,k)\n' );
  fprintf ( 1, '\n' );

  seed = 123456789;

  for i = 1 : 20

    [ phi, seed ] = r8_uniform_ab ( 0.0, 0.5 * pi, seed );
    nbound = ( 1.0 / sin ( phi ) ) ^ 2;
    [ n, seed ] = r8_uniform_ab ( - nbound, + nbound, seed );
    kbound = abs ( 1.0 / sin ( phi ) );
    [ k, seed ] = r8_uniform_ab ( - kbound, + kbound, seed );
    fx1 = ellipticPi ( n, phi, k^2 );
    fx2 = elliptic_inc_pik ( phi, n, k );

    fprintf ( 1, '  %14.6f  %14.6f  %14.6f  %24.16g  %24.16g\n', phi, n, k, fx1, fx2 );

  end

  return
end
