function vertex_contains_ijk_test ( )

%*****************************************************************************80
%
%% vertex_contains_ijk_test tests vertex_contains_ijk().
%
%  Discussion:
%
%    5           x x x x-x x
%                 /|  / /
%    4         x-x x x x x
%             /    |/ /
%    3       x-x x x x x
%              |    /
%    2     x x x x x x
%             /    |
%    1   x-x-x x x x
%       /         /
%    0 x-x-x-x-x-x
%
%      0 1 2 3 4 5
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    03 November 2020
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'vertex_contains_ijk_test\n' );
  fprintf ( 1, '  MATLAB/Octave version %s\n', version ( ) );
  fprintf ( 1, '  vertex_contains_ijk() reports whether an IJK triangle is\n' );
  fprintf ( 1, '  inside a polygonal region with given vertex list.\n' );

  imin = 0;
  imax = 5;
  jmin = 0;
  jmax = 5;

  v = vertex_gort ( );
  name = '  Gort vertex list';
  vertex_print ( v, name );

  fprintf ( 1, '\n' );
  fprintf ( 1, '   i   j   k  inside?\n' );
  fprintf ( 1, '\n' );

  for j = jmin + 1 : jmax
    for i = imin + 1: imax
      for k = 0 : 1
        inside = vertex_contains_ijk ( v, [ i, j, k ] );
        tf = boolean_to_string ( inside );
        fprintf ( 1, '  %2d  %2d  %2d  %s\n', i, j, k, tf );
      end
    end
  end
%
%  Terminate.
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'vertex_contains_ijk_test\n' );
  fprintf ( 1, '  Normal end of execution.\n' );

  return
end
