function triangulation_test21 ( )

%*****************************************************************************80
%
%% TEST21 tests TRIANGULATION_ORDER3_PRINT.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    24 August 2006
%
%  Author:
%
%    John Burkardt
%
  node_num = 9;
  triangle_num = 12;

  node_xy(1:2,1:node_num) = [ ...
       0.0, 0.0; ...
       0.0, 1.0; ...
       0.2, 0.5; ...
       0.3, 0.6; ...
       0.4, 0.5; ...
       0.6, 0.4; ...
       0.6, 0.5; ...
       1.0, 0.0; ...
       1.0, 1.0 ]';
  triangle_node(1:3,1:triangle_num) = [ ...
       2, 1, 3; ...
       3, 1, 6; ...
       2, 3, 4; ...
       4, 3, 5; ...
       7, 4, 5; ...
       5, 3, 6; ...
       7, 5, 6; ...
       9, 4, 7; ...
       6, 1, 8; ...
       7, 6, 8; ...
       7, 8, 9; ...
       2, 4, 9 ]';
  triangle_neighbor(1:3,1:triangle_num) = [ ...
       -28,   2,  3; ...
         1,   9,  6; ...
         1,   4, 12; ...
         3,   6,  5; ...
         8,   4,  7; ...
         4,   2,  7; ...
         5,   6, 10; ...
        12,   5, 11; ...
         2, -34, 10; ...
         7,   9, 11; ...
        10, -38,  8; ...
         3,   8, -3 ]';

  fprintf ( 1, '\n' );
  fprintf ( 1, 'TEST21\n' );
  fprintf ( 1, '  TRIANGULATION_ORDER3_PRINT prints out a triangulation.\n' );

  triangulation_order3_print ( node_num, triangle_num, node_xy, ...
    triangle_node, triangle_neighbor );

  return
end
