function pythag_triple_ijk_test ( )

%*****************************************************************************80
%
%% pythag_triple_ijk_test() tests pythag_triple_ijk().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    17 October 2020
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'pythag_triple_ijk_test():\n' );
  fprintf ( 1, '  pythag_triple_ijk() computes a primitive Pythagorean triple\n' );
  fprintf ( 1, '  from i, j, k "coordinates".\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '   I   J   K   A   B   C  A^2+B^2   C^2\n' );
  fprintf ( 1, '\n' );

  for i = 0 : 3
    for j = 0 : 3
      for k = 0 : 3
        [ a, b, c ] = pythag_triple_ijk ( i, j, k );
        d = a^2 + b^2;
        e = c^2;
        fprintf ( 1, '%4d  %4d  %4d  %6d  %6d  %6d  %11d  %11d\n', ...
          i, j, k, a, b, c, d, e );
      end
    end
  end

  return
end
